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Abstract 



The central topic of this thesis is two dimensional Quantum Gravity and its properties. 
The term Quantum Gravity itself is ambiguous as there are many proposals for its correct 
formulation and none of them have been verified experimentally. In this thesis we consider 
a number of closely related approaches to two dimensional quantum gravity that share 
the property that they may be formulated in terms of random graphs. In one such 
approach known as Causal Dynamical Triangulations, numerical computations suggest 
an interesting phenomenon in which the effective spacetime dimension is reduced in the 
UV. In this thesis we first address whether such a dynamical reduction in the number of 
dimensions may be understood in a simplified model. We introduce a continuum limit 
where this simplified model exhibits a reduction in the effective dimension of spacetime 
in the UV, in addition to having rich cross-over behaviour. 

In the second part of this thesis we consider an approach closely related to causal 
dynamical triangulation; namely dynamical triangulation. Although this theory is less 
well-behaved than causal dynamical triangulation, it is known how to couple it to matter, 
therefore allowing for potentially multiple boundary states to appear in the theory. We 
address the conjecture of Seiberg and Shih which states that all these boundary states 
are degenerate and may be constructed from a single principal boundary state. By use 
of the random graph formulation of the theory we compute the higher genus amplitudes 
with a single boundary and find that they violate the Seiberg-Shih conjecture. Finally we 
discuss whether this result prevents the replacement of boundary states by local operators 
as proposed by Seiberg. 



i 



ii 



Acknowledgments 



I would like to thank my supervisor, John Wheater, for the significant time he has spent 
helping me and the numerous interesting projects he suggested. This thesis would not 
exist without his help. I would also like to give a very large thank you to Stefan Zohren 
who has provided significant encouragement and motivation in the final year of my DPhil. 
Both Stefan and John have been of invaluable help during the writing of this thesis, both 
suggesting numerous improvements and raising interesting questions. 

I would like to acknowledge the support of STFC studentship PfPA/S/S/2006/04507 
in providing me with the funding to undertake this research. 

I would like to thank the senior tutors at Christ Church college, Axel Kuhn, Derek 
Stacey and Alan Merchant for providing me with the opportunity to teach and partake 
in admissions. This was a life-saver; especially when funds were short! My enjoyment of 
my time at Oxford was also greatly enhanced by my office mates, Sesh, Shaun, Seung Joo 
and Tom, who I thank for providing constant interesting conversation, whether related to 
physics or not. 

My family have also been a source of encouragement throughout my DPhil and have 
always been there to offer support. Finally, I want to thank Lauren for being there 
throughout it all and her endless willingness to lend an ear. It is to my family and Lauren 
to which I dedicate this thesis. 



iii 



iv 



Contents 



Abstract Q 

Acknowledgments 

Contents [hj] 

List of Figures Ivjj] 

1 Introduction |l 

1.1 Thesis Outline |4 

2 Discrete Approaches to 2D Quantum Gravity 3 

2.1 Dynamical Triangulation ^ 

2.1.1 Random Graphs and Matrix Models n| 

2.1.2 Solving the Matrix Model: Loop equations 

2.1.3 The Continuum Limit 

2.1.4 Combinatorial Interpretation of the Loop Equation 

2.2 Causal Dynamical Triangulation 

3 The Dimension of Spacetime [23 

3.1 The Hausdorff Dimension dn [24 

3.1.1 Computing d H for DT and CDT 2J 

3.2 The Spectral Dimension d s [28 

3.2.1 Dimensional Reduction [30 

3.2.2 The Spectral Dimension in CDT [31 



v 



vi CONTENTS 



4 A Toy Model of Dimensional Reduction 

4.1 Combs and Walks 

4.1.1 Definitions 

4.1.2 Random Walks 

4.1.3 Two point functions 

4.1.4 Spectral dimension and the continuum limit 

4.2 A simple comb 

4.3 Combs with Power Law Measures 

4.4 Multiple Scales 

4.5 Generic Distributions 

4.6 Discussion 

5 DT with Matter, String Theory and Branes 

5.1 Adding Matter using Matrix Models 

5.1.1 The Continuum Limit and Multicritical Points 

5.2 Conformal Field Theory and Minimal Models 

5.3 Boundary States of Minimal Models 

5.4 String Theory 

5.4.1 Liouville Theory and Non-Critical String Theory 

5.4.2 Wave-functions and Correlation functions in Liouville Theory 

5.4.3 Minimal String Theory 

5.4.4 FZZT Branes and the Seiberg-Shih Relation 

5.4.5 Cylinder Amplitudes and the Seiberg-Shih Relation 

6 FZZT Branes in the (3, 4) Minimal String at Higher Genus 

6.1 Testing the Seiberg-Shih relations using a Matrix Model 

6.1.1 Loop equations for SX Matrix Model 

6.1.2 Disc, Cylinders and Disc-with-Handles 

6.1.3 Dual Branes 

6.2 The Seiberg-Shih Relation and Local Operators 

6.2.1 Extracting Local Operators From Loops 

6.3 Discussion 



101 



101 



105 



107 



7 Summary 



111 



CONTENTS vii 



A Mellin Transforms and Asymptotics 1115 

A.l Standard generating functions 115 

A. 2 Bump functions 116 

A. 3 Asymptotic Series and Dirichlet Series Il IT 



B Liouville Cylinder Amplitudes 

C Loop equations for Disc with Handle 



121 



125 



Bibliography 



127 



CONTENTS 



List of Figures 



2.7 A typical geometry in CDT 



3.1 By drawing a baby universe in the lattice we see how this gives rise to a 
fractal structure [45] 

3.2 A causal triangulation drawn in the plane. Note we have added an extra 
vertex with links to all vertices in the inner most spatial slice [16] 

3.3 Casual triangulations have a bijection with trees. In this case the tree 
associated to the triangulation is shown by thick lines. [16] 



4.1 A comb 

4.2 An illustration of the strip E in which T>^{s\ Wi) satisfies property (1) given 
given in the text. The poles of V^(s;Wi) are denoted by crosses and the 
horizontal strips in which the growth condition holds are indicated by the 
dark grey regions 

ix 



19 



2.1 A triangulation (grey) together with its dual graph (black) [10 

2.2 Move 1 relates an unrestricted triangulation (left), with a mark denoted by 
the white circle, to another unrestricted triangulation with less triangles. 

2.3 Move 2 relates an unrestricted triangulation (left), with a mark denoted by 
the white circle, to two disconnected unrestricted triangulations 

2.4 A causal triangulation 

2.5 A triangulation of a baby universe. If we successively remove triangles, 
then the final link between the baby universe and the parent will be a 
double link 

2.6 A typical geometry in dynamical triangulation exhibiting the lack of a 
smooth continuum limit [55] 



21 



21 



22 



27 



32 



32 



37 



54 



X 



LIST OF FIGURES 



5.1 The worldsheet corresponding to the interaction of two strings may be 
transformed via a conformal transformation to a worldsheet of equivalent 
topology but with all asymptotic string states transformed into local oper- 
ators 72 



CHAPTER 1 



Introduction 



Perhaps the most ambitious question that one may try to answer is "What is reality?", 
where by reality we mean the sensation of being alive, of experience. A priori there are 
many possibilities; a subscriber to idealism would believe one's experiences are entirely 
generated by one's mind. If, on the other hand, one subscribes to materialism, then what 
constitutes reality is a representation, generated by one's brain, of an independent physical 
universe. The quest to understand reality then becomes a matter of understanding the 
true nature of this physical universe. The 20th century saw many advancements in our 
understanding of the universe in this respect. The development of quantum mechanics 
fundamentally altered what is considered the true nature of the matter and forces, while 
the introduction of general relativity, which describes the dynamics of spacetime and how 
this manifests itself as gravity, banished the perception of space and time as being a static 
background stage on which dynamical processes took place. 

However, the revolution in our view of the universe caused by quantum mechanics and 
general relativity was only half completed. General relativity assumes that all matter re- 
siding in spacetime is classical. However, we know such matter is fundamentally quantum 
mechanical and so we must, at the very least, construct a theory of gravity in which the 
matter is allowed to assume its true quantum mechanical form. 

There have been attempts to pursue the most conservative line of attack; keep the 
dynamics of spacetime as they are in general relativity, or at least classical, but modify 
how matter couples to gravity in order to allow for quantum mechanical rather than 
classical sources. However, there are a number of compelling arguments to suggest this 
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does not produce a consistent theory [35]. A more natural possibility is that gravity itself 
is quantum mechanical in nature. This avoids the aforementioned problems of coupling 
classical and quantum mechanical systems in addition to the philosophical appeal of 
all phenomenon in the universe being fundamentally quantum. Such a theory and the 
attempts to construct it are known as quantum gravity. 

The are further compelling reasons that suggest a quantum mechanical formulation 
of general relativity is necessary. Most prominently, general relativity inevitably forms 
singularities in the form of black holes or the initial big bang singularity [47J, at such 
times the effect of quantum mechanics is expected to become important at the very least 
in the matter sector and most likely in the gravitational sector too. 

Given that quantum mechanics provides the recipe of quantisation to produce a quan- 
tum mechanical theory starting with a classical one, the obvious way to construct quantum 
gravity is to apply this procedure to the theory of general relativity. In doing so one is 
immediately faced with a problem; what degrees of freedom should be quantised? The 
answer depends on what one considers to be the fundamental quantity whose dynamics 
general relativity describes. Furthermore, once we have decided upon the form the fun- 
damental degrees of freedom take we must also choose a method of quantisation. For our 
purposes there two possible methods of quantisation, canonical and path-integral, and 
each may be applied to the theory resulting from our choice of degrees of freedom. As one 
might imagine this leads to a large number of different approaches to quantum gravity; 
all of which have been pursued at one point or another. We now consider some common 
answers to the above choices. 

The spacetime metric g is fundamental. In this approach we quantise the metric of 
spacetime. In a canonical approach this involves rewriting general relativity in terms 
of a spatial metric together with functions that describe how the spatial metric evolves 
between each spatial hypersurface. The resulting theory may be quantised leading to 
the Wheeler-de-Witt equation [27J 08]. A path integral approach to the same problem 
requires one to make sense of integrals of the following form, 





where 
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where fi is a constant known as the cosmological constant and k is a constant satisfying 
k oc 1/G, where G is Newton's constant. Additionally one must decide exactly which 
geometries the integral sums over. 

The perturbation to Eat space, h defined by g = rj + h, is fundamental. In this view, 
the vacuum of general relativity is Minkowski space and we quantise the perturbations 
around it. This reduces the problem of quantum gravity to a problem of quantising 
a spin-2 quantum field theory in flat space. One can approach this using canonical or 
path-integral methods, however in both cases one finds that the resulting theory is pertur- 
batively non-renormalisable [1]. Although this would appear to rule out treating general 
relativity as anything other than an effective theory one suggestion due to Weinberg is 
that the coupling of the theory may in fact run to a non-Gaussian fixed point in the UV 
[2]. This would render the theory UV complete while still appearing non-renormalisable 
at the level of perturbation theory. A second possibility arises in the program of string 
theory in which all fundamental objects are treated as quantised strings. In this approach 
the spectrum of the string contains a spin-2 state and the S-matrix for the scattering of 
this state reduces to that of the S-matrix obtained for graviton scattering when treating 
h as fundamental pf9l loTj . 

No variables in general relativity are fundamental. A more radical possibility is that 
general relativity doesn't contain the correct microscopic degrees of freedom at all. If this 
is the case then general relativity only becomes a description of spacetime in the long 
distance limit and attempting to quantise general relativity would be as misguided as 
attempting to understand the short distance regime of water by quantising the Navier- 
Stokes equation! Historically the problem with this possibility is that it was very hard to 
know where to begin as one must invent a theory in which spacetime and gravity appear 
as low energy concepts; in effect solving the problem of quantum gravity in a moment of 
brilliant insight. However, the advent of string theory has changed this situation. String 
theory provides a way to construct order by order a perturbative expansion in which each 
term is finite and correctly reproduces the low energy scattering of gravitons found in the 
linearised theory. Furthermore, this is achieved by quantising a theory which makes no 
assumptions about the classical dynamics of spacetime. The main problem however is the 
question of how to define the theory from which the perturbative series of various string 
theories arise. The solution to this problem has been tentatively named M-theory [50] . 
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however a full definition of it has yet to be given. 

Thesis Outline 

As one can see, the field of quantum gravity is vast and one can only hope, in a thesis such 
as this, to consider a small piece of it. In this thesis we will be particularly interested in the 
dimension of spacetime and the allowed boundary conditions of two dimensional quantum 
gravity. Quantum gravity in two dimensions is special for a number of reasons, one being 
that a number of the approaches listed above coincide. In particular the approach in 
which the metric g is fundamental coincides with the string theory approach and so by 
studying one of these theories we may study the other. 

In chapter [2] we introduce the theory of pure two dimensional quantum gravity in 
the metric-is-fundamental approach. We review how this theory may be obtained as a 
scaling limit of a discrete theory known as Dynamical Triangulation (DT) and introduce 
the matrix model and combinatorial techniques that may be used to compute certain 
observables such as the disc-function and string susceptibility. We also introduce Causal 
Dynamical Triangulation (CDT) models and review their relation to DT. 

In chapter [3] we consider the final observable not considered in chapter [2j the dimension 
of spacetime. In a theory of quantum gravity one expects that spacetime at very short 
distances is violently fluctuating. Furthermore, in string theory one might expect to also 
begin to probe compact extra dimensions. Both of these effects would lead the dimension 
of spacetime to deviate from its infrared value. We consider this question in the context 
of the approach in which the metric is fundamental. We argue that a-priori it is not 
necessarily equal to the dimension of the underlying discretisation and introduce the 
concept of the Hausdorff dimension as a tool to characterise fractal structures. We review 
the calculation showing that the Hausdorff dimension of DT is four and then argue that 
the Hausdorff dimension in CDT is two. We then argue that the Hausdorff dimension is 
inadequate to fully distinguish between smooth and fractal spaces and so we introduce a 
second definition of dimension, known as the spectral dimension, which is based on the 
properties of the Laplacian operator on the fractal structure. We argue that the spectral 
dimension can be obtained by considering random walks on the triangulation. Finally we 
review the recently observed phenomenon of dimensional reduction, in which the spectral 
dimension of the UV is lower then IR. By considering random graphs derived from CDT 
we then begin the search for a simple random graph model in which we might capture 
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the phenomenon of dimensional reduction. 

Chapter H] consists mainly of work published in [66J . We briefly introduce random 
combs and review some known results. We give a definition of the spectral dimension and 
then explain how this definition can be extended to show different spectral dimensions at 
long and short distance scales. We then introduce a simple model which we prove does in 
fact exhibit a spectral dimension that is different in the UV and IR. We generalise these 
results to combs in which teeth of any length may appear with a probability governed 
by a power law and examine the possibility of intermediate scales in which the spectral 
dimension differs from both its UV and IR values. Finally we analyse the case of a comb in 
which the tooth lengths are controlled by an arbitrary probability distribution and show 
that continuum limits exist in which the short distance spectral dimension is one while 
the long distance spectral dimension can assume values in one-to-one correspondence with 
the positions of the real poles of the Dirichlet series generating function for the probability 
distribution. 

In chapter [5] we first review how matter coupled to gravity may be realised using 
matrix models. We introduce the Ising matrix model and also consider multicritical 
points. We then move on to quickly review minimal models, string theory, Liouville 
theory and minimal string theory. The notion of a boundary condition in string theory 
is promoted to a dynamical object known as a brane. The different boundary conditions 
then correspond to different types of brane. If one wants to obtain a non-perturbative 
description of string theory then the spectrum of branes is important. We pursue the 
question, that arises due to a conjecture of Seiberg and Shih [TF], of how many distinct 
branes exist in (p, q) minimal string. We review the fact that this conjecture fails for 
cylinder amplitudes as found in [36] and conjecture a way in which it could be fixed. The 
final portion of this chapter appears in [67] . 

Chatper [6] consists of the remainder of the work in [67] in which we apply matrix 
model techniques to the question of whether the conjecture of the previous chapter holds. 
In particular we consider the matrix model which describes a (3, 4) minimal string theory 
in the presence of a boundary magnetic field. We show how a general class of matrix 
models, which includes the one describing the (3, 4) string may be solved, and by varying 
the boundary magnetic field we reproduce the conformal boundary states of the (3, 4) 
model. Using our general solution to this model we then compute the disc-with-handle 
amplitude for all conformal boundary conditions and compute their deviation from the 
Seiberg-Shih relation. We argue that these results show that our previous conjecture does 
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not hold. This makes it very difficult to see how the Seiberg-Shih relations could possibly 
be true. Finally we consider a different approach to testing the Seiberg-Shih relation based 
on expanding boundary states in terms of local operators. We find that the expansion in 
local operators appears to be valid for only certain boundary types. We support this by 
reproducing a recent calculation of Belavin in which the one-point function on the torus 
was computed. 



CHAPTER 2 



Discrete Approaches to 2D Quantum Gravity 



In the previous chapter we saw that there were a variety of approaches to quantum gravity 
which could be classified according to which degrees of freedom were considered to be 
fundamental. In two dimensions a number of these approaches coincide; in particular 
string theory and approaches in which the metric is considered fundamental. We will 
therefore restrict our attention to these in the remainder of this thesis. In this chapter 
we will consider the simplest case of the metric-is-fundamental approach, in which no 
matter is present and the only dynamical quantity is that of gravity itself. We first review 
why in two dimensions some of the problems of higher dimensional theories are alleviated 
before going on to review the few simple observables that exist in this theory. We then 
motivate a different approach to computing the path integral based on a discretistion of 
the contributing geometries and show that there exists a scaling limit which allows us to 
compute some of the observables. 

In the metric-is-fundamental approach we must evaluate the partition function, 



where we have divided by the volume of the group of diffeomorphisms of the spacetime 
manifold, since this is a gauge symmetry in general relativity. The usual strategy to 
evaluate path integrals is to Wick rotate to a Euclidean theory. This has the effect of 
converting the oscillatory measure to one which is exponentially damped. An immediate 
problem arises in the case of gravity, in that the resulting Euclidean action is not bounded 
from below and so it appears the integral does not exist. A further problem is that it is 
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unclear that a Wick rotation from the Euclidean theory back to Lorentzian theory even 
exists. 

Thankfully, when the spacetime At is two dimensional some of these problems are 
alleviated. Due to the Gauss-Bonnet theorem the curvature term in the action is entirely 
topological and so the problem of the unboundedness of the action is no longer present. 
However, restricting our attention to two dimensions does not solve the problem of what 
level of causality violations to allow i.e. how to Wick rotate between the two theories. 
We will come back to this point later in this chapter. For now we will put no causality 
constraint on the geometries appearing in the integral and hence we shall sum over all 
Euclidean geometries. Application of the Gauss-Bonnet theorem to ( 11.21) then gives, 

Z = e 4 ™* J [ M exp (- J d'a^iaj^ , (2.2) 

where x is the Euler characteristic of AX and is given by x = 2 — 2h, where h, the genus, 
is the number of handles of Ai. In order to perform this integral we must define the 
measure [dg\. This is a difficult technical problem whose solution results in a theory in 
which the overall scale factor of the measure acquires dynamics. The theory describing 
the dynamics of the scale factor is itself a difficult theory to solve. We therefore postpone 
our discussion of this approach until Chapter [51 However, it is still worth considering the 
observables, within the continuum framework, that one could measure in a theory of pure 
quantum gravity. There are a number of observables that have been commonly considered 
in the literature. These consist of, 

1. Integrated correlation functions of local operators are the most obvious observables, 
although these also proved the most difficult to compute. In the most general case, 
when the manifold Ai has n boundaries, the action given in (I2.2p must be extended 
to include a boundary action, 

Z = e 4 ™* J M exp (- - £ dsVhfx B }j , (2.3) 

where h is metric induced on the boundary curve and are boundary cosmo- 
logical constants and \ = 2 — 2h — b, where b is the number of boundaries. By 
utilising the symmetries of the theory any integrated correlation function has the 
form ^ 2 -^l 2 F{ii Bi /^p) PEUE2]. The quantity 7 is known as the string suscepti- 
bility jl2] and is an observable that has been used to compare different approaches 
to two dimensional quantum gravity. 
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2. The Hartle-Hawking wave function. This is the partition function for the appearance 
of a universe "from nothing" and evolving into some prescribed state. Geometrically 
this corresponds to summing over all geometries with a disc topology with a fixed 
boundary condition. Throughout this thesis we will refer to this quantity as the 
disc-function. More generally one can consider multi-loop amplitudes in which the 
spacetime manifold has a number of boundary loops; a disc would be a 1-loop 
amplitude and a 2-loop amplitude is often known as a cylinder amplitude; 

3. The dimension of spacetime. Living inside a particular universe one may probe the 
dimension of the spacetime by performing scattering experiments or even simply 
measuring the electro-static force between two charges. If there exist changes to 
the dimension of spacetime on the length scales such experiments probe, then we 
should be able to measure them. This observable may seem trivial as it appears we 
set the dimension spacetime in our theory at the outset, however, when probing the 
short distance scale of spacetime we expect that spacetime itself will become less 
smooth due to quantum gravity effects (such as microscopic black holes, possibly 
wormholes etc). This could easily alter the effective dimension of spacetime at these 
scales even if the underlying theory is of fixed dimension [12] . 

In order to compute these observables we will consider an alternative method of com- 
puting the above integral in which observables are calculated in a different, but related, 
theory with a far simpler measure. This will come at the expense of having a non-trivial 
dictionary relating the observables in the original theory to the one with a simple measure. 
We will review such an approach in the next section. 

Dynamical Triangulation 

The related theory which we will consider is motivated by considering lattice regulari- 
sations of flat space quantum field theories. If we take the lattice regularised theory as 
our starting point then one can obtain a continuum theory, corresponding to the original 
quantum field theory, by taking the lattice spacing to zero while also appropriately scal- 
ing the other parameters in the theory. In the case of gravity there exists a number of 
approaches to defining a lattice regularisation, however we shall utilise only one here. It 
is known as Dynamical Triangulations (DT) and it has a long history which is reviewed in 
p[2| [53] . In this case the spacetime geometry is encoded in the lattice, which is canonically 
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composed of triangles. The lattice links themselves are always of fixed length q For an 
example of a triangulation see fig 12.11 The benefit of this discrete model is that if we wish 
to compute a partition function, then rather than having to integrate over geometries we 
must sum over the triangulations. Because the latter is a sum, the issue of defining an 
integration measure on the space of geometries does not arise. If we are able to compute 
the partition function of the DT model then the hope is that there exists a continuum 
limit in which we obtain the original partition function (12. 2h . In order to compute the 
partition function we will set each link to be of length a meaning each triangle has an 
area proportional to a 2 . The discretised version of the integral (12. 2p is then, 

oo 

z dt = £ g— aa(t) = £ e -a* m = ^N h {n)g\ (2.4) 

TeTri h TeTri h n=0 

where Tri^ is the set of triangulations of genus h, \T\ is the number of triangles in the 
triangulation T, Mh{n) is the number of triangulations of a genus h manifold containing n 
triangles and g = e~ a2A . From the final equality we see that the evaluation of the partition 
function has been reduced to computing the generating function for A/h(n), which is a 
graph counting problem. However, the graphs we must count are rather complicated; in 

1 Links of fixed length is in contrast with the earlier approach of Regge calculus in which the link 
lengths were integrated over while keeping the triangulation fixed. 
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particular the order of each vertex is unbounded and there is a non-local restriction that 
each face of the graph is triangular. Both of these problems may be solved in a simple 
way by considering the dual graph of the triangulation. The dual graph is constructed by 
associating a vertex to each face of the triangulation and a link joining any two vertices 
whose corresponding faces share an edge. The resulting dual graphs have the simple 
property that each vertex is of order three with no non-local restrictions. Furthermore 
since each graph has a unique dual we need only count the dual graphs. The dual graph 
construction is illustrated in fig 12.11 

From the figure it is amusing to note that dual graph looks very much like a Feynman 
diagram composed of many vertices. In fact there exists a very useful class of tools, 
known as matrix models, which exploit this fact and are able to compute generating 
functions for the graph counting function Nhiji) for a variety of graphs, including those 
with labels. Recall that in quantum field theory the perturbative computation of the 
partition function is organised such that the contribution of each process to the overall 
amplitude is encoded in its Feynman graph; in particular the overall power of the coupling 
associated to each graph is given by the number of vertices it contains. If one were to 
find a way of computing the partition function then the perturbative expansion could 
be recovered by expanding the partition function around zero coupling. If the quantum 
field theory was sufficiently simple then one could use the resulting expansion to extract 
information about the number of graphs contributing at each order of the coupling. 

.1 1 Random Graphs and Matrix Models 

If we are to develop a quantum field theory that will solve our graph counting problem 
then it must also have the property that its perturbative expansion is organised by the 
topology of the graph. This requirement restricts the possible candidates to matrix valued 
field theories as such theories will produce a perturbative series, in the inverse size of the 
matrix, in which only graphs of a particular topology contribute at a given order. 

We must also make the theory simple enough that one can extract the information 
concerning the graphs. The simplest matrix valued theory is one in which the spacetime 
dimension is zero; such a theory is known as a matrix model. An example of such a model, 
which also has the property that it produces graphs in which all vertices have order four 
is. 




(2.5) 
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where the metric field M is an N x N hermitian matrix and the measure [dM] is defined 
to be the flat Lebesgue measure, 



N 



Y[ dM u Y[ dRe [M^] dim [M tj ] . (2.6) 

i=0 i<j 

If we treat the quartic term in the action as a perturbation to the free theory then we may 
construct a Feynman graph expansion for Z MM . The contribution to the total amplitude 
from a given graph is N 2 ~ 2h g n , where h is the genus of the graph and n is the number of 
vertices it contains. We therefore see that the expansion of Z MM around g = coupling 
will yield the series, 

oo 

Z MM {g) = Y, N2 ~ 2hz(h) (9)' ( 2 - 7 ) 

h=0 

oo 

ZF>{g) = J2rth(n)g n . (2.8) 



n=0 



The function that counts the number of graphs with n vertices, A//i(ra), differs from the 
one that appeared in the expression for Z DT since the perturbative expansion includes 
disconnected graphs. In order to compute A4(n) we must instead consider the free energy, 

oo 

\ogZ^ h \g) = Y,^h{n)g n - (2.9) 



n=0 



We see that we have actually done much better than merely being able to compute J\fh{n); 
we have reproduced the Z DT partition function exactly! 



ely 

f 



2.1.2 | Solving the Matrix Model: Loop equations 

There exist a number of approaches to computing Z MM (g), among them the saddle point 
method [HU H2], [53], orthogonal polynomials [121 E3] and loop equations [53], El]. It is 

2 The observant reader will note that the relation of log Z^ h \g) to Z DT is not, in fact, exact, due to 
the above matrix model producing not triangulations but quadrangulations. We could have considered a 
matrix model with a cubic potential and so obtained an exact equivalence with Z BT . However because 
the matrix model with cubic and quartic potentials have the same continuum limit we will continue to 
consider the quartic model for reasons of pedagogy; the quartic model is simpler to solve. 
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the last of these three that we will review in this section. In order to keep our discussion 
general we will consider a generalised one hermitian-matrix model, given by 

Z MM = j [dM]e- NT ^ M \ (2.10) 

where V is a polynomial potential. A loop equation is defined as any equation that 
may be derived from an infinitesimal change in the integration variable which preserves 
the domain of integration. For the matrix model defined previously, the most arbitrary 
infinitesimal change in the integration variable would be M — > M + ef(M) where / is 
an arbitrary function satisfying /(M') = f{My . In the case of (12.101) the effect of this 
change of variable is, 

Z MM = J [dM](l + eJ{M) + 0( e 2)) e -ivw(M)-^(Af) + o(^) ) (2n) 

where J is the trace of the Jacobian matrix for the transformation and K is the change 
in the potential, i.e. in this case K = NTrV (M)/(M). By equating the above equation 
to H 2 . 101) and considering the order e terms we obtain, 

(J) = (K) . (2.12) 

A central role in the loop equation method is played by the function defined by 

N ' n=0 

where w n (g) = (TrM n ) and Wm{x) is known as the resolvent. Its importance is due 
to the fact that the free energy may be computed from it. Furthermore, a very definite 
interpretation may attached to it. Note that the Feynman graphs contributing to w n (g) 
are those which have n external legs. Such graphs may be interpreted as corresponding 
to a discretisation of a disc and so we see the resolvent acts as a generating function for 
such amplitudes. We therefore have, 

oo 

w k (g) = Y t tf{n,k)g n , (2.14) 

n=0 

where JV(n, k) denotes the number of triangulations containing n triangles and k boundary 
links. Obviously such graphs are the ones which contribute to the disc function observable 
discussed in the previous section and therefore one might suspect that the continuum disc 
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function can be extracted from the resolvent. This indeed turns out to be the case and 
we will compute the disc function in the next section. Finally it is worth emphasising 
at this point that although the resolvent function appears naturally in the matrix model 
approach, the technique of constructing a generating function for particular classes of 
amplitudes, for example in this case disc-amplitudes, will be of more general use, as we 
will see when considering the combinatorial approach to loop equations. 

In order to use (12.12p to calculate Wm we first need to specify the change of variable 
we will use. We make the choice f(M) = 1/ (x — M). By analysing the effect of the change 
of variable on the measure (12. 6p . one can show [21] that if f{M) = l/(x — M) then, 

1 „ 1 



J = Tr 



:Tr- 



x- M x-M 



Using this result in (12. 12ft we get, 

— r (Tr Tr 

N 2 \ x-M x-M 



1 

N 
V'(x 



Tr 



iV 



V\M) 



Tr 



1 



x-M 



x 



M 



Q(x) 



(2.15) 

(2.16) 
(2.17) 



where Q(x) = -k (Tr[ — x ~1m ]) is a polynomial in x. We now note that the function on 
the left hand side has contributions from connected and disconnected Feynman diagrams. 
For such amplitudes we can in general write (TtATtB) = (TrA) (TrB) + (TyATyB) c , 
where the c stands for "connected". We therefore get, 

1 



:Tr- 



V'(x)W M (x)-Q(x). 



(2.18) 



M x-M, 

This equation is the loop equation associated with the change of variable defined by 
f(M) = l/(x — M). Let us define the two-loop function as, 

1 „ 1 



W M ;m(x U X 2 ) 



Tr- 



-Tr- 



x 



(2.19) 



M x-M, 

We see that since the diagrams contributing to this amplitude are cylindrical then it goes 
as ~ 1 as iV — > oo. Furthermore, we expect that both Wm and Wm-,m wm have a large 
iV expansion of the form, 



W M (x) 



Wm;m(xi,X 2 ) 



J2N' 2n W^(x), 



n=0 

oo 



(2.20) 
(2.21) 



n=0 
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If we substitute these expansions into the loop equation then we find, 



wS\x) 2 = V'(x)wS\x)-Q(°\x). 



(2.22) 



We see that in the large N limit we end up with a quadratic equation for W M \x) in 
which Q^(x) contains a number of unknown constants. Note that this equation defines 
an algebraic curve in the variables x and W$ . This is known as the spectral curve of the 
matrix model. 

Up to this point we have made an assumption that the partition function for the 
matrix model admits an expansion in inverse powers of N. This assumption is in fact 
sometimes incorrect, as oscillatory terms appear in the large N expansion. The exact 
condition under which a large N expansion and therefore a topological expansion exists 
is given in [21] to be when the spectral curve is genus zero. We see that for the spectral 
curve of the matrix model considered here, the genus zero condition is equivalent to there 
being a single branch cut. By enforcing this condition we may compute the unknown 
constants appearing in Q^(x). The solution to ( 12.22ft is 



Note the discriminant of this is a polynomial in x. We may therefore enforce the genus 
zero condition by requiring that all but two of the roots of the discriminant are double 
roots. For the particular potential in ( 12. 5p . we therefore must be able to write, 



where a(g), A(g) and B(g) may be determined by requiring that the expression in (12.240 
reproduces the large x expansion of ( 12.230 . Because Wm(x) was defined using a series 
valid for large x we know that there exists at least one branch of the solution which has 
the property that wj^\x) ~ 1/x as x — > oo. If one expands ( 12.240 about x = oo, then 
for generic values of a(g), A(g) and B(g) the leading order term will be of higher order 
than 1/x. Requiring that the leading order term is 1/x gives the following expressions for 
"(fiO, A (g) and B(g), 




(2.23) 



W^(x) = \ [V(x) - (a(g)x 2 - A{gf)yf* - B(gf 



(2.24) 



A(gf 



9 



1 + ^{gfg 



1 



B(g) 2 (A + 3B(g) 2 g) 



(2.25) 



1. 
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2.1.3 | The Continuum Limit 

In order to obtain a result that corresponds to a continuum theory we must take what 
is known as a scaling limit of the above computation. This corresponds to tuning the 
coupling of the theory to the value at which it undergoes a second order phase transition. 
At such a phase transition the average number of triangles in the triangulation diverges 
and by rescaling their size along with the couplings in theory we can derive the continuum 
theory which describes the model at this critical point. We see from (I2.4p that the number 
of triangles in the partition will diverge when g reaches the radius of convergence of ( 12. 411 . 
This radius of convergence can be easily computed by considering the quantity 102(g) 
which can be obtained from the coefficient of x~ 3 in the large x expansion of Wm\x). 
The radius of convergence for 102(g) can be read off from the position of its branch points 
as a function of g. 

Having found the critical point of g, g c , it remains to fix how we scale the parameters 
in our theory. If we refer back to ( 12. 4ft then we expect the bare cosmological constant A 
to be related to the continuum cosmological constant \x by, 



which implies g = g c e~ a M , where g c = e . As it stands, if we were to substitute this 
scaling relation in to the equation for the resolvent then this would correspond to taking 
a limit in which the number of triangles in the bulk diverges while also shrinking each 
triangle to zero size. Such a limit would yield a continuous geometry in the bulk. However, 
note that we have not scaled the parameter x, which controls the number of links on the 
boundary of the disc. Leaving the parameter x unsealed would produce geometries which, 
although continuous in the bulk, have a boundary of finite size in lattice units. From the 
perspective of the bulk such a boundary would appear infinitesimal and therefore look 
like a local operator insertion. We will denote the geometry which has a boundary of 
length n in lattice units by w n (g). 

The resolvent provides one way of finding a scaling limit in which the number of 
boundary links can be made to diverge. By noting that the resolvent has a radius of 
convergence in x determined by the position of the branch point B(g), we can see that 
by tuning x to the value B(g c ), the average number of boundary links will diverge. We 
therefore enforce a scaling of x = x c e aC - , where ( will become the boundary cosmological 
constant. Inserting these scaling forms for g and x into the solution for W^\x) and 
taking a — > 0, we obtain, 



A = a 




(2.26) 
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z + Vz 2 - 1 



3/2 



+ (z- V^ 2 - 1 



') 



3/2' 



+ C(a 3/2 ), (2.27) 



where we have introduced the dimensionless quantity z = C/y/J^- 

A peculiarity of this scaling limit is that the universal quantity that corresponds to 
the continuum theory amplitude is not the leading contribution to Wj^ as a goes to zero. 
The obvious way to identity the universal result is to compute the same quantity using 
different potentials in the matrix model and then note which term in the scaling limit 
remains unchanged. A way to short-cut this procedure is to conjecture that the universal 
part of the above expansion corresponds to the first term non-analytic in the bulk and 
boundary cosmo logical constant. One can motivate this by arguing that non-analyticities 
only appear in the thermodynamic limit, i.e. when the surface and boundary are of infinite 
area and length in lattice units. This conjecture is borne out by comparison to continuum 
calculations. We therefore have 



where we have introduced the notation of adding a tilde to a quantity in order to denote 
its continuum version. Note that this continuum quantity is related to the term appearing 
in (I2.27P by a wavefunction renormalisation which removes the overall factor of a 3//2 . 

This completes our computation of the continuum resolvent, and hence the disc func- 
tion observable, using the matrix model formulation of DT. In particular we have shown 
that there exists a non-trivial limit of the discrete equations. From this we may also 
extract the string susceptibility in the case of a disc by removing the marked point on 
the boundary. The mark may be removed by integrating the above expression w.r.t ( to 
obtain, 7 = — |. Another question one might ask is whether there exists a different disc 
function obtained by scaling k in the equation Wk{g) while tuning g to its critical value. 
This would yield a disc function in which the boundary length, rather than boundary 
cosmological constant was fixed. This quantity will turn out to be useful in the future 
and so it is important to understand how it is related to the disc function. Given that the 
boundary cosmological constant has dimensions of a -1 we should expect the boundary 
length to scale like k = A/a, we therefore define the continuum disc-function, 




(2.28) 



W$\\,fi) = \iviaPx- c xla w x/a (g c e- 



(2.29) 
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where r] and x c are chosen to give a non-trivial limit. With this definition it is reasonably 
straightforward to see that wj^\\, jj) is related to W^\(,/j,) by a Laplace transform, 



Jo 

Note that the above relation applies to boundaries which carry a marked point. If we 
were to remove the marked point then the relationship would become, 



where we have used an oo to denote amplitudes with non-marked boundaries. 

The matrix model approach is particularly useful if one wants to compute higher genus 
amplitudes such as discs with handles. This is due to the fact that the loop equations of the 
matrix model produces relations between quantities which include all genus contributions 
to an amplitude. This was seen in f )2.18p . which was then used to derive an equation 
satisfied by the genus zero contribution (I2.22p . Obviously, the loops equation ("12. 18[) may 
also be used to obtain expressions for the higher genus corrections to and this will 
be of great interest in later chapters. 

The question of whether this continuum theory we have uncovered is at all related 
to the original continuum theory (12.21) . is a question we leave until Chapter [51 How- 
ever, it is worth jumping to the punchline; matrix model computations do indeed match 
the continuum formulation whenever two quantities have been computed in both. It is 
therefore worthwhile spending some time investigating what we may learn in the discrete 
formulations of two dimensional quantum gravity. 

.4 | Combinatorial Interpretation of the Loop Equation 

We have seen that the loop equations as obtained from the matrix model are an efficient 
way to compute the disc function and that the matrix model is particularly suited to 
computing higher genus amplitudes. However, it is useful to consider a distinct method 
of obtaining the loop equations which, although it is not as useful for computing the 
higher genus amplitudes, gives greater understanding of the content of the equations. 
This method consists of performing a more direct combinatorial analysis of the trian- 
gulations by considering the way one may build up a triangulation of a given size from 
smaller triangulations. This eventually leads to a recursive equation for the number of 
triangulations of a given size which is equivalent to the loop equations. 




(2.30) 




(2.31) 
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Figure 2.2: Move 1 relates an unrestricted triangulation (left), with a mark denoted by 
the white circle, to another unrestricted triangulation with less triangles. 




Figure 2.3: Move 2 relates an unrestricted triangulation (left), with a mark denoted by 
the white circle, to two disconnected unrestricted triangulations. 

Rather than trying to compute the combinatorics of the triangulations exactly it is 
useful to consider a slightly generalised class of triangulations known as unrestricted 
triangulations in which double links are allowed and triangles do not necessarily have to 
join along an edge. We also require one of the boundary edges to carry a mark. Such 
a triangulation is shown in fig 12.21 Following the arguments of jl3] we will see that the 
continuum limit of these triangulations fall into the same universality class as the usual 
triangulations in DT and indeed give rise to the same loop equations. 

We can classify the graphs with n triangles and k boundary edges into two sets; those 
in which the marked point occurs on a link which forms a side of a triangle and those on 
which the marked point resides on one of a pair of double links. All graphs in the first 
class may be put into a one to one correspondence with graphs with n — 1 triangles and 
k + 1 edges. This is achieved by the map defined by adding double links to each side 
of the marked triangle which is an unmarked boundary and then removing the marked 
triangle. The mark is then transferred to the most anti-clockwise new boundary edge. 
This is shown in fig 12.21 

All graphs in the second class may be put in one to one correspondence with pairs of 
triangulations which are obtained by cutting the triangulation along the marked double 
link and marking the new triangulations as shown in fig 12.31 

This the means that the number of graphs with n > 1 triangles and k > 1 edges 
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Figure 2.4: A causal triangulation. 



satisfies 

n k-2 

Af(n, k) = Af{n - 1, k + 1) + X^( m > W« - m,k — 2 — I), (2.32) 

m=0 1=0 

which when substituted into (12.141) and (12.131) yields the equation, 

W^\zf + (gz 2 - z)W$\z) + 1 - g( Wl (g) + z) = 0. (2.33) 

This equation is the same as ("12.221) with a particular cubic potential. We therefore see that 
the combinatorial approach has yielded the correct loop equations. In the next chapter 
we will use the combinatorial approach to obtain an expression for quantities not easily 
computable in the matrix model approach, for instance the dimension of spacetime. 

Causal Dynamical Triangulation 

Moving to two dimensional gravity alleviated a number of technical problems in the 
metric-is-fundamental approach, however the problem of what level of causality violations 
to allow remained. In the last section we considered the case when no causality conditions 
were enforced, technically this was because we summed over all Euclidean metrics before 
continuing back to a Lorentzian metric. We found that we could compute the partition 
function by considering a discretised version of the path integral. In this section we will 
consider the possibility of computing the Lorentzian path integral directly by means of 
triangulations that are inherently Lorentzian. This will also allow us to explore the role 
causality violations play in the theory. 



2.2 Causal Dynamical Triangulation 
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Figure 2.5: A triangulation of a baby 

universe. If we successively remove tri- pigure 2 g, A typic&1 geometry in 
angles, then the final link between the dynamical triangulation exhibiting the 
baby universe and the parent will be a ^ of & gmooth continuum limit m 
double link. 

Following [4j we will construct a Lorentzian triangulation by requiring all triangles 
composing it to be flat Minkowski space, with one spacelike and two timelike edges. Fur- 
thermore we restrict the topology of spacetime to be cylindrical with spacelike boundaries. 
We define the set of triangulations that will be summed over by giving an algorithm for 
their construction. Suppose the initial boundary is composed of / vertices joined by space- 
like links, then we may attach future pointing timelike links to each of these vertices to 
obtain the set of vertices that, by joining them with spacelike links, form the next spatial 
slice. Furthermore, we require that each vertex has at least one timelike link ending on 
the same vertex as the rightmost timelike link of the previous vertex. This result of this 
procedure is shown in fig 12.41 

Again the problem of computing a disc function can be solved using a matrix model 
technique [15] or by the combinatorial approach. Although results from CDT will motivate 
later work in this thesis we will not have cause to compute any actual CDT observables. 
For this reason we will merely review the results of these computations and refer the 
interested reader to the literature for the technical details. 

One can obtain a qualitative understanding of the difference between DT and CDT 
by considering which geometries are excluded from the path integral in CDT. From the 
construction for the causal triangulations given above we can see that there is no point at 
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Figure 2.7: A typical geometry in CDT [72]. 



which the spatial slice can change topology. In terms of the combinatorial interpretation 
this means that the casual triangulations contain no double links. In the triangulations 
discussed in DT the double links were able to join two disconnected triangulations to form 
a larger one. If one were to view the evolution of a DT universe such a joining of two 
disconnected triangulations would correspond to part of the universe budding off from the 
main one to form a baby universe fig 12.51 Together with topology changes, baby universe 
production is the only type of geometry excluded in CDT which is present in DT [I] . 

The difference between DT and CDT then depends on how much difference baby 
universe production makes to the resulting continuum theory. It can be shown that the 
effect of baby universe production is dramatic; it completely dominates the continuum 
limit of DT, with a baby universe forming at every point in the space [I] . The resulting 
spacetime is highly non-smooth and fractal like as can be seen in fig 12.61 In contrast in 
CDT there exist no baby universes in the continuum limit. This leads to spacetimes such 
as the one shown in fig 12.71 

These result suggest we should consider some method by which the "fractalness" of 
space can be measured. In fact the observable we have yet to consider, the dimension of 
spacetime, is perfect for this role and will be discussed in detail in the next chapter. 



CHAPTER 3 



The Dimension of Spacetime 



In the previous chapter we introduced the matrix model method, by which the disc func- 
tion and string susceptibility of DT may be computed. Of course there was one observ- 
able discussed in the previous chapter that we did not analyse using the matrix model 
approach; the dimension of the spacetime. This observable is in some sense the most 
interesting, as if we were a being living in such a universe the dimension is a far more 
accessible observable than the disc function. 

Finding the dimensionality of the spacetime may seem trivial as we defined the theory 
using two-dimensional building blocks; indeed we even refer to the theory as two dimen- 
sional quantum gravity. However, we will see that the situation is more subtle than this. 
To give an idea of why this is the case one can consider the paths contributing to the 
propagator of a particle in the path integral formulation of quantum mechanics. Although 
all paths are built piecewise from one dimensional lines the entire path becomes "fuzzy" 
as we take the continuum limit. Is the resulting fuzzy path a one or two dimensional 
object? One would be tempted to answer that it somehow has a dimension in-between 
these two values. This notion of fractional dimension can be made precise in a number of 
ways and in this chapter we will see how such concepts may be applied to the quantum 
gravity models defined in the previous chapter. 
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1 | The Hausdorff Dimension dn 

One manner by which the dimensionality of a space can be defined is by observing how 
the volume of a ball scales with its radius. In smooth n-dimensional space we have 
V(r) oc r n , which motivates the definition of the Hausdorff dimension for a metric space. 
The Hausdorff dimension is defined such that if the Hausdorff dimension of the space is 
dn, then for all points the following holds; 

\B(r)\~r dH r^O, (3.1) 

where |-B(r)| is the volume of the ball of radius r. In the models of quantum gravity we 
are considering we are not interested in the dimensionality of a single space but rather 
the expectation value for the dimensionality of spacetime. We therefore must consider, 

(\B(r)\) ~ r dH r^O. (3.2) 

The quantity B(r) in the above equation may be computed for a given geometry i.e. 
metric, contributing to the path integral and the above statement is a statement about 
its mean behaviour after summing over all geometries. We will now review a number of 
ways that this quantity may be computed in DT and CDT models. 

.1 1 Computing dn for DT and CDT 

We will begin by assuming that the Hausdorff dimension is not affected by the global 
topology of the spacetime; this allows one to consider a particular amplitude for which the 
calculation is easier. In order to compute the Hausdorff dimension we need to understand 
how to compute the average area of a spacetime of fixed topology in addition to restricting 
the geometries included in the amplitude to be of a certain "radius". To compute the 
mean area we need only differentiate with respect to the bulk cosmo logical constant, as 
this will insert the identity operator into the bulk. To enforce a fixed "radius" requires 
finer control of the computation than allowed by matrix models and therefore we must 
instead turn to the combinatorial approach. In order to make the notion of "radius" 
precise we define the geodesic distance between two vertices to be the number of links 
on the shortest path through the graph connecting the vertices. The distance of a point 
from a boundary is defined as the minimum geodesic distance between that point and a 
vertex on the boundary and the distance between two boundaries is defined if all points 
on one boundary have the same geodesic distance to the other boundary. 



3.1 The Hausdorff Dimension dn 



25 



The canonical amplitude to consider when computing the Hausdorff dimension is the 
two point correlation function on the sphere with the restriction that the geodesic distance 
between the operators is fixed to t (HI |56]. This means our computation differs slightly 
from the one that would be performed if we wanted to reproduce the definition of dn 
given in ( 13. 2D . Instead we will compute how the total volume of the spacetime depends 
on the distance between two maximally separated points as the distance is made large. 

Following [731 US] , this amplitude may be obtained by considering the more general 
amplitude, G(h,l 2 ,t), corresponding to a cylinder amplitude in which the geodesic dis- 
tance between the entrance and exit loops, of length l\ and l 2 respectively, is fixed to 
be t. The two point amplitude may then be obtained by considering the case when both 
entrance and exit loops shrink to zero size. Once we obtain the continuum two-point func- 
tion G(t, fx) for operators separated by a distance r we can then compute the Hausdorff 
dimension defined implicitly by, 

<y(T))=Km- 1 aa ^ ~ T *', (3.3) 

where V(r) is the total volume of spacetime containing two points separated by a distance 
of r. For triangulations containing n triangles in which the entrance and exit loops are of 
length li and l 2 respectively and are separated by a distance t we may derive a recursion 
relation for the number, J\f{n,l\,l 2 ,t), of such triangulations by the same combinatorial 
methods used for the disc function [44J. There are two possible ways to construct such a 
triangulation from smaller triangulations depending on whether we add a new triangle or 
double link, as detailed in fig 12.21 and fig 12.31 This leads to the recursion relation, 

n h-2 

n — m,lx — 2 — k,l 2 ,t), (3.4) 

m=0 k=0 

where the factor of two in the final sum is due to the fact that the exit loop may appear 
in either of the two disjoint triangulations being joined. In terms of J\f(n, h, l 2 , t) the two 
loop amplitude G(li,l 2 ,t) is given by, 

oo 

G(h,l 2 ,t) = ^2Af(n,h,h,t)g n . (3.5) 

n=0 

Consider a graph with n triangles and l\ boundary links in the entrance loop. If we were 
to add li more triangles to the entrance loop evenly along the boundary, thereby adding 
a new layer to the triangulation, we would have increased the distance separating the 
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entrance and exit loops by one. Therefore, if we were to add only a single triangle to the 
graph then we would have increased t by roughly 1/1%. Combining this with ( I3.4p . leads 
to, 

G(hJ 2 ,t + l/h) = G(h,l 2 ,t) + h t G{h,l 2 ,t) (3.6) 

h 

oo 

= gG(h + l,l 2 ,t) + 2j2^n(g)G(h-n-2,l 2 ,t). (3.7) 

n=0 

Finally, if we now introduce the 2-loop resolvent defined as, 

oo oo oo 

W M . M (z u Z2,t) = J2J2J2 Af ^ l ^ l ^ t ^ nz i h ' lz 2 h ~\ (3-8) 

n=0 ii=0 l 2 =0 

we have 

^W M . M (zi,z 2 ,t) = A [(z 1 -gzf-2W$(z 1 ))W M M z ii z 2i t )] , ( 3 - 9 ) 

where Wm{z\) is the disc function as obtain using the combinatorial approach in the 
last chapter. To find the appropriate initial condition for Wm-,m{zi, z 2 ,t) we first con- 
sider the initial condition for G(li,l 2 ,t). We require G(li,l 2 ,0) = 5^, which implies 
Wm-,m{zi, z 2 , 0) = (-Zi-22 — l) -1 . To obtain a equation for continuum quantities we must 
take a scaling limit. This requires us to decide how t will scale in addition to the sec- 
ond boundary parameter z 2 . One can see that in order for the initial condition to have 
a non-trivial scaling limit we must set Z\ = z c exp(a£i), as usual, but for z 2 we set 
z 2 = z~ l exp(aC 2 )- This gives, 

WWG,C 2 ,0) = -^— , (3.10) 

U + S2 

where we have set Wm-,m(zi, z 2 ,t) = a~ 1 W / M ; M(Ci) (2, t) in order for the scaling dimension 
of both sides of the equation to agree. Using this scaling behaviour for Zi, z 2 and Wm-,m 
in (13.91) we find that the only non-trivial limit occurs when t = resulting in the 

equation, 

d T W M M(uC2,T) = ~d Cl (W$\()WmM(i,(2,t)) . (3.11) 

This equation can be solved for Wm;m((i, (2,^), which together with the initial condition 
fl3TT0l) gives, 

W M M(i,C2,r) = Ci + (2 , (3.12) 



3.1 The Hausdorff Dimension dn 



27 




where / satisfies, 



df(C,r) 
dr 



-W M (f((,r)). 



(3.13) 



The boundaries may be shrunk to points by taking the two cosmological constants to 
infinity and extracting the leading order coefficient, giving the result, 



G(t,h) = {ail) 



3/4 



cosh 




sinh 3 


"(a//) 1/4 r/2 



(3.14) 



where a is an unimportant numerical constant. Using the above equation in (13. 3p gives 
d}j = 4 0. This is quite a striking result as it shows the effective dimension of the space 
is far from what might be expected to arise from two dimensional building blocks. We 
should also emphasise that the above result does not depend on considering a long or 
short distance limit and therefore that dn = 4 on all scales. 

One can avoid some of the above calculation and obtain the Hausdorff dimension 
via more qualitative means by noting that once the scaling dimension of the geodesic 
distance has been decided, then, given that we expect a power law relation between the 



1 A naive application of (|3.3[) will not in fact produce dn — 4 due the presence of a separate length scale 
in the problem; that associated with the cosmological constant. A more transparent calculation of du 
would involve determining the maximal distance between points as a function of /i and then eliminating // 
from an expression of the unconstrained total volume. In our case we simply note that from dimensional 
analysis /i ~ r -4 yielding the stated result for du ■ 
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geodesic distance and area, we can obtain the exponent from dimensional analysis. In the 
above computation we found that r = a 1 * 2 t and since area has dimensions of a 2 this means 
A ~ r 4 . It is interesting to note that the scaling relation r = a 1//2 t indicates that each link 
traversed has an effective size of a 1//2 . This is reminiscent of a random walk and indicates 
that the internal links of the graph form a highly random lattice, such as that shown in fig 
13.11 Intuitively, this random structure is caused by there being baby universes everywhere 
on the lattice, meaning any path through the lattice must pass through many of them. 
This insight allows one to give a reasonable conjecture for the Hausdorff dimension in 
CDT. If one considers the triangulations arising in CDT then due to the sliced structure 
we do not have any baby universes and we expect that the shortest path through the 
triangulation will scale in the same way as the boundary length, i.e. with dimensions of 
a. We therefore expect that the area and geodesic distance will be related by A ~ r 2 , 
giving d,H — 2. This qualitative argument is confirmed by an exact computation using 
combinatorial methods PEJUS]. 

The Spectral Dimension d s 

We saw in the last chapter that due to baby universe production we might intuitively 
expect to obtain a continuum theory from DT that describes the dynamics of a fractal 
spacetime. The last section gave concrete evidence for this in that the mean Hausdorff 
dimension of the continuum spacetime was four despite the lattice theory being defined 
using two dimensional simplexes. However, the fact that the Hausdorff dimension is 
exactly four does raise the question of how we can be sure that the continuum theory 
is indeed fractal; perhaps it is a smooth four dimensional manifold. This question arises 
because the Hausdorff dimension is in fact not that sensitive to the fractal structure of a 
space. One way to see this is to note it is the same for different graphs with greatly differing 
structure, for example, the square lattice and the comb graph obtained by removing all 
its horizontal links except one. The difference is due to the decreased connectivity of the 
comb graph. 

In this section we will introduce another method by which the dimension of a space 
can be defined based on the properties of a propagating particle. This will allow us to 
distinguish the DT continuum limit from a smooth manifold. 

Consider the rate at which the force of an electric charge decreases with distance, it 
clearly depends on dimension. This is related to the form the particle propagator takes 



3.2 The Spectral Dimension d, 



29 



in the corresponding field theory. Both of these things are determined by the Laplacian 
operator of the spacetime on which they are defined and so one would expect that the 
dimension of the space could be extracted from the Laplacian. A simple way to do this is 
to look at some dynamical process whose equations of motion involves the Laplacian and 
look for properties of its solutions that are dependent on the dimension of the manifold. If 
the dynamical process can be defined on spaces other than manifolds then these dimension 
dependent observables can be used to define the dimension for the space. 

The usual choice for the dynamical process is that of heat diffusion. We define the 
heat kernel K(x, x'; a) by, 

(d a -A x )K(x,x';a) = 0, (3.15) 

together with the boundary condition a > and K{x,x'\ 0) = 8(x — x'). On a smooth d 
dimensional manifold the heat kernel has the property that as a — > 0, 

K ( x > x > a )~ (4^p' (3 ' 16) 

from which we can read off the dimension of the manifold. The dimension computed in 
this way is referred to as the spectral dimension of the space. 

The question of generalising the spectral dimension to less smooth spaces now amounts 
to generalising the process of heat diffusion. This is in fact why the heat diffusion process 
was chosen; qualitatively, the heat kernel K(x, x'\ a) can be thought of as a "scaling limit" 
of the probability of a random walk moving between the points x and x' in t steps. By 
a scaling limit we mean a process in which the number of steps in the walk is taken to 
infinity while reducing the size of each step. 

The process of a random walk is obviously something that can be naturally defined 
on a graph, which allows the spectral dimension to be defined for each triangulation 
contributing to the path integral. In particular, for triangulations, we need to compute 
the return-probability pit) for the random walk to leave a given point and return to it in 
t steps. In the scaling limit this quantity will become K(x, x; a) and we then can study 
its behaviour as a — > to find the spectral dimension. 

Obviously the above procedure depends on giving a precise definition of the scaling 
limit. In fact we may extend the definition of the spectral dimension to discrete structures, 
without taking a scaling limit, by considering the t — )■ oo behaviour of pit). The reason 
for studying the t — > oo limit in the discrete case can be understood from considering the 
process of taking the scaling limit, in which the link length of the discrete space is reduced 
to zero. A walk whose length is of finite size measured in link lengths would produce a 
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measurement of the spectral dimension contaminated by cutoff effects. In contrast, a 
short walk in the continuum space will still correspond to a walk whose length is much 
greater than the link length. We therefore expect the large t behaviour of p(t) to fall off 
with the same exponent as the small a behaviour of K(x,x;a). For discrete spaces we 
therefore define the spectral dimension by, 

d , = - 2 , im !«, (3.17) 

t^oo log t 

if the limit exists. 

As we saw with the computation of the Hausdorff dimension, the matrix model ap- 
proach has difficulty computing quantities that are not disc-functions or their generalisa- 
tions. Instead, we had to augment the computation with a combinatorial approach. This 
problem is in fact worse for computing the spectral dimension; this is not surprising as the 
spectral dimension relates to a much more complicated phenomenon than the Hausdorff 
dimension, which simply measures volume. In DT there exists an approach to computing 
the spectral dimension, based on the continuum formulation [57], although we will not 
review this approach in detail here. The result of the computation in |57J is that the spec- 
tral dimension in DT is two. Because the Hausdorff dimension differs from the spectral 
dimension, this confirms our expectation that the continuum limit in two dimensional DT 
is not describing a smooth manifold but a fractal space. 

When it comes to CDT we are mostly stuck for the moment with merely exploring 
the spectral dimension numerically. The result of such numerical simulations is that 
in (i-dimensional CDT (where d = 3,4 has been considered) the spectral dimension is 
consistent with d s = d [68] . For d s = 2 the situation is not so problematic and there 
has been progress, which we will review shortly, in computing the spectral dimension 
analytically. The result of such work suggests that unlike DT, the long distance behaviour 
of the resulting spacetime is that of a smooth d- dimensional manifold. 

3.2.1 1 Dimensional Reduction 

In addition to providing a measure of the dimension of the space, more information can 
be extracted from the analysis of a random walk by controlling the walk length. Since the 
size of the space explored by a random walk of t steps goes as \/i we can probe shorter 
distances by reducing the length of the walk. We therefore can gain some insight into 
the Planck scale behaviour of gravity in CDT and other models. The effective spectral 
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dimension as a function of the length probe has been investigated numerically in three 
and four dimensional CDT with the surprising result that although the dimension is three 
and four respectively on long distance scales, on very short scales the spectral dimension 
in both cases is measured to be two [?]• This phenomenon is known as dimensional 
reduction. 

The phenomenon of dimensional reduction of the spectral dimension is compelling, as 
similar behaviour has been observed in a variety of other approaches to quantum gravity 
in four dimensions such as the non-Gaussian fixed point possibility in the perturbation- 
is-fundamental program [3] and also Horava-Lifshitz gravity [HJ [TD] . Other evidence has 
been discussed in (TTJ [12] . This is very suggestive of dimensional reduction being a robust 
feature of quantum gravity. It will be our task in the next chapter to attempt to construct 
a toy model in which dimensional reduction of the spectral dimension occurs. 

3.2.2 | The Spectral Dimension in CDT 

The triangulations in DT and CDT are formulated mathematically as a set of graphs 
together with a measure on that set. We will refer to a set of graphs together with a 
measure as a random graph or an ensemble of graphs. To construct a toy model exhibiting 
dimensional reduction we will want to consider a random graph in which a scaling limit 
can be taken. Furthermore, we will want to keep the ensemble of graphs as closely related 
to the CDT graphs as possible. Before doing this it is worth pausing to draw inspiration 
from CDT. We will briefly review one attempt to derive the spectral dimension in two 
dimensional CDT. This will provide clues at to what sort of graphs should be used to 
construct our ensemble. 

As previously discussed, matrix model and combinatorial methods are not appropriate 
for computing the spectral dimension. An alternative approach is to directly analyse the 
random walk. One method that proves useful in this approach is to encode the first- 
return-probability for a random walk in a generating function [31 HHJ HE]. By deriving 
sufficiently tight bounds on this generating function we can obtain an exact value for the 
spectral dimension. This method has the advantage that, when it works, it can provide 
stronger statements than merely averages; it can show that a particular value for the 
spectral dimension holds occurs with probability one [El [16] . It also has the advantage 
of being entirely mathematically rigourous. 

To date it has proven too difficult to carry out the above strategy directly on the tri- 
angulation. Instead attempts have focused on considering ensembles of simpler graphs for 
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Figure 3.2: A causal triangulation 

Figure 3.3: Casual triangulations have 

drawn in the plane. Note we have 

a bijection with trees. In this case the 

added an extra vertex with links to all 

tree associated to the triangulation is 

vertices in the inner most spatial slice 

„ , shown by thick lines. 16 . 

HB]. 1—1 

two reasons; firstly, because they serve as toy examples in which to develop the necessary 
techniques and secondly because a number of simpler graph ensembles can be obtained 
from CDT while preserving some of its properties. In fact it can be shown that certain 
simpler graphs have spectral dimensions that bound the spectral dimension of CDT. In 
particular the graphs obtained by collapsing all vertices of a given height to a single vertex 
have been used to show that d s < 2 almost surely for CDTs [TB]. On the other hand one 
can relate a causal triangulation to a tree graph T by the following procedure [T6]. First 
we recognise that casual triangulations are planar graphs and so can be represented as 
lying in the plane with the spacelike slices forming concentric circles as shown in fig 13.21 
We denote the vertices in the kth spacelike slice as Sk, we define So as containing a single 
new vertex connected to all vertices in Si. The algorithm then maps the triangulation 
into a tree T: 

1. All vertices in the triangulation are in T in addition to a new root vertex that is 
only linked to S . 

2. All links from Sq to Si are in T. 

3. All links from a vertex in Sk to vertices in Sk+i, apart from the clockwise- most link, 
are in T. 
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This algorithm is shown in fig 13.31 In fact it can be shown that the map produced by 
this algorithm is a bijection and furthermore preserves the measure on the CDTs. The 
spectral dimension of the resulting random tree has long been known to be 4/3 and this 
gives a lower bound on the spectral dimension of the CDTs. One tool to obtain results 
about such trees is another random graph known as a comb. It is this random graph that 
we will use in the next chapter to construct a toy example of dimensional reduction. 
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CHAPTER 4 



A Toy Model of Dimensional Reduction 



In the last chapter we saw that both DT and CDT exhibit fractal behaviour in some 
regimes of the theory. However DT appears to be far more pathological as fractal space- 
times dominate the continuum limit. Furthermore we saw that the situation in CDT 
is much improved, with both Hausdorff and spectral dimension agreeing and agreeing 
with the dimension, on long distance scales, of the underlying simplexes. At short dis- 
tances however it exhibits a reduction in the spacetime dimension; a phenomenon that 
has been observed in other approaches to quantum gravity and has been given the name 
dimensional reduction. 

In this chapter we will construct a definition for a scale dependent spectral dimension 
and show that there are models which do indeed exhibit scale dependent spectral dimen- 
sions defined in this way. In particular we develop a simple model based on previous work 
on random combs [3]. These are a family of simple geometrical models which share some 
of the properties of the CDT model; instead of an ensemble of triangulations we have an 
ensemble of graphs consisting of an infinite spine with teeth of identically independently 
distributed length hanging off (we define these graphs precisely in Section 2). It was 
shown in [3] that the spectral dimension is determined by the probability distribution for 
the length of the teeth. In this chapter we show that it is possible to extend the work 
of [3] by taking a continuum limit thus ensuring that the cut-off scale is much shorter 
than all physical distance scales. We find that the spectral dimension is one if we take 
the physical distance explored by the random walk to zero and there exists a number 
of continuum limits in which the long distance spectral dimension differs from its short 
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distance counterpart. As a by-product of this work we also extend some of the proofs 
given in [3] to a wider class of probability distributions. 

This chapter is organized as follows. In Section 14.11 we briefly review some known 
results for combs and their spectral dimension and then explain how in principle these 
can be extended to show different spectral dimensions at long and short distance scales. 
In Section I4T21 we introduce a simple model which we prove does in fact exhibit a spectral 
dimension that is different in the UV and IR. This model forms the basis of all later 
generalisations. In Section 14.31 we generalise the results of Section 14.21 to combs in which 
teeth of any length may appear with a probability governed by a power law. In Section 
14.41 we examine the possibility of intermediate scales in which the spectral dimension 
differs from both its UV and IR values. In Section 14.51 we analyse the case of a comb in 
which the tooth lengths are controlled by an arbitrary probability distribution and show 
that continuum limits exist in which the short distance spectral dimension is one while 
the long distance spectral dimension can assume values in one-to-one correspondence with 
the positions of the real poles of the Dirichlet series generating function for the probability 
distribution. We then show how these techniques can be used to extend the results of [3] . 
In Section 14.61 we discuss our results and possible directions for future work. 

4.1 1 Combs and Walks 

In this section we review some basic facts about random combs and random walks. As 
much as possible we use the same notation and conventions as [3] and refer to that paper 
for proofs omitted here. 

4.1.1 1 Definitions 

We use the definition of a comb given in [3]. Consider the nonnegative integers regarded 
as a graph, which we denote Noo, so that n has the neighbours nil except for which 
only has 1 as a neighbour. Furthermore, let Ng be the integers 0,1, ... ,£ regarded as a 
graph so that each integer n£iV< has two neighbours nil except for and I which only 
have one neighbour, 1 and £ — 1, respectively. A comb C is an infinite rooted tree-graph 
with a special subgraph S called the spine which is isomorphic to with the root at 0. 
At each vertex of S, except the root 0, there is attached one of the graphs Ng or N^. We 
adopt the convention that these linear graphs which are glued to the spine are attached 
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Figure 4.1: A comb. 



at their endpoint 0. The linear graphs attached to the spine are called the teeth of the 
comb, see figure 14.11 We will find it convenient to say that a vertex on the spine with no 
tooth has a tooth of length 0. We will denote by T n the tooth attached to the vertex n 
on S, and by Ck the comb obtained by removing the links (0, 1), . . . , (k — 1, k), the teeth 
Ti, . . . , Tfc and relabelling the remaining vertices on the spine in the obvious way. The 
number of nearest neighbours of a vertex v will be denoted <j(v). 

It is convenient to give names to some special combs which occur frequently. We 
denote by C = * the full comb in which every vertex on the spine is attached to an 
infinite tooth, and by C = oo the empty comb in which the spine has no teeth (so an 
infinite tooth is itself an example of C — oo). 

Now let C denote the collection of all combs and define a probability measure v on 
C by letting the length of the teeth be identically and independently distributed by the 
measure /i. We will refer to the set C equipped with the probability measure v as a 
random comb. Measurable subsets A of C are called events and v(A) is the probability 
of the event A. The measure of the set of combs A with teeth at ni,n 2 , . . . , having 
lengths 4,4, ■ ■ ■ ,4 is 

k 

u(A) = Hii(ij)- (4-1) 

i=i 

For any z/-integrable function F defined on C we define the expectation value 

(F(C)) = J F{C)du. (4.2) 

We will often use the shorthand F for (F(C)). 

1 Since the space of combs is discrete, the integration should be replaced by a sum. However, we will 
continue to use an integral sign, with the understanding that it refers to a sum over discrete structures. 
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.2 | Random Walks 

We consider simple random walks on the comb C and count the time t in integer steps. 
At each time step the walker moves from its present location at vertex v to one of the 
neighbours of v chosen with equal probabilities a(v Unless otherwise stated the walker 
always starts at the root at time t — 0. 

The generating function for the probability pc{t) that the walker is at the root at time 
t, having left it at t — 0, is defined by 

oo 

Q c (x) = J> - x)^p c (t) (4.3) 

t=0 

and we denote by Pc(x) the corresponding generating function for the probability that 
the walker returns to the root for the first time, excluding the trivial walk of length 0. 
Since walks returning to the root can be decomposed into walks returning for the 1st, 2nd 
etc time we have 

= t=w (4 ' 4 > 

It is convenient to consider contributions to Pc(x) and Qc(x) from walks which are 
restricted. Let P^\x) denote the contribution to Pc( x ) from walks whose maximal 
distance along the spine from the root is n and define 

n-l 

P^(x) = J2 P c k) (*) (4-5) 

k=0 

which is the contribution from all walks which do not reach the point n on the spine. 
Similarly we define 

oo 

Pt~ 1) (*) = J2 P c ) W- (4-6) 

k=n 

Clearly Pc(x) can be recovered from P^ <n \x) by setting n — > oo. We define the cor- 
responding restricted contributions to Qc(x) in the same way. By decomposing walks 
contributing to P^ n \x) into a step to 1, walks returning to 1 without visiting the root, 
and finally a step back to the root it is straightforward to show that 

Pc n \*) = l ~ X <, Tl . ( 4 - 7 ) 

Z-P^{x)-P^{x) 

where we have adopted the convention that for the empty tooth, T = 0, 

P $ (x) = 1. (4.8) 
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The relation (14. 7p can be used to compute the generating function explicitly for any comb 
with a simple periodic structure and we list some standard results in IA.1I 

There are a number of elementary lemmas which characterise the dependence of Pc{%) 
on the length of the teeth and the spacing between them [3]. We state them here in a 
slightly generalized form which is useful for our subsequent manipulations. 

Lemma 1 The function Pq (x) is a monotonic increasing function of Px k (x) and P^ n k \x) 
for any n > k > 1 . 

Lemma 2 Pj^ n \x) is a decreasing function of the length, i^, of the tooth for any 
n > k > 1. 

Lemma 3 Let C be the comb obtained from C by swapping the teeth and T^+i, k < 
n — 1 . Then P^ T ^ (x) > P^f (x) if and only if P^ h (x) > Pr k+1 (%) ■ 

The proofs use (14.71) and follow those given in [3] for the case n = oo. 
An important corollary, valid for any comb, of these lemmas is that 



which we will refer to as the trivial upper and lower bounds on Qc{x)- The result follows 
from Lemma [2] with n = oo, which gives 



and the explicit expressions for P*(x) and -Poot^) given in IA.1I 

4.1.3 | Two point functions 

Two point correlation functions on the comb correspond to the probability of a walk 
beginning at the root being at a particular vertex on the spine at time t. In particular, 
let pc{t', n) denote the probability that a random walk that starts at the root at time zero 
is at the vertex n on the spine at time t having not visited the root in the intervening 
period. We will refer to the generating function for these probabilities as the two point 
function, Gc(x; n), and define it by 



x 4 < Qc( x ) < x 2 1 



(4.9) 



P*{X) < P C {X) < Poo{x), 



(4.10) 




(4.11) 
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Gc(x;n) may be expressed as 

n-l 

G c (x; n) = a(n)(l - x)^ 2 J] P Ck (x) (4.12) 

fc=0 

which may be used in conjunction with Lemma [2] to obtain the bounds, 

G*(x;n) G c (x;n) Gj^yn) 

3 - a{n) 2 1 ' } 

Now let rc(t; n) denote the probability that a random walk that starts at the root at time 
zero is at the vertex n on the spine for the first time at time t having not visited the root 
in the intervening time. We define the modified two point function, G^(x;n), by, 

oo 

G° c (x-n) = Y,(±-xY /2 rc(t ] n) (4.14) 
and note the following lemmas; 

Lemma 4 The contribution Pq (x) to Pc{x) from walks whose maximal distance from 
the root is N or greater satisfies 

P^ >N - 1] (x) < 3x-^G° c {x; AO 2 . (4.15) 

The proof is given in section 2.4 of [3J. 

Lemma 5 The modified two point function satisfies 

G { °\x;n) < G ( c\x;n) < G^(x;n). (4.16) 

To prove this note that 

/-, \—( n —2)/2 ™~ 2 

G ^">= ( '(Li) Il P c?-%) (4-17) 

1 ' k=0 



and use Lemma [2] . 
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.4 | Spectral dimension and the continuum limit 

As was discussed in the previous chapter, we may characterise the dimension of spaces 
which are not manifolds by use of the spectral dimension. For spaces obtained as the 
continuum limit of a graph or ensemble of graphs the probability that a random walk 
returns to its initial position generalises the heart kernel. However, given that it is easier 
to work with the generating function for the return probability, Q(x), it is more convenient 
to define the spectral dimension in terms of the behaviour of Q(x) implied by (13.171) . 

Q c {x)~x- l+d °l\, (4.18) 

where by f(x) ~ g(x) we mean that 

cg(x) < f(x) < cg(x), < x < x , (4.19) 

where c, d and x are positive constants. The property (14.181) was adopted in [3] as 
the definition of spectral dimension, assuming it exists. The spectral dimension of an 
ensemble average is defined in the same way, simply replacing pc(t) and Qc(x) by their 
respective expectation values. 

Our present goal is to extend the definition of the spectral dimension to give a math- 
ematical meaning to the notion of a scale dependent spectral dimension as required to 
describe the phenomenon of dimensional reduction. In order for the spectral dimension 
to be scale dependent there must exist a length scale L in the system besides the cutoff. 
Furthermore, in order for this length scale to survive in the continuum, it must be scaled 
as the continuum limit is taken. We will therefore introduce a second parameter into the 
tooth-length probability distribution fi(£), which defines a length scale for some structure 
of the comb. One would hope that in more realistic models such a length scale could be 
generated dynamically. 

We now give a precise meaning to the term continuum limit. We assign the value a to 
the distance between adjacent vertices in the graph and take the limit a — > and L — > oo 
in such a way that the scaled combs have a finite characteristic distance scale; it is this 
limit to which we give the name continuum limit and quantities which exist in this limit 
we call continuum quantities. Walks much longer than L will probe different structure 
from walks much shorter than L but nonetheless both can be very long in units of the 
underlying cut-off scale a. 

In the following sections we will denote dependence of a function on a number of 
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variables Lj □, i = 1, . . . , iV, by Lj passed as one of the function arguments. Given a 
random comb ensemble specified by Li) and the corresponding Q(x; Li) we define the 
continuum limit of Q by, 



where the scaling dimensions A M and A; are chosen to ensure a non-trivial limit and the 
combinations £Aj are dimensionless. As we shall see, the choice for A M and Aj is unique 
given mild assumptions. The function Q can be used to define the spectral dimension at 
short and long distances. 

The variable £ of the generating function plays the role of a fugacity for the walk 
length. Indeed it is the Laplace conjugate variable to walk length, this is completely 
analogous to the way in which the variable for the resolvent in Chapter [2] plays the role of 
the continuum boundary cosmological constant and is related to the boundary length by 
a Laplace transform. We therefore expect £ — > oo to correspond to very short continuum 
walks and £ — > to correspond to long continuum walks. Let us for a moment suppose 
that the long and short distance behaviour of Q(C,; A) is that of a power law and further 
that there exists only one length scale A. If a diffusion experiment were performed in 
which the diffusion time was much greater than A ^> 1 then an experimenter would see 
the power law behaviour associated with the £ — > oo limit of A). Conversely, if the 
walk length was much less then A then an experimenter would observe the power law 
behaviour of the £ — > limit of A). If the exponent of the power law differs in these 
two limits then the experimenter would conclude the spectral dimension differed on short 
and long distance scales. 

We now give an explicit definition for the long and short distance spectral dimension 
and show that the above long and short distance limits only contain contributions from 
long and short walks respectively. We will assume for simplicity that there is just one 
scale L and that the spectral dimension in the sense of ( 13 . 1 71) exists for (pc{t)) which 



Q(frAi) = lim<A 



Q«;a- A *Af'), 



(4.20) 



2 In the following, the parameters Li will becomes length scales related to structures in the comb. 
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implies that there exists a constant to such that (pc(t + 1)) < (pc{t)) , i > to- Note that 



x 



t/2 



Q(x;L)- J2 (Pc(t))(l-x) 



t/2 



t=0 



t=T+l 

oo 



Q(x; L) - (1 - x)( T+1 )/ 2 £ (pc(t + T + 1)) (1 



— a: 



t/2 



f=0 

oo 



> Q(x; L) - (1 - x) (T+1)/2 £ (pc(0> (1 " *) 



t/2 



t=0 



to 



+(1 - ^) (T+1)/2 £(pc(f) - Pc(t + T + 1 ))( 1 



,t/2 



t=0 



Since the final sum is finite it will have an analytic dependence on a. Now choose 

1 



(4.21) 



(4.22) 



£io g (i + ^; 

and set x = a£ and L = a~ A A A in (I4.2ip to get 

cAQ«; a" A A A ) (1 - exp (— £A)) + a A "i2exp (— £A) 

T 

<a A »J2(Pc(t))(l-&Y /2 
t=o 

<a A Q(ae;a~ A A A ), (4.23) 

where R is a constant arising form finite sum in ( 14.2 ip . Provided that the limit in (I4.20p 
exists we see that the behaviour of Q(£; A) as £ — > oo characterizes the properties of walks 
of continuum time duration less than 

1 



lim 



A. 



^oo£l g(l + i) 

and we define the spectral dimension d® at short distances by 

rf : = 2 | 1+lim i2iM|A))' 

^oo log 4 



(4.24) 



(4.25) 



provided this limit exists. 

We can define the spectral dimension at long distances in a similar way. First note 
that by g^J) 

oo , x t/2 T x ,\t/2 / -I \ t r 

v^>E (po(*)> 1 - y ) > E <M*)> f i - y ) > f 1 - t ) E (^( t ))( 4 - 26 ) 

t=o V ^ t=o V / V / t=0 
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so that 

1 oo 

Q{x- L) - VT(1 - -Y T < (Pc(t)) (1 - xf 2 < Q{x- L). (4.27) 

t=T+l 

This time letting T = 1 l°g(l + £A)J — 1 we get 

oo 

cAQ«; cT a A a ) - e^ 1 log(l + £A) < cA ^ (p c (t)) (1 - £a)'/ 2 < a A "Q«; cT A A A ). 

i=T+l 

(4.28 

If the above limit exists and furthermore Q(£;\) diverges as £ — > then we see that in this 
limit A) only receives contributions from walks whose length exceeds lim^o log(l+ 
£A) = A. We define the spectral dimension df at long distances to be 

C = 2fl + lim 1 ° g( 1 ^J A)) V (4.29) 
\ log£ J 

provided this limit exists. 

At this point it is unclear whether there exists any examples for which the short and 
long distance spectral dimensions as in (I4.25P and (14.291) exist and furthermore differ in 
value. We will now discuss some examples based on random combs in which we realise 
such examples. In all the examples given in this chapter it turns out that the exponent 
is A, = \. 



A simple comb 

We now introduce a random comb whose spectral dimension differs on long and short 
length scales and thus illustrates that the behaviour described in section W. 1.41 can actually 
occur. This comb is defined by the measure, 



i_ 

L ■ 

0. 



£ = 0, 
£ = oo, 
otherwise. 



(4.30) 



This random comb has infinite teeth and they occur with an average separation of L. 
Intuitively we would expect that if a random walker did not move further than a distance 
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of order L from its starting position it would not see the teeth and therefore would 
measure a spectral dimension of one. If however it were allowed to explore the entire 
comb it would see something roughly equivalent to a full comb and so feel a much larger 
spectral dimension. To prove this intuition correct we proceed by computing upper and 
lower bounds for Q which are uniform in L and for < x < xq, where the constant xq is 
equal to one unless otherwise stated, and then take the continuum limit to obtain bounds 
for Q . 

With complete generality we may obtain a lower bound on Q(x) by use of Jensen's 
inequality which takes the form, 

Lemma 6 Let Pt{x; Li) be the mean first return probability generating function of the 
teeth of the comb defined by fi(£; Li), then 

Q(x;Li)>(l + x-P T (x;Li))-l (4.31) 

The proof is given in [3J. For the comb f l4.30p we have 



which implies 



P T {x- L) = 1 - 1(1 - P^x)) = 1 - Y ( 4 - 32 ) 



Q(x;L)>(^ + x) . (4.33) 



Letting x = a£ and L = a 2 A 2 gives 



Q(£;\) = lima5Q(a£;a-5A5) > ( v ~h + l) \ (4.34) 

a^O V / 

where we have introduced the dimensionless variable v = £A. 

To find an upper bound on Q(x; L) we follow [3] and use Lemmas [TJ [2] and [3] to compare 
a typical comb in the ensemble with the comb consisting of a finite number of infinite 
teeth at regular intervals. First we define the event 

A(D, k) = {C:Di<D:i = 0, k}. (4.35) 

where Di is the distance between the i and i + 1 teeth and then write, 

Q{x;L) = J Q c (x;L)du 

= [ Q c (x;L)du+ [ Q c (x;L)du. (4.36) 

JC/A{D,k) JA(D,k) 
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Since the Di are independently distributed 

u{A(D,k)) = (l-{l-l/L) D ) k . (4.37) 

Consider a comb C G A(D, fc); then by Lemmas HJ [2] and El 

P c (x;L) <P c >(x), (4.38) 

where C is the comb obtained by removing all teeth beyond the k tooth and moving the 
remaining teeth so that the spacing between each is D. Now we can write 

P ol {x) = P£ m \x) + P£ Dk - 1] (x). (4.39) 

Since the walks contributing to Pjy Dk \x) do not go beyond the last tooth we have 

P^ Dk \x) < P* D (x), (4.40) 

where *D denotes the comb consisting of infinite teeth regularly spaced and separated by 
a distance D. Using (14.401) . Lemmas 0] and [5] we have, 

P c (x; L) < P* D (x) + 3x-*GW(x; Dk) 2 (4.41) 

uniformly in A. P*d{x) and G^(x;n) are given in Appendix A. Now set D = [D\ and 
k = \k], where, 

D = 2L\ logxL 2 |, k = (xL 2 )~ 1/2 . (4.42) 

Since G^(x; n) is manifestly a monotonic decreasing function of n and P*d{x) an increas- 
ing function of D, 

Q(x; L) < x^ 2 (l - (1 - (1 - l/Lf- 1 )'^ 1 ) + Qu(x)(l - (1 - l/L) b ) k (4.43) 
where we have used (14.91) and 

Qu(x) = [l - P^{x) - 3x^G^(x; (D - l)fc) 2 ] _1 . (4.44) 
Taking the continuum limit of (14. 43ft and using the results of IA. II then gives 

Q(Z;\)<C 1/2 Fm, (4.45) 

where 

1 + ofv^ 1 ), v — > oo, , 

F(v) = < 1 ^ j (4.46) 

I v * yj logf 2 1 + o(v 2 ), v — > 0. 

It follows from (14251) . ( 14^91) . (Q4l and (1Q5]) that 

d° s = 1, d°° = |- (4.47) 
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Combs with Power Law Measures 



We now consider slightly more general combs in which the measure on the teeth is a power 
law of the form, 



L > 



= o, 

>0, 



(4.48) 



where C a is a normalisation constant and as before L plays the role of a distance scale. 
We consider laws in the range 2 > a > 1 as it is known that for a > 2 the comb has 
spectral dimension d s — 1 in the sense of (14. 18f) [3] and therefore it is not possible to get 
a spectral dimension deviating from 1 on any scale. 

To compute a lower bound on the return probability generating function for the above 
distribution we apply Lemma and reduce the problem to computing an upper bound on 
1 — Pt{x). The first return generating function Pg(x) for a tooth of length I is recorded 
in ( 1A.5I) ; bounding tanh(-u) above by the function f(u) = u for u < 1 and f(u) = 1 for 
u > 1 Gfive^j 



l-P T (x;Li) < y/x~ 



e=1 i=[m£]+l 



I TTloo 0*0 / ^ ^ 1 

< y/x - rrioo^y/x / y j n(£; Lj) du. 



(4.49) 



To obtain the second inequality we have applied the Abel summation formula. We there- 
fore have, 

Lemma 7 For a random comb defined by the measure fi, 

r i 

1 — Pt(x; Li) < y/x — m 00 {x)\/x I x{ u \Li)du (4.50) 

Jo 

where the cumulative probability function x{u; Li) is defined by xi u ', Li) = X/l=o M^S ^*)- 



3 In this particular case we could in fact compute 1 — Pt(x) exactly by the Abel summation formula. 
However the bound we use is good enough to give the desired result with the advantage that the calculation 
can be done with elementary functions. 
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We will see shortly that all behaviour of the spectral dimension of the continuum comb is 
encoded in the asymptotic expansion of x(tt; Lj) as u goes to infinity @. In the present case 
x{u; L) is trivially related to the partial sum of the Riemann (^-function whose leading 
asymptotic behaviour is well known and we find 

X (u;L) = l-^^—L + 5(u), (4.51) 
L a — 1 

where 

\8(u)\ < yU~ Q , u>2. (4.52) 
It follows that for x < Xo, where m^xo) = |, 

1 - P T (x; L) < moo (x)^ (j-m^xy- 2 + ^m^x)*" 1 + j-^j , (4.53) 

with &i,2,3 being constants depending only on a and &i > 0. Choosing L = a~ A> X A ' with 
A' = 1 — a/2 yields a lower bound on the continuum return generating function, 

A) > T 1/2 (1 + b^X)-^)- 1/2 . (4.54) 

To obtain a comparable upper bound we need 

Lemma 8 For any random comb and positive integers H, D and k, the return probability 
generating function is bounded above by 

Q(x; Li) < x^l\l - (1 - (1 - p) D ) k ) + Qu(x)(l - (1 - p) D )\ (4.55) 

where 

oo 

e=H+i 

Qu(x) = [l-PnM^-^G^&Dkyy 1 , (4.56) 

and Ph,*d is the first return probability generating function for the comb with teeth of 
length H + 1 equally spaced at intervals of D. 



4 In general it is not obvious that this asymptotic expansion exists due to the discontinuous nature of 
X- We will address this issue later when we consider generic measures. 
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The proof is a slight modification of the upper bound argument used in Section 14.21 First 
define a long tooth to be one whose length is greater than H; then the probability that a 
tooth at a particular vertex is long is 

oo 

p= (4-57) 

i=H+l 

Define the event 

A(D,k) = {C : Di < D : i = 0, ...,&} (4.58) 
where now Di is the distance between the % and i + 1 long teeth so that 
Q(x;Li) = J Q c (x;Li)du 

Qc{x; Li)dv + / Q c (x;Li)du. (4.59) 

C/A(D,k) JA(D,k) 

Since the Di are independently distributed 

u(A(D,k)) = (l-(l-p) D ) k . (4.60) 

Now use Lemmas [2] and [3] in turn to note that for 

PceA(D,k)(x,L) <P c ,(x,L) (4.61) 

where C is the comb in which all teeth but the first k long teeth have been removed 
and the remaining long teeth have been arranged so that they have length H + 1 and a 
constant inter-tooth distance D. By the same arguments as we used in Section H~2l to get 
(14.411) we obtain the bound 

PceA(D,k)(x, L) < P H+li , D + 3 X - 1/2 G^(x, Dkf . (4.62) 

Lemma [8] then follows from (14 .9p , (I4.60P and (14.621) . We now specialise to the power law 
measure (14.481) and set H = \_H\ , D = \_D\ and k = \k~\, where 

H = x- 1 ' 2 

£, = (A' + l)?—±x A '- 1 / 2 L\\ogxL 1 / A '\ (4.63) 



(xL 1 ^')-*' 
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Using Lemma [H the scaling expressions for Ph,*d and given in (IA.7j) . and taking 
the continuum limit, gives, after a substantial amount of algebra, 

Q(e;A)<r 1/2 ^A), (4.64) 

where 

F {v) J l + ^ V ^°°> (4.65) 

[cv^ 2 - a ^y/\hgv 2 \+0{v A '), v^O. 

The main result of this section is 

Theorem 1 The comb with the power law measure (I4.48P for the tooth length has 

d° s = 1, d? = 2 - -. (4.66) 
The result follows immediately from ffl~25j) . f l4~29j) . ffl~5lj) and f l4~64j) . 



Multiple Scales 

Given the results for the power law distribution it is natural to investigate the behaviour 
for a random comb that has a hierarchy of length scales. The easiest way to achieve such 
a comb is through a double power law distribution, 

f 1 - L7 l - Lr\ ^ = 0, 

Ai (£; U) = 1 2 (4.67) 

We may assume without loss of generality that the length scales Li scale in the continuum 
limit to lengths A» such that Ai < A2 and that 1 < a» < 2. 

Following the procedure of previous sections a lower bound on <3(£;Aj) is obtained 
by using Lemma [6] and noting that x{ x 'i Li) for this comb is essentially the sum of the 
cumulative probability functions for each power law. This gives 

1 - P T (x; Li) < rrUx)^ V (^. moo ( x )^ + b -^ moD ( x )^ + h f 

■ , V l^i Li Li 

i=l x 

(4.68) 



4.4 Multiple Scales 



51 



Choosing Lj to scale like Lj = a~ *\ \ where = 1 — ajj/2, gives a bound on the 
continuum return generating function of, 

T 1/2 (co + PiCeAx)-^- 01 ^ + c 2 (Z\ 2 y^^)- 1/2 < Q(£). (4.69) 

An upper-bound on Q(£) is obtained by application of Lemma [8] in which we set H = \_H\ , 
D = [D\ and k — \k], where 



H 
D 

k 



-1/2 



/3x~ 1 ^ 2 G(xL]^ A ' 1 , xL 1 ^' 2 )' 1 1 log xL 1 ^' 1 1 



(4.70) 



and for convenience we have introduced the function, 



G(vi,i> 2 ) 



«! — 1 «2 — 1 

Using Lemma [8] and the scaling expressions in I A. II then gives 



g(x; Ao<r 1/2 



i-(i-v~ s T + 



3coseclT(| logt£|) - 7 + W7 2 + 1 + 2 7 coth(| logufl/G) 



(4.71) 



(4.72) 



where t>j = £Aj, s = sgn(logt>i), 7 = tanh(l) and we have suppressed the arguments of 
G(vi,v 2 ) in order to maintain readability. 

We can now examine ( I4.69P and f)4.72p to see what they tell us about the behaviour 
of Aj) on various length scales. 

• When £ ^> A^ 1 both upper and lower bounds of Q(£; Aj) are dominated by the 
behaviour so taking the £ — > 00 limit leads to d° = 1 as in the previous sections. 

• If oti < a 2 then when £ <C A^ 1 both upper and lower bounds of Aj) are dom- 
inated by the £~ ai / 4 behaviour so taking the £ — >■ limit leads to = 2 — a;x/2. 
There is no regime in which a 2 controls the behaviour. 

• If a 2 < ol\ then when £ <C hr l where 



A- 1 = X ( { 



(2— oi)/(ai— 02) \^(2-02)/(o2-oi) 



(4.73) 
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both upper and lower bounds of Q(£] Aj) are dominated by the £~ a2 / 4 behaviour so 
taking the £ — > limit leads to df = 2 — 0:2/2. However there is an intermediate 
regime A -1 <C £ <C A^ 1 where the £~ ai / 4 behaviour dominates and Q(£; Aj) lies in 
the envelope given by 

Cl r Qi/4 < m \) < c 2 r ai/4 V\^i¥\, (4.74) 

where the upper and lower bounds will have corrections suppressed by powers of 
£A2 and the upper bound will also have corrections of order Both A2 and j3 may 
be chosen to make the corrections arbitrarily small in this scale range. The system 
therefore appears to have spectral dimension 5s = 2 — ai/2 in this regime. This 
is a fairly weak statement because Q(£) could in principle exhibit a wide variety 
of behaviours between its upper and lower bounds; this region is just a part of the 
crossover regime from d s to d™. However, as we are free to chose A 2 to be as large 
as we like compared to Ai, this regime can exist over a scale range of arbitrarily 
large size. We therefore can force the leading behaviour of Aj) in this range to 
be as close to a power law with exponent 5 s = 2 — a±/2 as we like. This is what 
might be observed, for example, in a numerical simulation; if the difference between 
the scales Ai and A 2 is large then there will be a substantial range of walk lengths in 
which the data will indicate a spectral dimension of 5s- We will refer to a spectral 
dimension that appears in this weaker way as an apparent spectral dimension and 
denote it by 5s rather than ds- 

Generic Distributions 

So far we have considered combs in which the distribution of tooth lengths has been 
governed by power laws or double power laws. In this section we extend the results of 
the previous sections to the case were the form of the tooth length distribution is left 
arbitrary. The most general situation is that the measure on the combs is a continuous 
function of some parameters u>j. The continuum limit of such a comb is obtained in the 
usual way but with parameters Wi scaling in a non-trivial way; Wi = w Ci + a di Ui- Given a 
random comb with such a measure we would like to know how many distinct continuum 
limits exist and for each compute how the spectral dimension depends on the length scale. 

The approach we adopt here closely mimics the arguments of the preceding sections, 
indeed the main complication is technical. As we have seen the properties of the continuum 
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comb are controlled by the asymptotic expansion of x{ u ) as u goes to infinity The 
main difference in the generic case is that we may arrange matters so that the scaling 
dimensions of the coefficients in the asymptotic expansion are such that sub-leading terms 
appear in the continuum. For the generic case we obviously have no way of knowing the 
full asymptotic expansion of x{ u 'i w i)- However, we will see that for a large class of 
measures, the form of the asymptotic expansion is encoded in the asymptotic expansion 
of a particular generating function for V0j). 

Our first task is to introduce this generating function and relate it to the asymptotic 
expansion of x{ u \ w i)- To this end we introduce the notion of a smoothed sum [13J, 

oo oo 

X±(u;Wi) = ^2fi{£;wi)r]±(e/u) = fi{0;Wi) + Wi)r] ± (£/u) = /x(0) + x±\u;w i ), 

£=0 i=\ 

(4.75) 

where r\± is the smooth cut-off function introduced in IA.2l and u controls where the cut-off 
occurs. Such smoothed sums are related to x{ u \ w i) by 

X-(u; < x(u; Wi) < x+(u; w { ). (4.76) 

The reason for introducing the smoothed sums is that we may use powerful techniques 
from complex analysis to compute their asymptotic expansion (see eg [2]). The gener- 
ating function to which the asymptotic expansion of the smoothed sums is related is the 
Dirichlet series generating function of 

e=i 1 

We now introduce a number of results and notations, 

1. For a strip in the complex plane £(&, a) = {z : a < Ke[z] < b} where b > a, we say 
V^(s;Wi) has slow growth in E(6, a) if for all s G £ we have V^(s;Wi) ~ 0(s r ) for 
some r > as Im[s] — > ±oo. We say V^s; wi) has weak slow growth if the above 
property only holds for a countable number of horizontal regions across the strip. 
See figure PI 

2. Define S± to be the set containing the triples (— o + ir, —k, —r±) such that the 
Laurent expansion of V fJj (s;Wi)J^i[^/± t ](s + l)/s about the point —o + ir contains 
the term r±/(s + o — ir) k where k > and — a + ir / 0. Here Ai denotes the 
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Figure 4.2: An illustration of the strip £ in which P M (s;Wj) satisfies property (1) given 
given in the text. The poles of T>^(s\ Wi) are denoted by crosses and the horizontal strips 
in which the growth condition holds are indicated by the dark grey regions. 

Mellin transformation and ^-i-,, is introduced in IA.2I We will often find it useful to 
refer to the elements of S± using an index j and denote the jth element of S± by 
(— (Tj + tTj, —kj, —rj t ±). Note that since A^[^± e ](s + 1) is analytic in £ then the 
positions of the poles are determined only by V^s; Wi). 

3. Define the indexing sets Jr and Jc, such that if j G Jr then Tj = whereas if 
j G Jc then Tj ^ and for both — <jj + irj G £ i.e. they index the poles in E which 
lie on the real line and off the real line respectively. 

4. For a Dirichlet series with positive coefficients the abscissa of absolute and condi- 
tional convergence coincide. Furthermore, since ^^ / x (^' u '«) = 1 ^ ne abscissa of 
convergence is less than zero. 

5. Landau's theorem: For a Dirichlet series with positive coefficients there exists a pole 
at the abscissa of convergence. A corollary of this is that the coefficient of the most 
singular term in the Laurent expansion about the abscissa of convergence is positive 
and furthermore it is the right-most pole of the Dirichlet series in the complex plane. 

If T)^ is of slow growth in E(a, b) and no pole occurs to the right of E(a, b), then by 
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using the expression (1A.26j) the asymptotic expansion of x±{ u '-> w i) is> 

X± { u \ w i) = M°5 W i) + 77— / U S VJs\ Wi) M ^ ±e ^ S + L rfg _ Wi 

Z7TZ 



c 



(logu) 



where C is the rectangular contour introduced in lA.slpl . rj ± are the coefficients of the 
Laurent expansions of V^s; Wi)Ai[^/± e ](s + l)/s, and R{u) is a remainder function which 
satisfies R{u) ~ M _iV , where iV > <jj, as a goes to infinity. Note that the difference 
between the coefficients r Jj+ and r 3 - is of order e, the parameter introduced in IA.2[ which 
can be taken to be arbitrarily small and so r Ji+ and are for all purposes equal. 

We now are in a position to prove the main result of this section. Note that in the 
following we require the continuum comb is such that A M = 1/2 and so has spectral 
dimension one on the smallest scales. It seems very likely that this restriction could be 
lifted but we will not pursue such a generalisation here. 

Theorem 2 For a comb in which the teeth are distributed according to the measure /i 
then ifT>^ has slow growth in the strip £ = {z : — 1 < Re[z] < 0} and has no poles on its 
line of convergence besides at the abscissa then continuum limits of the comb exist with 
d® — 1 and df taking values in the set {"^p- '■ < o~j < l,j G Jr}- 

Proof The proof proceeds in much the same way as the proofs in the previous sections; 
we first derive upper and lower bounds on Qc and then use these bounds to deduce the 
behaviour of Qc in different scale ranges. 

We begin with the lower bound which by Jensen's inequality amounts to finding an 
upper bound on 1 — Py. Note that (I4.3ip implies that the scaling dimension of 1 — Pt 
must be greater than or equal to one in order for A M = 1/2 and only contributes to the 
continuum limit if it has scaling dimension one. From Lemma [7] we have, 



« 1 

— / moo(x) 

I ~ P T {x]Wi) <m 00 (x)y/x I (1 - X-W Wi))du 

Jo 





I mac (x) 

< y/xm^{x)C{K]Wi) +m o0 {x)\fx j (1 - X-W Wi))du, (4.79) 

J K 



5 Since we have assumed there are no poles to the right of S we may choose c = b, where c is the 
constant appearing in the definition of the contour. 
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where C(K) = J (1 — X-i u i Wi))du is a constant independent of x. It is important to 
recall that since < x±( M ; w i) — 1 f° r au values of u and wi, . . . , %, in particular when 
Wi, . . . ,Wm assume their critical values, that the scaling dimensions of the coefficients rj t ± 
appearing in (14.781) must be positive as otherwise x±( u 'i w i) would diverge if we were to 
set Wi, . . . ,wm to their critical values. Upon performing the integration in ( 14. 79ft we will 
find that a given r^± now is the coefficient for a number of terms of increasing scaling 
dimension. Since we require the scaling dimension of 1 — Pt to be greater than or equal to 
one, only the term with smallest scaling dimension can appear in the continuum limit and 
we will drop any term that does not appear in the continuum limit. If we now substitute 
(SITED into (14791) we find 

1 - Pt(x; Wi) < v ^m 00 (x)C(K; «,,) + ^ V (-'og^" 1 ^ + 




— log m r - 1 " 



\fc,— i 

oo J 



(kj - 1)! 



m%, (4.80) 



where the remainder term R disappears since it cannot appear in the scaling limit. We 
now suppose we may choose the critical values and scaling dimensions of the parameters 
Wi such that Tj t ± has a scaling form that can be written as, 

rj,±= (^-) ' (-h°Z a ) ' , ( 4 - 81 ) 



where 9j, 9j and Xj t ± are constants and we have included the factor of | for later con- 
venience. From ( 14.80(1 one can see that since we require A M = 1/2 the only terms which 
appear in the continuum limit are those for which 9j = (1 — (Tj)/2 and 9j = —(kj — 1) and 
so it is useful to define a restricted indexing set 

jR = {jeJ R :9 j = {l-a j )/2 and 9 d = -(% - 1)} (4.82) 

and an equivalent one Jq for Jq- The continuum limit is thus, 

1 - Pt(x; Wi ) <aCZ+J2 {k ._ 1) \ 1 _ a ^ K V)"^ + ( 4 - 83 ) 
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where 4>j = (\\j-\/\j-)( 1 ~ aj ^ 2 (l — Oj + ir^)~ x . The lower bound on Qc is then obtained 
from Lemma [6l 

Before deriving an upper bound on Qc, we must analyse the consequences of any 
of the oscillatory terms in the second sum appearing in the continuum limit (i.e. if Jc 
is non-empty). In the case of the double power law measure we saw that one could 
obtain intermediate behaviour in which an effective spectral dimension was measured 
that differed from the UV and IR spectral dimensions. A similar phenomenon occurs in 
the generic case when the various length scales of the continuum limit are well separated. 
If we scale the coefficients of the oscillatory terms such that these terms appear in the 
continuum limit then we are lead to an inconsistency by the following argument, 

1. Let the term associated with an oscillatory term have index jo- Consider choosing 
the scaling form of x and Tj- such that £|Aj 0i _| <C 1 but ^lAj,-! 1 if |Aj,_| ^> 



2. If we were to take the scaling limit in such a scenario then the term that dominates 
the size of 1 — Pt is the one associated with the length scale |A JOi _|. This term will 
be oscillating around a mean of zero and hence must be going negative infinitely 
often as we approach the continuum. 

3. Since Pt is a probability generating function it cannot be negative, hence showing 
we have an unphysical limit. 

The cause of this behaviour is that we declared by fiat that the parameters Wi must be 
scaled to give (I4.81P ; however the parameters must also satisfy the constraints Wi) > 
and //(^; Wi) = 1 for all values of the parameters and we have not ensured these 
constraints are compatible with (14.811) . It is sufficient for our purposes to understand that 
the oscillatory terms prevent the continuum limit being taken for certain walk lengths so 
they are certainly unphysical; we therefore must scale them so that they disappear in the 
continuum. We therefore have as a lower bound on Qc, 



We now consider the upper bound. Without loss of generality we may arrange that 





(4.84) 
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the indexing set Jr has the property that if ji < then Xj lt ± < Xj 2i ±. Furthermore define 
H(x) = [H(x)\ = [aT 1/2 ], 

D{x) = [D(x)} = [^^(1 - x (^(x)))- 1 |log(r / , + a;^ i (-iloga;) fc - 1 )|], (4.85) 
k(x) = [k(x)] = [x- 1 / 2 (l- X (H(x)))], 

where I = inf Jr, i.e. \i t ± is the smallest length scale with a spectral dimension differing 
from one. Some modified versions of the above quantities will also be needed, 

D ± {x) = [D ± {x)\ = [^^(l- x± (if( a; )))- 1 |log(r /i+ x^ i (-ilogx)^- 1 )|l4.86) 

k ± (x) = [k ± (x)] = [x-^ 2 (l-x±(H(x)))}. 

It is clear that D^(x) < D(x) < D + (x) and k + (x) < k(x) < k_(x) which together with 
Lemma [8] allows us to conclude that the continuum limit is, 



Q(x;ui) < £ 



-1/2 



1 _ e G-io g (i-vy^) + 



,(G+-l)Iog(l-«7"') 

= — ^ 



3cosech 2 (|log^|(l- 1/G+)) - 7 + • v / 7 2 + l + 2 7 coth(|log^|/G' + )J 



J7) 



where we have defined G±(£, X±, . . .) = J2jeJ R v j,± 2 ~ -0' anc ^ arranged that no 
oscillatory terms appear. 

We are now in a position to prove Theorem [2] by analysing the behaviour of (I4.84p and 
(14.87P for various walk lengths. We first note that both the upper and lower bounds are 
controlled by the relative sizes of the quantities t>- 1+CTj//2 . In particular the largest of these 
quantities, T = sup{v- 1+a ^ 2 : j £ J i? } ) will determine the behaviour in a particular scale 
range; if we suppose T = vj 1+ * i/2 for some j then the leading contribution to both the 
upper and lower bound will be proportional to £ - ( 1+ °j)/ 4 . Using this we find the following 
behaviour, 

• On very short scales corresponding to £ ^> Xj 1 for all j then the lower bound of 0(0 
is dominated by the behaviour, which together with the trivial upper bound 
means that taking the £ — > oo limit leads to = 1 as in the previous sections. 

• On very long scales corresponding to £ <C Xj 1 for all j then there will exist a length 
scale A such that for < £ < A -1 , T = v~ + J ^ where Oj < <x,- for all j. Explicitly 
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A will be given by, 

A = sup{A A, :j6 Jr}. (4.88) 

j j 

If we take the limit £ — > we obtain df = (3 — er--)/2. 

We see that the spectral dimension increases monotonically on successive length scales 
and the spectral dimension has measured values given by d s = (3 — er--)/2 therefore proving 
Theorem [2l 

It is interesting to note that a constraint on the crossover behaviour exists for generic 
measures much as it did for the case of the double power law measure. In particular on 
intermediate scales there exists J G Jr such that £ <C A" 1 for j < J and £ ^> Xj 1 for 
j > J. Hence v, l+a ^ 2 <g i if j > J anc l so T = sup-fi^- 1+a ^ 2 ■ j g j Rj j < j}. This 

means there will exist a length scale Aj such that for Aj + i _1 < £ < A j^ 1 , T = v- * j(J) ^ 
where we have defined j(J) implicitly by ffjrn < &j for all j < J. Of course this does not 
ensure that Aj + i _1 < Aj -1 , indeed if the length scales Aj are not sufficiently separated 
this may not be true and we would not have a scale range in which this term dominated. 
The expression for Aj is, 

Aj = supjA^ A~. (J) n ' : j e Jr,J < J} (4.89) 

and so we may choose Aj + i independent of Aj thereby allowing the scale range over which 
this behaviour exists to be arbitrarily large. This would result in an apparent spectral 
dimension of 5 S = (3 — cr-.,n)/2. 

Finally, an interesting application of the techniques used to prove Theorem [2] is that 
they allow one to analyse a wider class of combs than the class for which the results of [3] 
are valid. In particular it was proven in [3] that for a random comb the spectral dimension 
is d s = (3 — 7o)/2 where 7 = sup{7 > : J 7 < oo} and J 7 = X]^o ^ ^ ■ This was P rove d 
subject to the assumption that there exists d > such that, 



J2 M*) ~ xd ( 4 - 9 °) 



-1/21 



as x goes to zero. 

Given the results in this section we see that —70 may be interpreted as the abscissa of 
convergence for the Dirichlet series generating function. Furthermore, it is clear from our 
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results that there are distributions where the assumption (14.901) does not hold and that 
we may use the techniques we have developed to analyse these cases. Recalling that for 
a random comb we may compute the spectral dimension using the relation (14. 18j) then 
we must perform a similar analysis to that done for the continuum comb but now only 
scaling x to zero. 

Due to Landau's theorem there always exists a pole at the abscissa. Suppose it is of 
order k and consider 1 — x±{ u )) 



1 - X±(u) 



r o,± + Re[ r j,±] cos(— Tj logw) + Imfr^-i-] sin(— Tj log 
fc-i 



u 



(logw 



,fc-l 



-u 



-70 



= n(«- A ; K±» 



(logw) 



-u 



-70 



(4.91) 



where j runs over the poles on the line of convergence and Q(y; {fj t ±}) is implicitly defined 
in the second line. We have not included any poles of order less than k since these will 
not be leading order as y goes to zero and we have assumed the there are no poles of 
order greater than k as then the above quantity would go negative due to the oscillating 
terms. Applying the same argument as we did for the continuum comb we obtain the 
lower bound, 



Q{x) > 



d{k 



[0(^/2; ri)± (l_ 7o±iTj .)-l) 



1/2 



-I logs 

2 I 5 



(i-*)/a 



x 



7Q + 1 

4 



(4.92) 



where d is a constant. To obtain an upper bound we apply a very similar argument as 
used for the continuum comb; the only difference being that we choose, 



(4.93) 



- fix 1//2 | log a; |, 

--^{\-xiM)Y x \\ogx\, 

x-^il-xiH)), 
= f3(l- X ±(H)r 1 \logx\, 
= x- 1 ' 2 {l- X ±{H)) 

and compute the behaviour as x goes to zero (x, (D — l)(k — 1)), Ph*d+ an d the 

size of the set A in (14.581) . The result is, 
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D± 
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where c" is a constant. We see that we reproduce the result of [5] when k — 1 and there 
are no poles on the line of convergence. If k ^ 1, the spectral dimension is the same as 
the k — 1 case as only logarithmic corrections are introduced. If we allow poles to appear 
on the line of convergence then Q(x; Tj^±) will have an oscillating x dependence which, if 
the bounds above are tight enough, would also imply that the functional form of Q is not 
a power law and so the spectral dimension does not exist. Whether the above bounds are 
tight enough to make this conclusion we leave for future work. 

Discussion 

In this chapter we demonstrated that there exist models in which a CDT-like scale de- 
pendent spectral dimension could be shown to exist analytically. That we could do this 
was important, as all evidence for a scale dependent spectral dimension in CDT has been 
numerical in nature and therefore does not provide any understanding of the mechanism 
causing the reduction in dimensionality. We hope the work begun in this chapter may be 
extended to shed some light on this question. 

We have given in Theorem [2] a reasonably complete classification of what behaviours 
one can find on a continuum comb and indeed it is likely that the cases not covered 
by the theorem do not have a well defined spectral dimension. The behaviour for the 
combs covered by Theorem [2] is fairly rich as there exists a cross-over region between the 
UV and IR behaviours in which a hierarchy of apparent spectral dimensions exist. One 
might expect that if one were able to extend the results of this chapter to more CDT like 
models one could find examples in which the spectral dimension again depends on the 
scale. Such an expectation is not unreasonable since the proof of the bijection between 
two-dimensional CDT and trees in [16] shows that any ensemble of trees is in bijection 
with a CDT-like theory. It is only when the random tree has the uniform measure that 
we obtain precisely CDT on the other side of the bijection. One could therefore imagine 
a variation of CDT, constructed from an ensemble of random trees with a measure that 
differs from the uniform measure and with a dependence on a length scale. This length 
scale could then be scaled while taking a continuum limit, as we have done in the work 
here. Such a model should display a spectral dimension that differs in the UV and IR. 

One may query how closely related a random comb is to an actual model of quantum 
gravity. For the case of both CDT and random trees one may show that the evolution of 
the spatial volume with time is generated by a Hamiltonian. A crucial difference when 
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one considers random combs is that the growth of a comb is not a Markovian process 
§; the probability of having a certain number of vertices at height h from the root may 
not be computed from the state of the system at height h — 1, meaning a Hamiltonian 
formulation of the evolution of a comb-like universe is not possible. An extension of the 
work described in this chapter, using the techniques of [15], to the case of random trees 
would be very interesting as this would constitute an example of CDT-like dimensional 
reduction for a model of quantum gravity which admits a Hamiltonian formulation. As 
was briefly mentioned earlier, a random graph that is even more closely related to the 
random graphs of CDT was introduced in [16] by identifying all verticies belong to the 
same spatial slice. This leads to a line like graph with multiple links between each node. 
The question of whether the above phenomenon can be extended to this case would also 
be interesting to pursue. 



6 In fact the growth process of a comb is Markovian if the tooth length distribution is of exponential 
form, however we have seen this does not lead to any interesting scale dependent behaviour of the spectral 
dimension. 



CHAPTER 5 



DT with Matter, String Theory and Branes 



In the previous chapters we considered the emergent dimension of spacetime in models 
of pure quantum gravity in addition to the disc function and the string susceptibility. 
This exhausts the list of observables that can be studied in such theories. We expect the 
theory to become richer upon the addition of matter and it is this topic we will address 
in this chapter. We will also discuss the relation of the metric-is-fundamental approach, 
including matter, to string theory. In fact the consideration of the string interpretation of 
these models will lead to a large class of two dimensional quantum gravity models known 
as (p, q) minimal gravity or (p, q) minimal string theory. 



Adding Matter using Matrix Models 

Of the approaches to discrete two dimensional quantum gravity considered in the preced- 
ing chapters the one which most easily admits a generalisation to models that describe 
the interaction of matter with gravity is the matrix model approach 0. An immediate 
consequence of this is that since the standard method used to compute amplitudes in 
CDT is the combinatorial method, it has so far not been possible to analytically investi- 
gate the consequences of adding matter to CDT. There have been a number of attempts 

^his does not mean the matrix model approach allows any matter to be coupled. In general it is 
restricted to CFTs arising from spin systems such as the Ising model. Furthermore, transfer matrix 
techniques have not been pursued to the same extent as matrix models in quantum gravity and therefore 
could yet yield a description of matter interacting with gravity. 
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to formulate CDT as a matrix model [551 [32J \5§\, however these have so far been too 
complicated to solve or resisted the generalisations necessary to include matter. Hence, 
for the rest of this thesis we will work within the DT approach. 

We begin by considering a naive method by which one can construct a matrix model 
that describes matter interacting with two dimensional gravity. In particular we know 
that the Ising model on a fixed lattice possesses a critical point at which the theory is 
described by a continuum quantum field theory. Therefore if we attach a spin to each 
triangular face in a triangulation and take a scaling limit in which both the lattice and 
the Ising model on it become critical, we expect we might recover the critical Ising model 
interacting with two dimensional gravity. For such a model, we expect the partition 
function for the discrete theory to be of the form, 

Z^= e-^> J ^- AA(T \ (5.1) 

TeTri h 

where J is the coupling between spins and <7j is the spin on the ith triangle in the tri- 
angulation, which can take on values of ±1. For purposes of illustration and to make 
contact with graph theory terminology it is usual to refer to the spins as "colours" so 
that an assignment of spins to a triangulation becomes a "colouring" of the graph. For 
a triangulation T, let n + _ (T) denote the number of edges shared by triangles differing in 
colour and let n p (T) be the total number of edges. The action may then be written as, 

J(n+_(T) - [n p (T) - n+_(T)]) - A|T| = 2Jn + _(T) - (A + 3J/2)\T\, (5.2) 

which implies that the partition function, assuming the couplings are small, may be 
written as, 

oo oo 

^ DT = EE^ m ^ fcm > ( 5 - 3 ) 

71=0 m=0 

where k = exp (—2 J), g = exp(— A — 3J/2) and A4(n; m) is the number of triangulation 
containing n triangles and m edges separating triangles of differing colours. It is this 
partition function we must realise using a matrix model. 

The key observation is that we can generate labelled graphs (i.e. containing colours) 
using a matrix model by considering models which contain more than one matrix. A 
matrix model containing more than one matrix can produce Feynman diagrams containing 
more than one type of "particle" species, which can then be considered as a label. As an 
example, consider the matrix model defined by, 

Z = j dM x dM 2 exp ^-JVTr (^MM + ^(M 2 + M 2 2 ) + ^(Mj 5 + M 2 3 )^ j . (5.4) 
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It is straightforward to see that this model will produce Feynman diagrams containing 
propagators for the Mi and M 2 matrices, which will be denoted by grey and white double 
lines respectively. Furthermore, due to the first term in the action, the model also pro- 
duces a mixed propagator that interpolates between grey and white. Finally, the cubic 
interaction terms will produce grey and white 3- vertex interactions. The various propaga- 



tors and vertices are shown in fig 5.1(a) together with their Feynman rules. An example 



of a Feynman diagram produced by this model is shown in fig 5.1(b) 

To simplify the algebra let the contribution of a colour changing propagator be 7 and 
that of a colour preserving propagator be S. The value of a Feynman diagram composed 
of n triangles with n + _ colour changing propagators and n p propagators in total is, 

where Q is the symmetry factor of the diagram. We can now compute the free energy 
for the matrix model by summing over all connected graphs and for each graph summing 
over If we let gS 3 ^ 2 = g and | = k then we obtain precisely (15. 3ft . 

5.1.1 1 The Continuum Limit and Multicritical Points 

The matrix model H5.4I) belongs to a general class of two matrix models of the form, 

Z = J dM 1 dM 2 exp (-iVTr ( y -cM 1 M 2 + V X {M X ) + V 2 (M 2 ))) . (5.6) 

These models have been solved [21] using loop equations yielding an equation for the 
resolvent W$](x) = (Tr(x - Mi) -1 ); 

(Vl(x) - y)(V 2 \y) -x)- P{x, y) + 1 = 0, (5.7) 
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where y = Vi(x) — (x) and P(x,y) is a polynomial in x and y with unknown co- 
efficients. The equation is sometimes known as the master equation. The equation for 
Wm (y) can be obtained by letting x = V^y) — (y) in the master equation. If we 
specialise to the Ising matrix model considered above we see that the loop equation re- 
sulting from (I5.7P is a 3rd order polynomial in W^(x). In order to compute the unknown 
constants appearing in P we enforce the condition that the master equation must be a 
genus zero curve. We may then compute the unknown constants by requiring that the 
resulting solution for the resolvent has the correct behaviour as x — > oo, i.e. there is a 
sheet of the resolvent in which W^ (x) ~ 1/x. In order to compute the critical point 
of the matrix model we again consider the quantity Wi(g,c) and look for the value of g 
at which w\ is non-analytic. In this case however we see that since we have a second 
variable, c, the matrix model will possess a critical line in g, c space rather than having a 
single critical point. In the presence of a critical line there is the possibility that there is 
a higher order critical point lying on the line. This is in fact the case for the Ising matrix 
model. 

We may qualitatively understand these critical lines and points by first letting c = 0; 
in this case the Ising spins do not interact and we have two copies of the one matrix model 
considered previously. We therefore expect that any critical point obtained by tuning g 
so that we hit the critical line, corresponds to pure gravity. At the critical point however, 
c has a very particular value and therefore we expect it is at this point that we obtain a 
theory of critical Ising matter interacting with gravity. Such higher order critical points 
are known as multi-critical points. 

One may wonder what the two matrix model with general potentials corresponds to. 
In this case there will be a number of free parameters contained in each potential and 
therefore we have the possibility of a large number of critical hypersurfaces of decreasing 
dimension. By tuning the coupling so that the model lies on some critical hypersurface 
we can obtain a variety of continuum models [121 153] . We shall see that in this way we 
may obtain a family of different field theories known as (p, q) minimal models interacting 
with gravity. The Ising model corresponds in this classification to the (3, 4) model. 

Conformal Field Theory and Minimal Models 

One manner in which one could specify a particular quantum field theory is to give the 
symmetries of the theory. One can then specify the dynamical objects in the theory, by 
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which we mean we specify which representation of the symmetry group they form. These 
two pieces of information largely fix the dynamics of the theory. In the usual formulation 
of quantum field theory these properties are introduced using the action and Feynman 
path integral, which can then be used to construct a perturbative series that computes 
correlation functions. However, in two dimensions the group of conformal transformations 
is infinite dimensional. If we have a QFT invariant under conformal transformations, this 
allows correlation functions to be computed without the use of an action or path integral. 
We will now review this construction following [60] . 

It is often customary when discussing two dimensional field theories to use complex 
coordinates, defined by z = x + iy and z = x — iy. A field <p(z, z) is defined to be primary 
if under a conformal transformation, z — > w — w(z), z — >• w = w(z), it transforms as, 

Fields which do not transform in this way are called secondary, of which an important 
example is the stress-energy tensor. The importance of the stress-energy tensor is that its 
Fourier coefficients generate conformal transformations. In terms of complex coordinates 
the stress energy tensor has two non-zero components that lie on the diagonal, we denote 
them T(z) and T(z). The Fourier expansion of the stress energy tensor may be written 
as T(z) = J2ne% z ~ n ~ 2 L n and T(z) = Tl, n& iZ n ~ 2 L n . We will refer to L n and L n as the 
generators of the conformal transformations. The algebra for the generators of conformal 
transformations is, 

[L n ,L m ] = (n 

[L n , An] = ( n 
L n , An] = 0, 

and is known as the Virasoro algebra, the real number c is known as the central charge. 

Since we work in a complexified Euclidean space the choice of how to slice the space 
into surfaces of constant time is arbitrary. This is equivalent to the freedom one has in 
splitting up a statistical mechanical system into subsystems in order to define a transfer 
matrix. It is useful in conformal field theory to take the surfaces of constant time to be 
concentric circles centred on the origin. The Hamiltonian for the system then generates 
the transformation z — > Xz, for A 6 R. The origin is taken to be at t = — oo and t = oo 
is at complex infinity. This setup is known as radial quantisation. Because the states at 



m)L n+m + 5 n+mfi —n{n 2 - 1), 
m)L n+m + 5 n+mfl — n{n 2 - 1), 



(5.9) 
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t = — oo are localised at the origin there is a one-to-one correspondence in conformal field 
theories between local operators and states known as the state-operator correspondence. 
Hence for every primary field there exists a state known as a primary state. The primary 
states are the highest weight states when constructing the highest weight representations 
of the Virasoro algebra. The states obtained from the primary state by use of the lowering 
ladder operators, L n and L n where n < 0, are referred to as descendant states and they 
in turn give rise to fields known as descendant fields. A primary field together with its 
descendants is known as a conformal family. 

A final important ingredient is the operator product algebra of the fields. This algebra 
arises by considering the correlation function of a product of local operators and studying 
its behaviour as two of the operators are brought towards the same point. In general the 
correlation function will diverge as the operators are brought together, however we can 
characterise this divergence by use of a Laurent like series in which the coefficients are 
other local operators. This series is known as the operator product expansion (OPE). 
Crucially which operators appear as coefficients in the OPE is independent of the other 
operators appearing in the correlation function, i.e., 



Since the OPE is independent of the other operators appearing in the correlation function, 
this relation is usually written as, 



The operator product algebra is defined as the collection of all OPEs of the fields in the 
theory. It turns out that the operator product algebra is entirely fixed by the value of the 
central charge together with the 3-point correlation functions of the theory. 

In general the OPE of two operators will include other local operators in the expansion 
besides the two original operators and their descendants. Generically this will mean that 
in order to have an operator product algebra that closes the theory must contain an 
infinite number of primary fields. However, for special values for the central charge there 
exist finite sets of primary fields for which the operator algebra closes. These special 
conformal field theories are known as minimal models. In particular, given two positive 
coprime integers p, q such that p < q then the conformal field theory with central charge 




(5.10) 




(5.11) 



n=— oo 



C = 



1-6 



(q-p) 2 



(5.12) 
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and whose primary fields have weights given by 

(qr-ps) 2 - (q-p) 2 
K ' s = 4p~q ' (5 3) 

where r and s are integers satisfying, 1 < r < p and 1 < s < q, is known as the (p, q) 
minimal model. 



Boundary States of Minimal Models 

To study minimal models in the presence of boundaries it is most useful to consider the 
case of the minimal model on a finite length cylinder. The boundaries of the cylinder carry 
boundary conditions we call a and (3. Using real coordinates x and t we will analyse the 
case when the boundaries of the cylinder are curves of constant t. A necessary condition 
that any theory must satisfy at a boundary is that there is no flow of energy or momentum 
across the boundary. In terms of the components this requirement corresponds to T tx = 
at the boundary. If we complexify the coordinates by changing to the coordinates z = t+ix 
and z — t — ix, then the condition T tx = becomes T(z) — T(z) = 0. Since the boundary 
is a curve of constant time, the system at the boundary should correspond to some 
state in the Hilbert space. We will refer to these boundary states as \a) and \/3). As 
was mentioned previously, calculations in conformal field theory are usually done in a 
coordinate system in which the slices on constant time form concentric circles centred on 

2-7T 

the origin. We can transform to such a coordinate system via the map w = e~^ z resulting 
in the cylinder being mapped to an annulus. The condition on the boundary states can 
then be written as, 

(L n - Z_ n ) \a) = 0. (5.14) 

The solutions to this equation are known as the Ishibashi states and have the explicit 
form, 

\j) = J 52\j;n)®\j;n) (5.15) 

n 

where the sum is over all states in the irreducible representation of the Virasoro algebra 
associated to the primary field j. Note that there is one Ishibashi state for each primary 
field in the theory. We therefore appear to have a space, with dimension equal to the 
number of primary fields in the theory, of boundary states. However, it turns out that 
there exist only a finite number of consistent boundary states, which are known as Cardy 
states and are in one-to-one correspondence with the primary fields in the theory. 
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String Theory 

The central hypothesis of string theory is that fundamental particles may be modelled 
as quantum strings §. The space in which the string is embedded is known as the target 
space. The simplest action for the string is known as the Polyakov action and takes the 
form, 



where X is the map embedding the string in the target space and g a b is the metric on the 
worldsheet, M^, which is a two-dimensional Riemannian manifold of genus h. We have 
taken both the worldsheet and target space to be Euclidean as then the corresponding 
quantum theory is well defined. Later we will show that the quantum theory is equivalent 
to a theory on a flat worldsheet and so may be Wick rotated back to a Lorentzian theory. 

The Polyakov action corresponds to a string whose excitations are purely bosonic and 
is therefore known as the bosonic string. The classical bosonic string has a number of 
symmetries that are manifest in the Polyakov action; 

Target space Poincare Symmetry: The action is invariant under the transforma- 
tion X^ t^ v X v + a M , where A M is a Lorentz transformation. 

Worldsheet Diffeomorphisms: The action is invariant under worldsheet diffeo- 
morphisms, a a — » f a (cr b ), where / is a smooth function. Under such transformations the 
embedding fields X behave as scalars. 

Weyl Symmetry: Finally the action is also invariant under Weyl rescalings of the 
metric, g a b{(j) — > e 2 ^ x ^g a b(cr), where <f>(x) is a smooth function of x. 

To quantise the string we integrate over the embedding fields X and the worldsheet 
metric g. We must only count each physically equivalent configuration once and it is 
therefore important to decide which symmetries of the action we consider gauge symme- 
tries. For now we will consider Weyl transformations and worldsheet diffeomorphism to 
be gauge symmetries. This has the consequence that we must divide by the volume of 
these symmetry groups in the path integral. The expression for the partition function is 
then, 



In order to perform this integral we must define the measures [dX] and [dg]. One might 
already note that the partition function for the string is that of the metric-is-fundamental 

2 This is unlikely to be the correct underlying principle but it is how the subject developed historically. 




(5.16) 




(5.17) 
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approach but with the addition of a set of scalar fields and a lack of cosmological constant. 

For now we will keep the measures in the integral as formal objects and instead work on 
gauge-fixing the integral. The diffeomorphism gauge symmetry allows g a b to be brought 
to the form e 2 ^g a b(Ti) where g depends only on a finite number of parameters Tj, known 
as the moduli. Furthermore, due to the extra Weyl symmetry in this theory the metric 
may be brought to the form g a b(ji). Inserting the Faddeev- Popov determinant necessary 
to achieve this change of variable we obtain, 

Z h = J dr J p|^|L j [db][dc] J [dx]e -sAxr 9 M-s s , M h ] i (5 . 18) 

where h and c are ghost fields. In order for the integrals over d(j) and <i£ to cancel the 
volume of the gauge group we need the integral, 

Z' h [g,4>,£] = J [db}[dc] J [^e-^^'^-^ 16 '^ 1 , (5.19) 

to be independent of and £. The action itself is invariant under Weyl transformations and 
diffeomorphisms as these are classical symmetries. Hence the only potential problems lie 
in the definition of the measures [dX], [db] and [dc\. If these measures have a dependence 
on <f) or £ then we will have a Weyl or diffeomorphism anomaly, respectively. In fact, 
it happens that the measures are invariant under diffeomorphisms but the case of Weyl 
transformations is more complicated and it turns out that there is indeed a Weyl anomaly 
present. The Weyl anomaly can be most easily seen by recalling that a field theory is 
invariant under conformal and hence Weyl transformations if the trace of the stress energy 
tensor vanishes i.e. (Tjf) = g^ v (T^) = 0. 

If one computes the trace of stress energy tensor for the theory defined by ( 15. 19ft one 
obtains (T") = — D ~ 26 R[g), where D is the number of scalar embedding fields present and 
R is the Riemann tensor for the worldsheet. At this point it is easy to see how the critical 
dimension of the bosonic string arises; if the number of dimensions D is twenty-six then 
the Weyl anomaly is not present and the integral over <f>{x) trivially cancels the remaining 
gauge group volume. This also has a number of important consequences. Firstly, since the 
geometry of the worldsheet is entirely a gauge degree of freedom we need not worry about 
which particular class of worldsheet geometries, for example Euclidean of Lorentzian, we 
sum over. Furthermore, since the worldsheet metric can be made flat everywhere, the 
computation of worldsheet correlation functions, which are important in the calculation 
of S-matrix elements, are reduced to correlation functions in a flat space conformal field 
theory. 
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Figure 5.1: The worldsheet corresponding to the interaction of two strings may be trans- 
formed via a conformal transformation to a worldsheet of equivalent topology but with 
all asymptotic string states transformed into local operators. 

Having defined the partition function for the string, an ansatz for how to compute S- 
matrix elements is required. If one considers a diagram for the interaction of two strings 
one can note that, by utilising the conformal symmetry of the worldsheet theory, the 
strings propagating from infinity may be encoded in local worldsheet operators, which we 
denote by V. Hence the string interaction diagrams may be reduced to the computation 
of correlation functions in the worldsheet theory. This is shown in fig 15.11 Furthermore, 
the genus of the worldsheet determines the power of the string coupling appearing in 
the amplitude. The ansatz is therefore that S-matrix elements for string states may be 
computed perturbatively by computing 



where (Vi . . . Vjv) ft is the iV-point correlation function on a surface of genus h and g s is the 
string coupling constant. Clearly this construction fails when the string coupling becomes 
of order one. In such cases one would need a non-perturbative definition of string theory 
in order to proceed. Unfortunately, such a formulation is currently incomplete. 

Although, a full definition of string theory is still lacking there has been much progress 
in understanding some of the crucial ingredients of such a theory. In particular it was 
found that Dirichlet boundary conditions of open strings correspond to the interaction of 
a string with a dynamical object which has become known as a D-brane. An important 
question in any string background then is what is the spectrum of D-branes, as this will 
have consequences for the non-perturbative physics. 

Another aspect in which the above formulation is incomplete is that the target space 
was taken to be flat Minkowski space of dimension D. Obviously there are other space- 



oo 




(5.20) 



h=0 
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times of great importance that arise as solutions to general relativity; namely the black 
hole spacetimes such as the Schwarschild metric and time dependent solutions, relevant 
to cosmology, such as the FRW metric. In particular it has proved difficult to gain any in- 
sight into the Big Bang singularity due to the difficulty of computing in a time dependent 
background. 

Finally, a point related to the problem of backgrounds is that in the most straightfor- 
ward interpretation given in this section, the number of dimensions in the target space is 
determined by the Weyl anomaly. However, the Weyl anomaly is not as stringent as this; 
it only requires that the overall central charge of the gauge fixed theory is zero. The ghost 
sector will always arise during the gauge fixing process and comes with central charge of 
—26, we therefore must start with a worldsheet theory with a central charge of 26. There 
are many ways one could arrange this besides requiring 26 non-interacting scalar fields 
to live on the worldsheet and indeed, one must modify the string model in this section 
as there are clearly not 26 non-compact dimensions! The most studied way to resolve 
this problem is to postulate that some of the target space dimensions are compact; in the 
case of the superstring these compact dimensions are usually taken to be a Calabi-Yau 
manifold. However there is nothing preventing us from cancelling the anomaly instead by 
introducing a conformal field theory on the world sheet with no geometric interpretation 
at all; this would be interpreted as degrees of freedom internal to the string. This will 
form the topic of the next section. 

.1 1 Liouville Theory and Non-Critical String Theory 

Usually the dimension of the target space is defined to be equal to the number of world- 
sheet scalar fields and we can therefore formulate string theory in a lower number of 
dimensions by incorporating a CFT with a non-geometric interpretation into the world- 
sheet theory. There is a well studied CFT which has a tuneable central charge satisfying 
the constraint c > 25; it is known as Liouville theory. The action for Liouville theory 
takes the form, 



where Q = b + 1/b and b and /x are free parameters. The field 0(a) transforms as 
0(er) — > (f)(cr) — u(a)/(2b) under the Weyl transform g — > exp(2u(a))g. More generally 
under a conformal transformation z = z(w), we have <fi(z) = (f)'(w) — ® log \ j^\ 2 - Finally, 



>S , Liouville[0, h] — — 




(5.21) 
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the metric appearing in this expression is defined to be Euclidean. This theory may be 
quantised, while preserving the classical conformal symmetry it possesses, to produce a 
theory with central charge c = 1 + 6Q 2 and primary fields of the form V a = e 2a< ^. 

If we include the Liouville theory in the worldsheet theory then we are left with a 
central charge of c < 1 and so we still require some other sector in order to completely 
cancel the Weyl anomaly. In the case when c = 1 a natural choice is to include a single 
scalar field. This is known as the c = 1 string and has been studied quite extensively, see 
[52] and references therein. However, if we choose b in ( I5.2ip such that c < 1 then instead 
of scalar fields, the natural candidate for the matter sector is a (p, q) minimal model. The 
resulting string theory is known as (p, q) minimal string theory. It is these string theories 
that will be studied in the remainder of this thesis. 

One might wonder, since the (p, q) minimal string has no embedding fields, if we have 
lost all geometric interpretation of the target space. However, there is a geometric way 
of interpreting the (p, q) strings. The (p, q) minimal models often arise as the continuum 
theory describing a statistical mechanical lattice model [53| |4"2] , such as the Ising or Potts 
model, at its critical point. Such lattice models assign a degree of freedom to each lattice 
point which can take values in a finite set G of points. If we interpret this degree of 
freedom as an embedding field then these model describe the embedding of the lattice 
into G. Hence we can think of the (p, q) minimal strings as being strings embedded in a 
discrete target space. 

We may go further and show that the Liouville sector of the worldsheet theory also 
has a geometric interpretation. Let us turn our attention back to ( I5.18P in which we 
have treated both the Weyl transformations and worldsheet diffeomorphisms as gauge 
transformations. Because of this it is important that there is no anomaly in either of 
these symmetries. However, it is possible to construct a different theory if we drop the 
interpretation of the Weyl transformations as a gauge symmetry. In this case there is 
nothing preventing the appearance of a cosmological constant term, resulting in, 

Z h = jdr J J [db] [dc] J [rfX ]e-^ [X ^-^ [6 ' c '^ ] -^ (5.22) 

In analysing this theory it will be useful to give a concrete definition of the measures in 
the integral, which we do by following (HU |62} [63]. To define a measure on a space of 
functions we must introduce a metric in that space. For the two dimensional spacetime 
metric g this metric is unique if we require it to respect diffeomorphism invariance and 
contain no dependence on derivatives of g; a condition known as ultralocality. The metric 
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then takes the form, 

(5g, 5g f ) = [ d 2 a^ {u5g ab 5g' cd g ab g cd + 5g ab 6g' cd g ac g bd ) , (5.23) 

J M 

where u is an arbitrary constant. Since this theory lacks a Weyl symmetry the diffeomor- 
phism gauge symmetry only allows g ab to be brought to the form e 2 ^^g a b(T,i). Inserting 
the Faddeev-Popov determinant for this transformation we obtain, 

Z h = J dr J l^pl J [^[dclfdXje-^^'^-^^^WM. (5.24) 

We can integrate over diffeomorphisms £ as both the integrand and the measures are 
invariant under them. However, as we have seen, the measure for the ghosts and embed- 
ding fields still contains a dependence on 0. Let us consider how this works explicitly. 
The metric on the space of embedding functions X is unique if we require it to respect 
diffeomorphism invariance and be ultralocal. The metric then takes the form, 

(6X,6X')= I d 2 a^5X(a)5X'(a), (5.25) 

JM h 

with a similar expression for the ghosts. As can be seen, both of these metrics have a 
dependence on 0. Using this definition for the measures it is possible to show that the 
dependence of the measures is given by, 

J[dXU~ g = J[dX]se& s M 

J [dbdc]^ = J [dbdc]- g e^ SL[4,] (5.26) 

where Sl is the Liouville action with Q = 1. We therefore obtain, 

Z h = J dr J [d<j>] J [db\[dc\ J [rfX]e- Sp[ ^^ s ^ t[6 ' c ' 5 '' ll - i ^ 5i[ ^ 5] , (5.27) 

where the original cosmological constant has been absorbed into the cosmological constant 
appearing in the Liouville action. Naively, we appear to have solved the problem of how 
to integrate over 0, however this conclusion is premature. The measure for <fi is inherited 
from the measure on g ab and therefore depends on itself. It has the form, 

(50,50')=/ rfVe^v^5050'. (5.28) 

JMh 
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The dependence on makes it difficult to split the measure into a product of Fourier 
coefficients which can easily be integrated over. In order to complete the calculation of 
Z we must relate this <fi dependent measure to one that is independent of <fi. It was 
conjectured in [64] that the effect of changing from this metric to a flat one, given by, 



would only renormalise the constants appearing in the Liouville action. This conjecture is 
borne out by comparison to the matrix model approaches. We therefore see that although 
classical gravity is entirely topological, the measure in the path integral generates dynam- 
ics for the metric. Furthermore, the Liouville action we added to the worldsheet at the 
beginning of this section can be understood as describing the degrees of freedom internal 
to the worldsheet that correspond to the scale factor of the metric. The consequence of 
this on the statistical lattice model picture of the minimal models is that the Liouville 
field describes the configuration of the lattice itself. We therefore see how the program of 
DT naturally emerges from the path integral. 

.2 | Wave-functions and Correlation functions in Liouville Theory 

We will now discuss an approximation, known as the mini-superspace § approximation 
that has proved useful in the study of Liouville theory. This approximation is obtained 
by neglecting all oscillations of the string, leaving the zero modes as the only dynamical 
quantity. Beginning with the Liouville action we may obtain the Hamiltonian and by 
setting the transverse oscillations to zero we obtain the mini-superspace Hamiltonian, 



where 4>q is the zero mode of 4> and p is its conjugate momentum. It will turn out to be 
useful to study this Hamiltonian using the Schrodinger representation in which states in 
the Hilbert space becomes functionals of the field variable. The energy eigenstates then 
satisfy, 



3 The term "superspace" in this context refers to the space of geometries of the worldsheet rather than 
a supersymmetric space. 




(5.29) 




(5.30) 
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An interpretation of the wave function for the zero mode may be found by consideration 
of the state-operator correspondence. Although discussed previously in Section I5.2[ the 
state-operator correspondence can be given a more explicit realisation using path integrals 
[49] . Working in the radial quantisation picture we can associate any local operator to a 
wave-functional by the following map, 



where the path integral is taken over the interior of the unit disc and the field <j)(x, t) 
assumes the value <fib( x , 1) a t the boundary. We therefore expect the mini-superspace 
wave functions to be approximations to disc functions in which the Liouville field has 
some specified value on the boundary. In order to compare the mini-superspace wave- 
functions to later quantities it is useful to note that in the mini-superspace approximation 
the value of (fib is completely fixed if we specify the boundary length of the disc. Indeed, 
the two are related by / = 27re ^ and so we may rewrite (15.311) to give, 



where we have now dropped the subscript b on the field (fi. The solutions to this equation 
are the probability distribution for the boundary length of the disc and they have the 
form, 



where A is a constant. Furthermore, given the state operator correspondence we expect 
the parameter E to label which local state has been inserted in the disc. It is this quantity 
we will compare with exact results from Liouville and matrix model calculations. 

We now turn to the problem of computing correlation functions in Liouville theory. 
As was mentioned in section 15.21 the only input necessary to fix the operator algebra of 
a conformal field theory is the three point function coefficients. We therefore need only 
compute the three point correlation function. It is unnecessary for our purposes to review 
the computation of the three point function in detail, however it will be useful to consider 
the general approach. Historically, the Liouville three point function was first computed 
in [191 [20] using the method of [69J, however a more efficient method was introduced later 
in [70]. The method of [69J is based on the observation that for particular combinations 




(5.32) 




(5.33) 




(5.34) 
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of operators and worldsheet topology the correlation function may be expressed in terms 
of free field correlators. Explicitly, consider the following n-point correlation function, 

/n 
[d^JJe^e-^M. (5.35) 

We proceed by splitting the integral into an integral over the zero mode 0o an d the 
non-zero modes 0, 

id^He^e-^M / d <p Q e< Q ^ 2 ^ a ^e-^°^ d2x ^ 2b \ (5.36) 
i=o 




where is the free field action obtained by setting /i — in Sl{4>\ and x is the 

Euler characteristic of the worldsheet. We will now focus our attention on the zero mode 
integration, 



4 / dt 



o 



(5.37) 



2b 

where s = (Qx — 2 XT=o a i)/(2b) is known as the KPZ exponent and we have introduced 
t = exp(260o). It is clear that this integral is divergent for s > and therefore must be 
regularised in those cases. The method employed in [69] is to compute the zero mode 
integral for Re[s] < and then analytically continue, giving the result, 



d 2 x^e 2b *) TTe W ) > ( 5 - 38 ) 






where the average is taken using the free action [</>]. We now make the observation 
that, if we arrange the operators and topology of the worldsheet to be such that s G Z + , 
then the remaining correlation function in f!5.38j) may be computed as it is a free field 
theory correlator. One problem we must address is the prefactor of r[— s] which is infinite 
at precisely the points for which the correlator becomes a free field correlator. In [20J this 
property was interpreted as implying that the three point coefficients have a pole precisely 
when the relation s G Z + is satisfied. The result (15.381) therefore gives the residue at these 
poles. Knowing the residue of the three point function at these poles allowed [20] and [19] 
to guess the correct analytic continuation of these results to all s. 

The above regular isat ion of the integral in the case s > was achieved via analytic 
continuation. A more brutal regularisation would be to simply cut off the 0o integral. 
Unusually, the divergence of the zero mode has both a UV and IR interpretation. With 
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respect to the background metric g the divergence of the zero mode integral looks like an 
infrared divergence which should be solved by placing the system in a finite sized box. 
Alternatively, with respect to the original metric g = e 2pr ^g, the divergence of the integral 
due to contributions from <fi — > — oo corresponds to a divergence due to the contributions 
from worldsheets in which all distances shrink to zero, which is a UV problem. We 
therefore may think of placing a cutoff at = — h log A as being either a short or long 
distance cutoff. If we implement the cutoff in (I5.37p . we obtain, 

1 y°° dt [r^e-^v^] = ^^-fi 3 d 2 x^ge 2b *y + P(ji,A), (5.39) 

where -P(/i, A) is a cutoff dependent polynomial in /i. The property that the cutoff depen- 
dent part of the above integral is polynomial in /i prompted [71] to classify contributions 
to correlation functions that are analytic in /i to be non- universal. Another way to think 
of this is that since these contribution arise from zero area worldsheets, then if they were 
to arise in the matrix model formulation of the theory such contributions would feel the 
effects of the lattice and hence be non-universal. Due to the exponential interaction term 
in the Liouville action the theory does not probe far into the region <fi > — ^ log fi. Hence, 
in the presence of the cutoff, the volume of the Liouville direction is V<f, = — ^ log (jjj- 

Recall that in the presence of boundaries we could compute amplitudes with either 
the boundary cosmological constant or the boundary length fixed, with these amplitudes 
related to one another via the Laplace transform (12.311) . One use of the fixed boundary 
length amplitudes is that one expects worldsheets of zero extent to not contribute since 
a surface can not disappear if the boundary length is non-zero|^. 

Contributions to the amplitude from degenerate worldsheets could only reappear pre- 
cisely when the boundary lengths go to zero. We therefore have two ways to determine 
which terms are non- universal; they are analytic in fi and also they should only have 
support when all boundary lengths are zero. 



5.4.3 | Minimal String Theory 

Minimal string theories are obtained by interpreting Liouville theory coupled to a (p, q) 
minimal model as the worldsheet theory for a string. The requirement that the total 

We make a distinction between vanishing worldsheets and one with zero area. Branched polymers, 
as generated by the Gaussian matrix model, represent an example of worldsheets with zero area but 
non-zero boundary length and therefore have universal meaning. See |25j for further discussion. 
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central charge of the system is zero, together with the expression for the central charge for 
the minimal models, Liouville theory and ghosts leads to b 2 = -. The physical operators 
of the (p, q) minimal string were first given by [28J as U n = O n V a ( n \ where, 



a(n) 

a(n) = b n^O mod p n^O mod q n G Z, (5.40) 

2p 

and O n is an operator containing contributions from the ghost and matter sector in 
addition to derivatives of the Liouville field. An important physical operator is obtained 
by gravitationally dressing the primary fields O r>s of the minimal model; such operators 
are referred to as tachyons and take the form, 

%, s = ccO rjS V ariS , (5.41) 

p + q-\qr-ps\ 
a r ,s = — , (5.42) 



pq 

where c and c are ghost fields. The expression for a TjS arises from the requirement that 
the conformal dimension of the total tachyon operator is one as this allows for the integral 
of the operator over the worldsheet to be conformally invariant. 



5.4.4 | FZZT Branes and the Seiberg-Shih Relation 

Let us return briefly to consider the standard critical string. It is well known that string 
theory contains other degrees of freedom which do not appear in the perturbative string 
expansion. The most well known of these are D-branes. Such objects correspond to 
extended membranes thereby generalising the idea of a string. The D-branes are partic- 
ularly interesting as they do admit a description in terms of perturbative string theory; 
although they do not appear in the closed string expansion they can be described by 
open strings in which the boundaries carry Dirichlet boundary conditions. If 26 — n of 
the string embedding coordinates carry Dirichlet boundary conditions then this means 
the end of the string is confined to move in a n dimensional hyperplane, thereby giving 
rise to the membrane like structure of the D-brane, which in this case would be called 
an n-brane. Since D-branes are described by worldsheets with boundaries the worldsheet 
theory is known as a boundary conformal field theory. 

Given this motivation it is natural to enquire whether the (p, q) minimal strings con- 
tain any branes. The answer to such a question would involve classifying the possible 
conformally invariant boundary conditions in the worldsheet CFT. A review of the result 
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of such a program for the (p, q) minimal model not coupled to gravity was given in Section 
15.31 We now consider a similar question for Liouville theory and hence the (p, q) minimal 
string. The Liouville theory on a manifold with boundary is defined by the Lagrangian, 

S L = / d 2 x^ (g ab d a <pd b <j) + Q<pR{g] + Air^) + / p B e b *dx, (5.43) 

4?r JM JdM 

where \x& is known as the boundary cosmological constant. The full conformal bootstrap 
of Liouville theory on the disc was performed in [18] [21], resulting in the identification 
of a consistent boundary condition for the Liouville theory; the FZZT brane. As part of 
the conformal bootstrap, the bulk one-point function on a disc, i.e. the bulk one-point 
function with FZZT boundary, was calculated: 

V ai (P) = (v P \a) FZZT = (vr /U7 (6 2 ))- lP/6 r(l + 2ibP)T{l + 2iP/b) COS ^ aP ^ ( 5 .44) 

iP 

where a is defined implicitly by, /icosh 2 7r6cr = /ie 2 sin(7r& 2 ) and \vp) is the primary state 
associated to Vq/2+ip- 

In the minimal string the full FZZT brane is \a; k,l) = \(t)fzzt ® \k,l), where \k,l) 
is the (k, I) Cardy boundary state of the minimal model. We will find it convenient to 
follow pT] in rescaling the bulk and boundary cosmological constant so that 717x7 (6 2 ) — > \i 
and HB\Z^l(b 2 ) sin(7r6 2 ) — > /is- This modifies the relation between /i and a to /Ib — 
y / /Icosh(7r6cr). 

It was noted in [19] [20] that the results of the conformal bootstrap on the sphere 

i ~ 1 

are invariant under the duality transformation b — > t, /i — > jl = /i^ . This self duality 
is present in the bootstrap on the disc if the boundary cosmological constant transforms 
as, hb A-B) where ^| = cosh The transformed boundary condition is referred to as 
the dual FZZT brane and provides a physically equivalent description of the FZZT brane. 
Given that all computed amplitudes exhibit this duality it is in our opinion important 
for this symmetry to be maintained in any modifications of amplitudes, such as those we 
consider later. 

The amplitude in the minimal string for the insertion of a tachyon in the bulk of a 
disc is obtained using (I5.44p [T7] . 

(T r ,s\v,k,l) = A D sin ( — J sin ( — J T{2ba - b 2 )T{2b- 1 a - r 2 - 1) x 

f ttA ( irt\ j_ ( irt(r \ , 
C^_i I cos - Ui-i I cos - y.** cosh f -— 1 , (5.45) 
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where t = qr — ps, 2a = Q — -J=, A D is a constant independent of k, r, s and I and 
U n (x) is the nth Chebyshev polynomial of the second kind. This expression for the disc 
amplitudes motivated Seiberg and Shih to conjecture in (T7J that the construction just 
outlined over counts the number of FZZT branes and that, up to BRST null states, the 
following identification holds, 

M,l>= E E k + i^a.!)- MB) 

n=-(H),2m=-(I-l),2 V 

This relation, which we will refer to as the SS (for Seiberg-Shih) relation, was originally 
obtained by inspection of ( I5.45P but was later derived for discs using the ground ring in 
[23] . Essentially it states that there is only a single FZZT brane \a; 1, 1) and that all the 
others are related to it by complex shifts in the boundary cosmological constant which 
has somehow absorbed all the information about the matter boundary condition. 

It is worth pausing at this point to compare the above exact result for the disc function 
amplitude in the presence of FZZT boundary conditions, with the result obtained via the 
mini-superspace approximation. To compare the above result (15.451) to ( 15.34ft . we must 
take the inverse Laplace transform. This may be achieved by use of the identity, 

7rbc ° s(27r 2 P ; N ) = r d ±e~ M ^«) K2ip{Ml y (5.47) 
2Psinh(^) J I — KJ v 1 

Setting M = ^[Ji, we find the length distribution for the 1 point tachyon disc amplitude 
to be, 

CT x (T r , s \o-Xl) oc K t/p (yfjll), (5.48) 

hence identifying E/b = n/p, where n is the integer defining the operator in ( 15.401) . We see 
that the exact calculation of the Liouville disc function reproduces the mini-superspace 
result up to rescalings of / and the normalisation constant. This has in fact been observed 
to be a generic feature of the mini-superspace approximation in this setting; it will give the 
correct functional form for any amplitude. This implies that for general disc amplitudes 
we have, 

(U n \ a; k, I) oc cosh ( ] . (5.49) 
This is consistent with the application of (I5.44p with Q/2 + iP = a(n). 
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.5 | Cylinder Amplitudes and the Seiberg-Shih Relation 



The minimal string cylinder amplitude was first computed in [35] for the case of \a; a, 1) 
boundary conditions. The cylinder amplitude with arbitrary boundary conditions was 
given in [36J and we have reproduced the derivation in Appendix [B] A compact expression 
for the general cylinder amplitude is, 



Z(k,l;a 1 \r,s;a 2 ) 



2vr 2 



00 p cos(27rcr 1 P) cos(27RT2-P) sinh( 

dP 



2tt_P_\ 

'pql 



,2mP^ 



V2 pq J P 2 + e2 



sinh(27r6P) sinh(^) 



where e is an IR cutoff and, 

\ p (k,r) A 9 (Z,s) p-1 q—1 

*w*) = E E E E 

7j=|fc-r|+l,2p=|i-s|+l,2 a=l b=_( g _l) 



. 7Tt . 7Tt 

sm( — ) sinl — ) x 



f/^_i(cos — )U p -i(cos- 



k,l,r,s\ 



/pq 
(5.50) 



(5.51) 



P 



TXt 



q cos — — cos z 

pq 



where X p (k, r) = min(A; + r — 1, 2p — 1 — k — r) and t = qa + pb. It is the purpose of 
this section to investigate whether the SS relation is satisfied when applied to cylinder 
amplitudes. For this purpose it is useful to define the deviation from the SS relations for 
a particular amplitude A(k, I; cr\X) as 

fc-i i-i 

. mp + nq 



AA(k,l;a\X) = A(k,l;a\X) - E E Ml,l\(T + 

n=-(fc-l),2 m=-(Z-l),2 



pq 



■\X) (5.52) 



where the presence of other boundaries and operator insertions is denoted by X. Further- 
more, it will be useful to have at hand some explicit formulae for amplitudes resulting from 
the application of (15.501) . For the (3,4) minimal string the full set of cylinder amplitudes 
are!, 

z\ - z 2 s , 1 



Z(l, 1; axil, 1; a 2 ) = 2(2, 1; ^2, 1; a 2 ) = ln(- 



Xi - x 2 



2V3e 



(5.53) 



Z(2,l;t7i|l,l;t7 2 ) 
<ZT(1, 1;o-i|1,2;ct- 2 ) 
Z(l,2;<7i|l,2;<7 2 ) 



- ln(^i + z 2 ) + 



1 



4v/3e 



Z(2,l;a 1 |l,2; ( r 2 ) = - ln(-l + 2z\ + 2^z x z 2 + 2^ 2 2 ) + — ^. 

2v 3e 

11 (zi +2 2 )(Zl -x 2 ) 4^' 



5 These expressions are correct up to an unimportant common normalisation constant and additions 
of numerical constants. 
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where we have introduced the following notation, 

z = cosh ( — — ) , x = -^7= = cosh ( — ^— ) , x = t-j= = cosh ( — — ) . (5.54) 

The first check of the SS relation was performed in [22] in which it was applied to cylinder 
amplitudes in which all boundaries carried boundary states of the form \a; a, 1). It was 
found that the SS relations are not satisfied identically in this case, however it was argued 
that the deviation was non-universal and therefore should be discarded. The evidence 
for non-universality was two-fold, corresponding to the two properties of non-universal 
terms introduced in Section 15.4.21 Firstly it was noted that the deviation possessed the 
following properties; 

Property 1 (PI): The deviation could be made independent of the bulk cosmological 
constant, and hence analytic in /x, by setting the IR cutoff e in (I5.50p to the volume of 
the Liouville direction. 

As an example consider the second of the (3, 4) minimal string cylinder amplitudes above, 

AZ(2, 1; 1; a 2 ) = - logfo + x 2 ). (5.55) 

4e 

Choosing e equal to the volume of the Liouville direction, \ = \ log — where A is a constant 
and recalling that x = \xb) ' \JJi we find, 

AZ(2, 1; 1; <r 2 ) = -- log A - log^Bi + ^b 2 )- (5-56) 

Since this is independent of \i it was interpreted as non-universal and therefore set to zero. 
For a demonstration that all cylinder amplitudes with boundary conditions given by the 
states \a\ a, 1) have this property the reader is referred to the appendix of [22] 



Property 2 (P2): The inverse Laplace transform, with respect to all boundary cosmo- 
logical constants, of the deviation minus any regularisation dependent parts is zero almost 
everywhere i.e. it is supported only at points. 

This property was uncovered by considering f l5.50p . Using ( I5.47P in (I5.5(jp we find, 

r°° sinh(^) l-niP 

Cr\z\(k,i-MV,s-M)<x \ dPMLi)ML2) . , , 9 ^ M iW— ), (5.57) 

Jo smh(27rP6) Jpq 
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where ipp(L) = wP sinh(27r P /b) K 2 ip/b{L) ■ By noting that the integrand is symmetric, 
the integral may be computed by closing the integral along the entire real axis by a 
semicircular contour. Since the Seiberg-Shih deviation is merely a linear combination of 
amplitudes its double inverse Laplace transform with respect to each boundary cosmo- 
logical constant also exists and is given by 

r°° sinh(^) 2mP 

£- 1 [AZ](A;,/;L 1 |r, S ;L 2 )cx / dPM^M^) . r (C T* AF M , r , a (-— ) (5.58) 

where 

AF k) i,r )S {z) = F k)l)T)S (z) - U k -i(cosqz)Ui- 1 (cospz)F 1At r tS (z) (5.59) 

and we have assumed that the Seiberg-Shih transformation has been applied to the (k, I) 
boundary. An important property of AF^i >ra (z) is that it is an entire function and 
therefore the only contributions to the integral comes from the poles at P = |§ where 
n G Z. In the case that (k, I) = (a, 1) and (I, s) = (a', 1) the residues of these poles are 
zero and the above expression is zero. 

Although contrary to the original motivation of P2, we have allowed for the possibility 
of the amplitude being non-zero at points away from zero. We merely require the weaker 
condition of having point-like support. 

We now want to review the results of [36] in which all boundary conditions on the 
cylinder amplitudes were considered and the consequences of this on the SS relations. 
Before proceeding it is important to note that the sum in the Seiberg-Shih relation does 
not respect the symmetry of the Kac table and so for a given amplitude there are two 
possible Seiberg-Shih relations that might possess PI or P2. However, there are cases 
where both possibilities lead to deviations that do not possess either property PI or P2. 

For an example of a deviation that does not possess property PI consider, 



AZ(l,2;a 1 |l,l;a 2 ) = In 



-1 + 2xi 2 + 2^/2x!X 2 + 2x 2 2 



(5.60) 



2(xi + x 2 ) 

or, since Z(2, 2; a x \l, 1; a 2 ) = Z(l, 2; a x \l, 1; <x 2 ), 

AZ(2,2;a 1 |l,l;a 2 ) = In [(-1 + 2x x 2 - 2^2x x x 2 + 2x 2 2 )l , (5.61) 



where we have suppressed the e dependent term. It is clear that since the argument of 
both expressions is not a homogeneous polynomial in x the dependence on \i cannot be 
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factored out and placed in a separate term. Hence the IR cutoff cannot be chosen to 
cancel all dependency on the bulk cosmo logical constant. 

An example of an amplitude not possessing P2 need only have the property that the 
residues are non-zero (I5.58p . leading to a function with global support. As an example 
consider, 



AZ(l,2;a 1 |l,2;a 2 ) = In 



X\ — X 2 

2(xi - x 2 ) 



(5.62) 



or, since 2(2, 2; cxi|2, 2; cr 2 ) = 2(1, 2; <7i|l, 2; a 2 ) 



AZ(2,2;ai|2,2;a 2 ) = In 



[x X - X 2 )(x\ + x\ - \) 

2(x\ - x 2 ) 



(5.63) 



By applying ( I5.58P it is now easy to show PI and P2 are not satisfied for (I5.62p and (15. 63ft ; 
in both cases there exist poles when P = ^ and therefore when the integral contour is 
closed around them it will result in a sum of terms containing Bessel functions of the form 
Kin(L) which have global support in L. 

Having shown that the cylinder amplitudes possess deviations from the SS relations 
which do not fit the non-universal classification advanced in [22] we now wish to motivate 
a possible extension of the criteria by which non-universality is judged in the hope of 
rescuing the SS relations. The observation, is that, the terms that disappear in (I5.50p 
under an inverse Laplace transform are dual under the Liouville duality to the terms 
that produce the troublesome poles in the above test of P2 for arbitrary amplitudes. We 
therefore propose that non-universal deviations should possess, 



Property 3 (P3): The deviation may be written in terms of contributions for which 
either the inverse Laplace transform with respect to all boundary cosmological constants 
[Mb or all their duals, [l b , has point-like support. 



We devote the remainder of this section to a more careful motivation of P3. Consider 
the integral representation of the amplitude (15.501) ; it may be computed by extending the 
region of integration to the entire real line and then splitting the integrand up into terms 
for which the contour may be closed in either the upper or lower half plane. Assuming 
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that o\ > cr 2 , then upon substituting in (IB.16P and (1B.17P this results in the expression, 
Z(k,l;a 1 \r,s;a 2 ) oc F k i rs (0)(— - nax)-^— + (5.64) 



E^PQtt i nn \TT i nn \TT i 7Tn \rr I TTTh * - "°l . TITIU 2 
tVfc_i(cos — )tv r _i(cos — )Ui-i(cos — )(7 s _i(cos — )e cosh — ■ h 
! w p p g q y/pq 

n^Omodp 
n^Omodq 



El. . „ „ . 7r?2 . sin(7m/o) , , , 

-(-l)"F fc< (_) ■ 1 ^ cosh(7r6na 2 ) + 

n g sm(7rpn/gj 



71=1 

n^O mod g 



E2. . _ „ .7TO. sin(7rn/n) Egggi , . 
- -l "^,i, r , s - ■ , ^ cosh 7rna 2 /6 + 

n p smynqn/p) 



n=l 
n^O mod p 



£ J.^!,*^ WTO)e -=5P cosh ^ 



;1 rapg v //*/ 

The above expression may be understood as representing a cylinder amplitude as a sum 
over disc amplitudes with various local or boundary operators inserted. If we compute 
the deviation using this expression the first sum cancels while the others remain with 
the replacement Fkj >r ,s ^-Fk,i,r,s- The sum in which each term is proportional to 
cosh(7r6rao"2) arises from the poles in (I5.50p due to the factor of sinh(^^) in the denom- 
inator; it is precisely these terms which vanish when we take the double inverse Laplace 
transform with respect to the boundary cosmological constant. These terms have the 
interpretation of some form of descendant dual boundary length operator inserted on the 
boundary of a disc. On the other hand, the sum in which each term is proportional to 
cosh(7mo"2/6) arises from the poles in (I5.50P due to the factor of sinh(27r6P) in the denom- 
inator; it is these terms which, under a double inverse Laplace transform with respect to 
the boundary cosmological constant, become the terms proportional to Bessel functions 
of the form Ktnilsi). Similarly, if we were instead to take the double inverse Laplace 
transform with respect to the dual boundary cosmological constant, the terms propor- 
tional to cosh(7mo"2/&) would vanish and the terms proportional to cosh(7m&er 2 ) would 
become terms proportional to KanfLz)- We therefore see that the deviation computed 

4 

from (15.641) can be understood as a sum of terms which vanish under a double inverse 
Laplace transform with respect to the boundary cosmological constant or its dual. 
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CHAPTER 6 



FZZT Branes in the (3,4) Minimal String at Higher 
Genus 



Testing the Seiberg-Shih relations using a Matrix Model 

In the last chapter we argued that the deviations for all cylinder amplitudes have property 
P3. Obviously it would be interesting to know whether this is true for deviations of any 
amplitude. Unfortunately amplitudes more complicated than the disc and cylinder have 
never been computed using the Liouville approach. An alternative is to use the matrix 
model formulation of minimal string theory in which the computation of amplitudes with 
arbitrary numbers of boundaries and handles is straightforward. 

The disadvantage of using the matrix model is that it gives the amplitude in a fully 
integrated and summed form and so it is hard to understand the structure of the amplitude 
in terms of continuum concepts such as states circulating in a loop. Furthermore, generally 
it is obscure how the graph labels map to the conformal field theory degrees of freedom 
and even more importantly, the matrix model appears to not contain all the boundary 
states present in the minimal string. Indeed, this was used as evidence in [17] in support of 
the Seiberg-Shih relation. However, there are special cases where the relation is manifest; 
we saw that the (3,4) minimal string admits a formulation as a matrix model (15 Ah in 
which the matrices are directly related to the spin degrees of freedom. In this model it 
is straight forward to construct all conformal boundary states. This is the model we will 
study in this chapter. In recent work of J37J [38] a number of other boundary states of 
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the (p, q) minimal string have been constructed in the matrix model formulation and it 
would be interesting to extend our results to those cases. 

In this chapter we will find it useful to change the definition of the coupling constants 
from those considered in (15.41) . Instead we will use the matrix model defined by, 

Z = J dM l dM 2 exp (~ Tr f-cM 1 M 2 + i(M x 2 + M 2 ) + ^(Mj 5 + M 2 3 )^ ^ (6.1) 

The coupling between neighbouring spins is controlled by the parameter c whereas g 
controls the cost associated with adding more vertices to the graph. It is g and c we use 
to tune the matrix model to its critical point which for our matrix model we do by setting, 

9 = 9c(l ~ a 2 r]fi), c = c c , (6.2) 

where rj is some constant which can be determined by comparing with the Liouville theory 
and the critical point is achieved by letting a — >■ Co We will often refer to this limit as 
the scaling or continuum limit. 

The conformally invariant boundary conditions of minimal CFTs are in one-to-one 
correspondence with the primary fields of the theory. For a CFT that describes the critical 
point of some discrete model the boundary conditions can sometimes be understood as 
universality classes of boundary conditions in the discrete model. As we approach the 
critical point of a discrete model all boundary conditions in a given class will flow to the 
same boundary condition in the CFT. In particular for the (3,4) minimal string its three 
conformal boundary conditions are the continuum limit of the discrete configurations of 
boundary spins consisting of all spin +, all spin — , or free spins. It is these boundary 
conditions we want to implement in the matrix model. 

The usual way of inserting a boundary in the matrix model is to compute the resolvent 
Wm 1 as it is a generating function for triangulations of discs. We will adopt the following 
notation for more general quantities, 

W Q V 1 Q2 f 2 .., Q3 Gi Q4 g 2 ...,..(xi, x 2 , ...:/„./,,,...: 
/!_ / Q 1 (M 1 ,M 2 ) Q 2 (M U M 2 ) 
\N T \x 1 -F 1 (M u M 2 )x 2 -F 2 (M ll M 2 ) 
1 (, Q 3 (M U M 2 ) g 4 (Mx,M 2 ) 
N T \ [ y 1 -G 1 (M 1 ,M 2 )y 2 -G 2 (M u M 2 ) 

1 The constants g c and c c have the well known values g c — 10c^ and c c = ^(2V7 — 1). 
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By tuning the Xi,jji, ... as we take the matrix model to its critical point we can extract 
continuum quantities corresponding to amplitudes with macroscopic boundaries. By way 
of example, for Wm 1 (x) we use (I6.2p and set, 

x = x c (l — a d K[iB) (6.4) 

where d is chosen to produce a non-trivial limit and k is chosen to agree with Liouville 
theory. As we let a — > 0, Wm 1 will have an expansion in powers of a □, 

W M Axc{l-a d K^ B )) = ^a^W i {^B^) + a dw WM 1 {^B^) + h.o.t . (6.5) 

i 

In this expression W^b, /•*) is analytic in both arguments whereas PVmi is defined as being 
the first term non-analytic in fi B and fi. It is Wm 1 that corresponds to the continuum 
quantity and so in this case is the partition function of the continuum theory defined on 
a surface with a finite sized boundary with the boundary condition determined by which 
universality class the discrete quantities (TrM") belong to. This is slightly complicated 
by the fact that the resolvent actually corresponds to a boundary with a marked point and 
therefore must be integrated with respect to the boundary cosmological constant before 
comparing to Liouville theory. This motivates introducing the integrated quantity, 



(6.6) 



where we have allowed for more general amplitudes which have more than one boundary. 
The generalisation to amplitudes of the form (16.31) is obvious. 

For the matrix model (16.11) it is easy to construct resolvents which generate boundary 
conditions that flow to the three different boundary conditions in the CFT: 

W Ml (x) = (±=Tr- 1 



N x - M x 
, 1 

W Mi +m 2 \x) = { TrTr- 



t N x-(M 1 + M 2 ) 

The first two resolvents Wm* generate surfaces on which only Mi vertices appear on the 
boundary. Since the M, vertices map directly to the spin degrees of freedom these flow to 



2 When taking the scaling limit in the remainder of the chapter we will choose r/ = 1 and k = y g c /c c 
It should be kept in mind though that by changing rj and n we can renormalise [Ib and /x. 
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the fixed spin boundary conditions in the continuum limit. The final resolvent generates 
boundaries in which both types of vertices appear with equal weighting. This will therefore 
flow to the free spin boundary condition. We can parameterise these resolvents in the 
following way, 

W*) = (j!-*jh) ^^M.-<^M 2 -< 1 + 4 (,8) 

This generates graphs with weighting corresponding to the Ising model on a random 
lattice but with a boundary magnetic field, to which the parameter a is related. Clearly, 
we can choose a = ±1 to reproduce the resolvents Wm v By taking a to infinity we may 
also obtain the resolvent Wm 1 +m 2 - One method to compute the resolvent Wx is to make 
the following change of variable in (16. ip . 



S = M 1 + M 2 + l + c, I = L^M 1 -^M 2 -(1 + /. (6.9) 



The partition function in these new variables takes the form, 



Z = J [dSdX] exp 



-— Tr (X 2 S + aXS 2 + V(S)) 
9 



(6.10) 



where, 



V(S) = 1(1 + 3a 2 )S* - C -S 2 + ~(3c - 1)(1 + c)S. (6.11) 



Generic values of a (i.e. the Ising model on a random lattice in the presence of a boundary 
magnetic field) were investigated using matrix model techniques in [3T], in which the disc 
amplitudes were calculated. The method employed in [31] to solve the matrix model is 
combinatorial and so is not easily generalised to compute more complicated amplitudes 
and so we will not employ it here. However, it is clear that if we are only interested 
in amplitudes in which all the boundary conditions are of the free spin form we may set 
a = in 1 16. 10p to obtain the 0(1) model whose solution is well known. Similarly if we are 
only interested in amplitudes in which all the boundary conditions are of the fixed spin 
form we may simply use (16. ip . This is in fact the case for most of the results we require for 
the remainder of this section and so they can be obtained form the literature [M] [291 130] . 
For more complicated amplitudes involving both fixed and free spin boundary conditions 
it is necessary to solve the matrix model (I6.10p ; this may be done using loop equations 
and we give details of this in the next section. 
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.1 1 Loop equations for SX Matrix Model 

The resolvents for the matrix model ( 16.101) can be obtained by generating an appropriate 
set of loop equations. We begin by obtaining the resolvents Wx and Ws- Consider the 
following change of variables in the matrix model, 

X^X + e—^— (6.12) 

z o 

this gives the loop equation, 

W x s (z) = | - ™Wstz) (6.13) 

Now consider, 

X ^X + e(^^ + ^^) (6.14) 
S ^S + e(-^^ + ^^) (6.15) 
X ^X + e^^^ + h.c) (6.16) 

Transformation (I6.14p together with (16.151) and (I6.13P gives, 

tt x (z, x) = Q s (z, x)VL S x(z, x) + — {Wsx;s(z, x; z) - aW S x-,x(z, x; x)) (6.17) 

g N 2 

where we have introduced the following functions, 

Q s (z,x) = W s (z) - aW x (x) - -(x 2 - a 2 z 2 + Viz)) 

9 

I fy 7 

n sx = W sx (z,x) + -(x-—) (6.18) 
9 ^ 

Q x = W£(z, x) + f + l(x- f ){x 2 - a 2 z 2 + V'(z)) 

where P = Q - a 2 S + ax - ^ and Q = (V'(z) - V'{S))/{z - S). 
Finally, the third transformation (16. 1 5[) gives, 

1 / 1 a 2 z 2 \ 

Q sx (z,x)Q sx (-z,x) = - W s (z) + Ws(-z) + aW x (x) + -(x 2 - — ) (6.19) 

9 \ 9 4 J 

JH W sx;Sx(z,x;-z,x) 
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Eliminating ilgx(z,x) between this and (I6.17P gives, 

1 f 1 a 2 z 2 

-n x (z,x)Q x (-z,x) = n s (z,x)n s (-z,x) W s (z) + W s (-z) + aW x (x) + -(x 2 - — 
9 V #4 

~Jp \9^s(z,x)fl s {-z,x)Wsx;Sx(z,X] -z,x) + Q x (z,x)W S x;s(-z,x; -z) (6. 
-afl x (z, x)Wsx;x(-z, x; x) + 0£ (-2, x) (w S x;s{z, x; z) - aW S x-,x(z, X] x)jj 

By substituting in the large N expansion for all resolvents appearing in this expression 
we can obtain a recursive definition for the genus h resolvents. The genus zero resolvents 
satisfy: 

-Q x ( °\z,x)n x m (-z,x) = (6.21) 
9 

Q s (z,x) {0) Qs(-z,x) {0) (Ws m (z) + W s {0) (-z) + aW x i0) (x) + - g (x 2 - ^f) 

This is an important equation and all subsequent loop equations will be related to this 
through derivative-like operations such as the loop insertion operator. It is important to 
note that the LHS is polynomial in z and so if we expand the RHS in x about any point, 
the resulting Laurent expansion must have coefficients that are equal to a polynomial 
function of z. A convenient point to choose for this expansion is infinity as the definition 
of the resolvent then coincides with the Laurent expansion, whose coefficients are the yet 
to be determined quantities, (TrX n ). Thus we produce a number of equations containing 
W s i0) (z), 

<lr jtls^xf^ {Isi-z^x^x (6.22) 



2nix 



(w s {Q \z) + Ws i0) (-z) + aWx i0) (x) + \x 2 - = p r , 



where is a contour around infinity, p n (z) is a polynomial in z with a finite number of 
unknown coefficients, n is an integer and the integral merely picks out the nth coefficient 
of the expansion. From the large x behaviour of Wx{x) ~ - it is easy to show that for 
n > 6 the LHS is zero and for n > 4 the LHS yields an expression with no dependence on 
W s {0) . The non-trivial cases occur for n = and n = 2 which generate enough equations 
to solve for Ws <y0 \z). The large z expansion of equations obtained for n < give relations 
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among the quantities (TrS n ). The explicit equation obtained for Ws^°\z) is, 
W s {0 \zf + 1 -(^^ + V'(~z)-2V'(z))Ws i0 \z) 2 (6.23) 
+ -^(y'{-z)(3z 2 a 2 -4V'{zj) -3z 2 a 2 V'{z)+W{z) 2 + P 2 {z)}w s {0) 



+ -L(P ( z ) + ( z *a 2 -V'(z))P 2 (z)) = 



where, 



4# 



P^z) = -Pi(z) + (Coefficient of (W s i0) (z))° inE22]for n = i) (6.24) 



The unknown constants in p n (z) are not all independent and they may be found in terms 
of the constants (TrS n ) by expanding the RHS of (16.221) about z = oo and equating p n 
to the polynomial part of the Laurent expansion. The requirement that the singular part 
of the expansion vanishes again gives relations between the quantities (TiS n ). 

To find Wx^°\ we again consider (16.211) . however we now know Ws^ in terms of a 
finite number of unknown constants. On the LHS we have the function Q x (z, x) where 
P is a polynomial in z and S as defined in (I6.18p . whose highest power for both z and S 
is d — 1 and so it contains d unknown functions of the form W x (x) for < n < d. 

We can generate a system of equations for these resolvents by again expanding both 
sides of the loop equation about infinity but this time in terms of z. This results in, 



dz 

~n+l 
C z 



-n x {0 \z,x)n x (0 \-z,x) - (6.25) 
.9 



^(z,x) (0) Q 5 (-z,x) (0) (Ws (0) (z) + Ws (0 \-z) + aWx m {x) + -(x 2 - 



a 2 z 2 . 







9" 4 

Again we may vary n to generate different equations, however we now generate non- 
trivial equations for all < n < d, which is enough to solve for all unknown resolvent 
functions appearing in fl x (z, x). However, unlike for which we obtained an 

equation for general V(x), we were unable to do this for the equations for W x "(x). We 
also note that for the equation giving Wx , the highest power of Wx^ appearing depends 
on the order of V(x). 

The approach thus outlined gives an algorithmic solution to this matrix model with 
an arbitrary potential and results in algebraic equations for all resolvents. 

We now apply the above procedure to the case when the potential is (16.111) . There 
are two unknown constants in Pq(z) and P 2 {z) corresponding to (TrS 1 ) and (TrS 2 ). Their 
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values are determined by the requirement that (16.231) is a genus zero curve - which, in 
this case, is equivalent to the one-cut assumption. 

For Wx we can generate two equations from (16.251) . for n = and n — 1, containing 
Wx^ and . The equation for Wx^ resulting from the elimination of W§-^ has a 

number of interesting properties. Generically it is quartic in Wx^\ however it reduces to 
a cubic for a = ±1 and it reduces to a quadratic of (Wi' ' + polynomial in x) 2 for a = 0. 
Because of this the scaling limit for generic a does not hold for these particular values of 
a and must be taken separately. 

As was discussed before, the resolvents corresponding to the spin + and — Cardy 
states are Wm-l and Wm 2 ^ respectively, which may be computed directly from (16. ip 
and are clearly identical. However a faster method, since we already have Wx^\ is to set 
a = -1 implying X = M x + \{l + c). 

It is easy to generate loop equations for a general 2-loop amplitudes of the form Ws-h 
and Wx-h where if is a string of S and X matrices of length A/#. We can write H as, 
n =i Xn, where \ n is a matrix defined by, 

, X, n e ly , 
^{S.nel* (6 ' 26) 

and Ix and T$ are disjoint indexing sets. For such a product, the notation Hu j\ is 
defined as , 

3 

H {i , j ) = Hxi (6.27) 

n=i 

The loop equations for Ws-h, can be generated by the following changes of variables. 

X -> X + 6—^—TrH(S,X) (6.28) 

z o 

s - s + f (^^ + ^7 1 s) M(s - Jf) (6 - 30) 

X -^X + £ (^^^ +ft .c)lV//( 5 ,X) (6.31) 



The loop equations generated by X — > X + eF(X, S)TrH, may by obtained from the loop 
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equations generated by X — > X + eF(X, S) by the following procedure, 



Wg 1 (X,S)Wg 2 (X,S) -> W Gl (X,S)W G2 (X,S)-H + Wg 2 (X,S)Wg 1 (X,S);H 

—( ) —( 



(...)+ E^cm)W)%,^)(^) ( 6 - 32 ) 



where Gi(X, S) is any function of X and S 1 . An equivalent set of rules applies for loops 
equations generated by S — > S + eF(X, S)TrH. The resulting loop equations can then be 
solved by using the same method as used to solve the one-loop amplitude equations. For 
example, the two loop amplitude Ws-,s(z] u) m &y be obtained by choosing H = (u— S)^ 1 
in (I6T28D . flBT29|) . flOO]) and flBT3T]) . Applying the rules in f[Q2]) we get, 



ifi£ ;5 0,x;w) = Q S X;s{ z i x '> u )^s{z,x) + 

^sx(z, x)Q S; s(z, x; u) + d u 



Wsx(u)-W S x(z) 
u — z 



(6.33) 



and 



n sx (z,x)n S x;s(-z,x;u) + n sx -s(z,x,u)n sx (-z,x) = 

- (W S ;s(z; u) + Ws-s(-z; u) + aW x , s {x; u)) (6.34) 
9 

By eliminating Qsx-,s(z, x, u) we obtain an equation relating Ws-s(z; u), Ws-,s(— z] u) and 
W X; s(x;u). If we arrange the resulting equation so that the LHS is polynomial in z, 
we may expand the RHS in x about x = oo to generate a number of equations for 
Ws ; s(z;u) and Ws-s(—z;u). This is exactly the same procedure we used to compute 
the disc amplitudes in which the resulting expression contained a number of unknown 
constants, which for the equal potential case were (TrS) and (TrS* 2 ). However here we 
get an expression for Ws-s{z\ u ) m terms of unknown functions of the form W -s{u) and 
IV s2 -s(u). These unknown functions may be found by computing the expansion of Ws-s 
at z = oo and, by the symmetry of Ws-,s> equating it to the expansion at u = oo. The 
expressions for W s - t s(z) and W s2 -s( z ) then depend on a number of unknown constants of 
the form, (TrS' m TrS' ra ) c , not all of which are independent, as can be seen by considering 
the large z expansion of the W s -s(z) and W s2 -s(z). In the end the independent unknown 
quantities correspond to (Tr5'Tr5') c and (TtSTtS 2 ) c which can be fixed by requiring that 
W -s(z) has no singularities besides those at the branch points of Ws- The loop equations 
then give Wx s m terms of Ws-s- 
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In order to compute the 1/N corrections it is clear from (I6.2(jp that Wsx;S> Wsx-,x 
and Wsx-,sx need to be calculated. The calculation of Wsx-,x an d Wsx-,sx, follows a 
similar procedure however the resulting loop equations contain new amplitudes of the 
form Wsx..sx- The calculation of these new quantities does not pose a great difficulty; 
appropriate changes of variables give loop equations very similar to the ones presented 
previously. We include the necessary change of variables in Appendix 

Once these quantities have been calculated and substituted in (I6.20p . we may again 
expand both sides around x = oo and then z = oo to find an expression for Ws(z) 
and Wx(x) . Again this procedure introduces a new unknown constant corresponding 
to the J? correction to (TrS 1 ) which may be determined by requiring that W s (z) {1) does 
not possess any poles besides those at the branch points of Ws(z)^°\ 

6.1.2 | Disc, Cylinders and Disc-with-Handles. 

To take the scaling limit of Wx at a = —1 we may use the results of [24J or the results in 
the last section. The critical value of x is x = — c and the solution to the loop equation 
gives, in the scaling limit, 

(4 4 \ 

(f*B + ^I*b-p\ + \ Pb - J ■ ( 6 - 35 ) 

For Ws we may use [291 EQ] or solve f)6.23p . The critical value of x is at x = and 
computing the scaling limit of Ws^ we get, 




W Ml+Ma (jM B ,fi) = -^- c ( l^+V^-H +l^-\l^ L -»\ I • (6-36) 

Here we run into one final technicality. For a resolvent that produces graphs with some 
configuration of spins on the boundary Wf(m 1 ,m 2 ), there exist an entire family of resol- 
vents which produce the same configuration of boundary spins but with different overall 
weighting 

Wi F{Ml! M 2) (x) = XW f{Mu m 2 )(\x), (6.37) 

where A is a constant. Changing A has the affect of renormalising fis, however we want the 
same normalisation for all continuum quantities and so this fixes A to be \/2 for Wm 1 +m 2 - 
Therefore using (I6.36P and (16. 35ft . we have that 

V(M 1+ M2)/y/2 (JJ>b,P>) = V2u Mi (lib,v)- (6.38) 



6.1 Testing the Seiberg-Shih relations using a Matrix Model 
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This is precisely the relation between these two disc amplitudes found in the continuum 
calculation. 

One may wonder what the scaling limit of Wx^ is for generic values of a; given that 
there are only a finite number of Cardy states the expectation is that there should be no 
new non-trivial scaling limits. For generic a the critical values of x occur at -^7 and 

a— 1 a+1 

Taking the scaling limit, for both critical values of x, gives 

U-X/ai^B, A 1 ) = Wm^Hb, fi) (6.39) 

which confirms our expectation. We have therefore demonstrated that we can build 
all the boundary states of the (3,4) minimal string. This will be confirmed further by 
reproducing the cylinder amplitudes shortly. The nature of the scaling limit when a = 
will be addressed in the Section 16.1.31 

The one case in which we are interested which cannot be found in the literature is 
Ws-x- The computation, using the formalism of the preceding section, results in 

cu s .,x(ijbi, Hb 2 , ft) = -^\og(-l + 2zl + 2V2z lZ2 + 2zl) (6.40) 

for the scaling form of Ws-x, which is in exact agreement with the Liouville result (up to 
renormalisations of /i# and /x). Of course the other cylinder amplitudes may computed 
in this matrix model or using [24] |30] and we find the scaling forms to be, 

Ux;x(hbi, fj) = In — — — , (6.41) 
10c~ X\ — x 2 

UX;y(^Bi, VB2, AO = ~Yb^2 ^ ^ + * 2 ' 

W5;s(/iBi, Hb 2 , AO = ttt^ In 7 r, (6.43) 

10c^ (zi + z 2 ){x 1 - x 2 ) 

again in agreement with the Liouville calculations. The fact that these amplitudes agree 
with the Liouville calculations including terms that were classified in |22j as non-universal 
suggests they really should have a physical meaning (and shouldn't be thrown away). This 
issue was raised in [26J. 

In order to test if P3 is a property of all deviations we need to examine more com- 
plicated amplitudes. The 1/N correction to the disc amplitude is the amplitude for a 
disc-with-handle. These computations have already been done in [30] and [23], the result 
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being, 



z(7 - 24z 2 + 48z 4 
235§36c c (-1 +4z 2 y 



(i) , , V2 z(5-2Az 2 + 48z 4 ) _ 2 

V l+ *,/>'") = 2hh6Cc (-1+4^)3 " " (6 ' 45) 



If we now compute the deviations we find 

J M 2 



Au>$ a (ji B ,fi) = 0, (6.46) 



(1) 2e4 ,2(-3 + 41^- 208z 4 + 416,2 t, -512z 8 + 256z-\) r 



6 



Au (m 1+ m 2 )/V2^b,») 5 | 9cc (-1 + 22) 3 (1 + 2z) 3 (l-16^ 2 + 16z 4 ) 3 M 

We now want to consider if Aw|^ +M ^^(/iB, /i) possesses P3. For it to possess P3 
we require that it can be expressed as, 



where C~ l [A] and C% l [B\ have point-like support. Since the definition of the deviation 
is a linear combination of a number of amplitudes, the functions A and B will also be 
a linear combination of contributions from each amplitude. Since the theory is unitary 
the contributions to A and B from each amplitude will be positive. From the expressions 
in (I6.44p we see that each amplitude goes as at most ~ - for z — > oo and therefore any 
contribution to A(x) or B(x) must go as at most ~ x^ 1 ^ 3 or ~ £~ 1//4 respectively. We 
therefore may conclude that A(x) ~ a; -1 / 3 and B{x) ~ a;~ 1//4 as x — > oo and hence that 
their inverse Laplace transforms exist. Furthermore, given that as z — > oo, 

Auj( Sl 1+ M 2 )/V2^ B > A ~ J ~ ^73 ~ ' ( 6 ' 48 ) 

we see that A and -B can not fall off faster than = ax~ 7 ^ 3 + o(x^ 7 ^ 3 ) or B(x) = 

6x _7 / 4 + o(x -7 / 4 ). However it is a well known theorem that a function with only point-like 
support may be expressed as a sum of derivatives of ^-functions; 

oo oo 
m=l n=l 

Such a function will have a Laplace transform of 

oo oo 
m=l n=l 



5 We denote in the subscript the variable the inverse Laplace transform is taken with respect to. 
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Given that we require that and are point-like supported this implies A(x) 

and B{x) have an asymptotic behaviour as x — > oo inconsistent with (I6.48j) and hence 
the deviation can not possess P3. 



6.1.3 | Dual Branes 

In [T7j the two matrix model description of the (p, q) minimal string was considered and 
it was shown that the continuum limit of the resolvent for one of the matrices was the 
disc amplitude with a fixed spin boundary. The continuum limit of the resolvent for the 
other matrix gave the dual brane amplitude. In our matrix model (16.101) it is not obvious 
if we can construct a correlation function which gives the disc with dual brane boundary 
conditions. In fact it is present in the form, 

W *V = ( X- ( m|-m,) ) < 6 - 51 ' 

which may be obtained from the general X resolvent (16.81) by setting a = 0. The resulting 
loop equation is easily solved as it is quadratic in x 2 , however the scaling limit is obtained 
using a different scaling of x, 

x = x c + a^pLB\fgJc 1 (6.52) 

which results in 

W k{Ml -M2){^B^) = _ ^|~|~ (j^ B + ^P-B-fi) + ~ \fi*B-frJ J ' ( 6 - 53 ) 

where k is a numerical constant. If we identify Jxb in this expression with the dual 
boundary cosmological constant then this is precisely the amplitude expected for the dual 
brane obtained from the spin + or spin — boundary state. This explains the a = 
possibility we left unexplained earlier in Section 16.11 



6.2 | The Seiberg-Shih Relation and Local Operators 

Suppose that the boundaries considered in the previous sections could be expressed as a 
sum over local operators in the theory, such as, 

oo 

(k,l;cri\ ~J2fn(k,l,vi)(O n \, (6.54) 

n=l 
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where O n represents some generic local operator and ~ means equivalence up to terms that 
diverge as o\ — > oo. If this were true then the effects of the Seiberg-Shih transformations 
on a boundary would reduce to studying their effect on the above series. The expansion 
of a boundary in terms of local operators was investigated in |25j in which it was used to 
compute the one and two point correlation numbers of tachyon operators on the sphere 
for the (2, 2p + 1) minimal strings. Later it was shown [27] that such a technique could 
also compute the three point function correctly for the (p,p + 1) models. 

One technical point raised in [25] and [27] was that the boundary state could only 
be expressed as a sum over local operators if we included operators not in the BRST 
cohomology of [28]. That the operators in the cohomology of [28], which we refer to as 
LZ-operators, are not all the physical operators of the minimal string is easily verified by 
noting that the cohomology does not contain the boundary length operator § dM dxe h< ^. 
Such an operator can be inserted on the disc by differentiating with respect to the bound- 
ary cosmo logical constant, giving the result 



where A(p, q) is a constant of proportionality that only depends on p and q. Indeed there 
exists a whole family of such operators with Liouville charge b p+q 2 ~ nq where n E Z + which 
we will denote by U nq . 

In order to compute the coefficients f n (k,l,ai) appearing in (16.541) we need to know the 
amplitude for any local operator inserted in a disc with a (k, I) FZZT boundary condition. 
For the tachyon operators the amplitude is (15.451) . For the non-tachyon operators we have 
argued using the mini-superspace approximation and the exact conformal bootstrap on 
the disc that the functional form of the one-point disc amplitude is (15.491) . Finally, we 
have seen in (I6.55P that for the first of the non-LZ operators the disc function has the 
functional form expected from a direct application of the mini-superspace approximation 
and also the exact conformal bootstrap. We therefore have, 




(6.55) 




(6.56) 



If we apply the Liouville duality transformation to the above expression we get 




(6.57) 
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where we have introduced the operators U np which have Liouville charge b p+Q 2p np . The 
expansion (I6.54p can now be given a more concrete form, 

oo 

(k, I; (7i| ~£)/ n (M,*i)(Z4j- (6-58) 

n=l 

Having introduced the necessary technical results we now want to prove the following 
claim, 

If all the FZZT branes in the (p, q) minimal string can be replaced by local operators 
then the deviation from the Seiberg-Shih relations is caused only by the non-LZ opera- 
tors. 



Consider the cylinder amplitude 2(k, I; <Ji\r, s; 02). If we replace one of the boundaries by 
a sum of local operators using (I6.58P then we get, 

Z(k,l;ai\r,s;a 2 ) = f n (k,l,a 1 )(U n \r,s;a 2 ). (6.59) 

n=l 

By comparison of (I5.64p with (I6.59p . we find for the coefficient of the LZ operators, 

An 1 TTTl Tin n _ 7rn<J l ~ _ 

f n (k,l;ai) = -C4_i(cos — )^-i(cos — )u,~^e = f n (k,l)e vw, (6.60) 

A D [n) n p q 

( \ 1/2 

where Ac is just a numerical constant, Ar){n) = T(l — ~)r(— ^) f sin(^) sin( z ^ : ) j and 
n is not a multiple of p or q. Furthermore, note that the LZ operator coefficients satisfy 

fc-i 1-1 

(6.61) 

o=-(fc-l),2 6=— (i— 1),2 

This shows that if all boundary states admit an expansion in terms of local operators then 
the Seiberg-Shih relations transform the coefficients of the LZ-operators correctly. Hence 
any deviation from the Seiberg-Shih relations must come from the non-LZ operators. 
However, it is the non-LZ operators one would expect to give deviations compatible with 
P3 as they correspond to boundary operators and their duals. One might then suspect 
that in fact the boundary states may not be expressed as a sum over local operators. We 
will now see this suspicion is borne out by examination of the cylinder amplitudes. 



104 



Chapter 6. FZZT Branes in the (3, 4) Minimal String at Higher Genus 



If we note that for cylinder amplitudes, Z{k, Z; &\\r, s; 02) — Z(k, Z; cr 2 |r, s; ci) then 

together with (I6.58p . this implies, 

00 00 

) j f n (k, l, ax) (U n \ r, s; a 2 ) = }J n (r,s,ai)(U n \ k,l;a 2 ). (6.62) 

n=l n=l 

Since (U n \ r, s; a) = An'^ g n (a) where g n form a set of linearly independent functions then 
this implies 

f n (k,l,a)=A^h n (a), (6.63) 

where h n is some function. We conclude that if all states admit an expansion of the form 
(16.581) then cylinder amplitudes can be expressed as, 

00 

<Z(M;<Tl|r,a;<7 2 ) = J2 A n' l)A[ n S)h n(*l)9n(cT2). (6.64) 
n=l 

We now want to see if the cx-independent coefficients of each term in (15.641) has the above 
property. Consider the coefficients appearing in the second sum of (I5.64p . 

Tin sm(Tin/q) 

Fk,i,r,s{ — ) -t— ; j-t = (6.65) 

q sm{Ttpn/q) 

X p (k,r) A g (/,s) 



E E ^ru,J^Jf)U m -^)-A 

u ^1 i i o A— 17 „liio \ H / \ F / 



where X p (k, r) = min(A; + r — l,2p — 1 — k — r) and n ^ mod q. We want to know 
whether this can be written in the form An' 1 ^ An' s ^ ■ That this is not true in general is 
shown by considering the case (p, q) = (4, 5), for which we have, 

?m sin(7m/g) 

-^2,2,2,1 (. — r^r - 0, (6.66) 

g sm{TTpn/q) 

ttti sin(7rn/g) 1 

^1,1,2,1 — —7 = -, 6.67 

g sm(7T]9n/g) 2 

which is clearly incompatible with the factorisation required (I6.64p . This lack of factori- 
sation appears to be generic; for instance if we consider the (3, 4) model we find that there 
is not even a subspace of states which admit an expansion in terms of local operators. 
This can be seen by considering the matrix, 



( nn^ sin(7m/g) 

-T k,Lr,s\ 



q sm(7rpn/q) 



( 2 -t^r 1 2 C0S 3^ \_ 

1 2(-l)" 2/ i\n+l m „3fn 

3 3 3 \ l ) LOb 4 

\ fcos^f §(-l) n+1 cos^f §(-l) n (3 -4cos 2 ^f ) J 

(6.68) 
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where the columns and rows correspond to (k, I) = (1, 1), (2, 1), (2, 2). If there exist A n ' 
and An'^ such that F^i^si 2 ^) — An An for some subset of the states, then at least 
one of the cofactors of the matrix in (16.681) must be zero, which is not the case. 

The same conclusions can be reached by studying the other coefficients associated 
to the non-LZ operators in (15.641) . This leaves us in an awkward position; there is no 
subspace of states that can be expanded in terms of local operators unless we allow the 
subspace to be one- dimensional (in which case why choose a particular boundary as being 
the one that can be expanded in terms of local operators?), yet from the results of [25] 
and [27] it is clear information about insertions of operators can be extracted from loops. 
It is worth noting that from the examination of cylinder amplitudes in which one of the 
boundaries is a ZZ brane, it appears that ZZ branes do not couple to non-LZ states 
[36] [39] HQ] and so some of the above problem might be avoided if a ZZ brane is present. 

.1 1 Extracting Local Operators From Loops 

An interesting observation is that all correlation numbers were extracted from loops with 
fixed spin boundary conditions. We now will present further evidence that correlation 
numbers can be extracted from the fixed spin boundary by reproducing the results of 
[26] for the 1-point correlation numbers on the torus in the (2, 5) minimal string. This 
also serves as an independent check on the computation in [26] as our methods are quite 
different. 

A trick that makes this computation easy is that all amplitudes for the minimal string 
are easily expressible in terms of z, which is the uniformising parameter of the auxiliary 
Reimann surface [T7J or equivalently spectral curve. If we write these amplitudes in 
terms of a new variable, w, defined by z = \{w + w~ l ) we may easily compute the large 
w expansion of the amplitudes. Since z = cosh^p=, the large w expansion will be an 



expansion in terms of the functions e which is exactly what is required to compute the 
correlation numbers. Explicitly, if we have an amplitude (1, 1; ai\ X) then 



oo 




(6.69) 



n=l 



n=l 



but we also have 



oo 




(6.70) 
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giving, 

(U n \X) = „ An . (6.71) 

/n(l,l) 

We now have the tools to extract the correlation numbers on the torus beginning with 
the disc-with-handle amplitude. The calculation in [26] of the torus one-point correlation 
numbers was performed for the (2, 2p+l) minimal string. The disc-with-handle amplitude 
for the (2, 2p + 1) minimal string has been computed numerous times using matrix model 
techniques. We will restrict ourselves to the case of the (2, 5) model in which the disc- 
with-handle amplitude takes the form, 

, l + 12z 2 _7 . , 

wm(/^b,a0 = 5 A* 4 - ( 6 - 72 ) 

z 6 

The boundary condition on the disc is the equivalent of our spin up boundary as it is 
computed using the resolvent of M where M is the only matrix appearing in the matrix 
model. We will now use this to compute the torus one point correlation numbers, for 
tachyon operators, by replacing the boundary with a sum of local operators as in f !6.58p . 
That this computation produces results that match exactly the results of [26] is evidence 
that we may replace this boundary by a sum of local operators. Obviously, in order 
to compare results we have to renormalise the bulk cosmological constant, boundary 
cosmological constant and M to be consistent with [25]. In order to avoid doing this we 
note that the effect of such a renormalisation will to be to change the amplitude by a 
factor dependent only on p and q and so will not affect the ratio of correlation numbers. 
We therefore shall compare the ratio of correlation numbers. Computing the large w 
expansion of ( 16. 72ft we get, 

oo 

oo M (fiB, = ]T(-ir +1 (n - 2)(n + S)^ 2 "" 1 (6.73) 

n=0 

where the series converges if w > 1. We therefore find the torus 1-point correlation 
number for tachyon operators in the (2, 5) minimal string to be, 

W (1) = -^(sm — sin— J r(l--Jr(--J(-l) — ( n -5)(n + 5)//— (6.74) 

when n is odd and zero if n is even. The 1-point correlation numbers on the torus for the 
(2,2p+ 1) models were computed in [26], giving, 



lu) m_, u) m, n ^ w V rt 1 - w> r (i) (2p + 1 - »)(2p + 1 + ?) 
-< Ml> <- 1)J la^) r(i - ^jrfi) MiTT, " 



(6.75) 
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for n odd and zero when n is even. Using standard T-function identities it is then easy 
to show that the above equation is consistent with our computation (16.741) . 

The reproduction of the torus correlation numbers is evidence that this procedure is 
correct and that the fixed spin amplitude can be used to compute correlation numbers of 
local operators. Furthermore, we can obtain evidence that the non-fixed spin boundary 
states do not admit an expansion in terms of local operators by applying the above 
procedure to the disc-with-handle amplitudes we computed for the (3, 4) model. We have 
the expansions, 

"1 1 5 17 



12 7 
2 3 5 3 cA 6 



UJ 



(i) 

(Mi+M 2 )/V2 



V2 



.1.2 .7 
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+ 



+ 



72w 
7 



24w 72w 5 72w 7 72w 



ii 



(6.76) 
(6.77) 



If the expansion in terms of local operators is valid for both these boundary states then 
we would expect the coefficients of each term in the large w expansion to be related by a 
factor of \/2 as can be seen from (16.601) ; this is clearly not the case. The implication of 
these results is that the fixed spin boundary condition is special and is the only one that 
admits an expansion in terms of local operators. 



Discussion 

In this chapter we have examined the Seiberg-Shih relations ( 15.521) for cylinders and 
discs-with-handles. We saw in the previous chapter that for cylinders the identification 
of FZZT branes is naively spoiled by terms arising from the gravitational sector of the 
theory, as first seen by [22] [36] . We found that the terms spoiling the identification 
cannot always be associated with degenerate geometries of the worldsheet, in the sense 
of [22], as other geometries which are degenerate in a dual sense are also present. This 
lead us to conjecture that the terms spoiling the FZZT brane identification could always 
be interpreted as arising from degenerate and dual degenerate geometries. We checked 
this conjecture by computing the disc-with-handle amplitudes for free and fixed spin 
boundary states in the (3, 4) minimal string using matrix model techniques. We found 
that the deviation from the Seiberg-Shih relations in this case could not be written in 
the way we conjectured. Given that our conjecture was a very generous interpretation of 
what terms might be unphysical we conclude this is strong evidence that the FZZT brane 
identification conjectured in [17] does not hold at all levels of perturbation theory. 
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We also considered an alternative approach to testing the Seiberg-Shih relations by 
expanding boundaries in terms of local operators. If such an expansion were possible 
we showed that it would lead to deviations from the Seiberg Shih relations consistent 
with our conjecture. We then gave explicit examples for cylinder amplitudes where a 
local operator expansion of boundary states fails. This lead us to a paradox; how could 
correlation numbers of local operators be extracted from boundary states in [25] [27] 
if such boundary states cannot be expanded as local operators? To investigate this we 
computed the disc-with-handle amplitude in the (2, 5) model using matrix model methods 
and then showed that the fixed spin boundary condition yielded an expansion from which 
correlation numbers in agreement with the results of [26] could be extracted. This lead us 
to conjecture that the fixed spin boundaries are special in that they allow an expansion 
in terms of local operators. 

There are many obvious generalisations of the work done here, the most obvious being 
that it would be interesting to extend the above results to the (p, q) case. Since the 
Liouville techniques are not suited to higher-genus computations such a project would 
have to be tackled using matrix model methods. The main problem with this approach is 
that it is difficult to represent all the states for a given minimal string in a matrix model; 
in the above we chose the Ising model precisely because there was an obvious mapping 
between the spin degrees of freedom and the matrix fields. However this problem may 
now be less serious given recent work [371 EE] in which boundary states were constructed 
in the matrix model formulation for (2, 2p + 1) minimal strings. Another strategy for 
attacking the same problem may be the geometric recursion techniques developed by 
Eynard et al. in [65] as this provides an efficient way of computing many amplitudes 
for a given matrix model. Such methods would also allow a relatively easy computation 
of the cylinder-with-handle amplitudes; these would be interesting to study as one could 
then check if quantum corrections destroy the fact, noted in [361 EHl HQ], that FZZT-ZZ 
cylinder amplitudes are consistent with the Seiberg-Shih relations. 

Another generalisation would be to perform a similar investigation in other non-critical 
string models. An obvious choice would be non-critical superstrings as they perhaps 
have more physical relevance in addition to being much better behaved. A second and 
perhaps less obvious choice of string model would be the causal string theories developed 
in [321 E21 EI]- Currently these models have only been solved for a target space of zero 
dimensions and so such a project would require a generalisation of the models to include 
worldsheet matter. However, even in the zero dimensional case such causal string models 
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are better behaved and many of the odd features of the usual non-critical string are 
absent due to the lack of baby-universe production; in particular the gravity sector in these 
theories seems to be much weaker. Since the terms spoiling the FZZT brane identification 
are due to the gravity sector (at least for the cylinder amplitudes) might they have an 
interpretation in terms of baby-universe over production and if so might they be absent 
in these causal string models? 

Finally, we saw that the matrix model is able to reproduce all the boundary states 
found in the (3, 4) minimal string apart from the dual of the free spin boundary condition. 
In [T7] this was enough as they claimed that there is only one brane and hence only one 
dual brane in the theory. However, in light of our results there should exist other dual 
branes in the theory. Are such states present in the matrix model and if so can they be 
represented as some form of resolvent? This question and the others outlined we leave to 
future work. 
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CHAPTER 7 



Summary 



In this thesis we have reviewed how random graphs provide a very useful tool by which 
observables in a variety of two dimensional quantum gravity theories may be computed. 
Using these techniques we were able to investigate conjectures related to two of these 
observables; the dimension of spacetime and Hartle-Hawking wavefunctions, or more gen- 
erally spacetime containing a number of boundaries. 

In Chapter [2] after introducing the approach to two dimensional quantum gravity 
in which the metric is treated as the fundamental degree of freedom we reviewed how 
the resulting path integral could be defined via a lattice regularisation referred to as 
Dynamical Triangulation. This lattice regularisation can naturally be interpreted as a 
random graph and we reviewed how the disc-function observable could be computed 
in this approach via matrix model and combinatorial methods. We then explored the 
consequence of the ambiguity inherent in the metric-is-fundamental approach resulting 
from the necessity of choosing a class of metric to integrate over. In DT, the class of 
metrics integrated over were the Euclidean metrics on the manifold. We reviewed how 
the approach of CDT could be use to define path integrals in which only geometries with 
a well-defined casual structure where included in the integration. These theories differed 
from DT due to the propensity of DT to generate baby-universes throughout the space; 
such behaviour ultimately lead to a fractal geometry for the spacetime. 

The question of how to investigate the fractal nature of the spacetimes was taken up in 
Chapter [3] in which the concepts of Hausdorff and spectral dimensions were introduced as 
ways of generalising the notion of dimension beyond smooth manifolds. These measures 
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of the dimension were useful in characterising the properties of the space; we saw that 
DT gave rise to fractal spacetime since dh = 4 whereas d s = 2. For CDT on the other 
hand we saw that = d s = 2 signalling the emergence of a smooth spacetime in the 
continuum limit. We then reviewed the recent evidence obtained in a variety of approaches 
to quantum gravity, most notably CDT, for the spectral dimension being lower in the UV 
than IR. We referred to this phenomenon as dimensional reduction. 

In Chapter H] we constructed a toy model based on random comb graphs, in which 
dimensional reduction occurs. This was useful for a number of reasons; it shows that a 
scale dependent spectral dimension may be given a rigorous definition and that models 
can be found in which this rigorous definition results in a phenomenon similar to the 
dimensional reduction observed in CDT and other approaches. We gave a characterisation 
of possible continuum limits and the behaviour of the spectral dimension in each case via 
Theorem (J2J). As discussed in Section I4TB1 we hope that these results may be used as the 
basis for constructing more realistic models in which a similar dimensional reduction can 
be observed. In particular we hope the work in [16] , in which a type of random multigraph 
that shares many features with CDT were introduced, may be extended in this direction. 

In Chapter [5] we gave a brief review of minimal models, Liouville theory and string 
theory. This was motivated by the desire to couple matter to the two dimensional gravity 
models considered in the earlier chapters. After showing how the matrix model techniques 
are easily extended to describe DT models coupled to matter, we reviewed the equivalence 
between these models and certain string theories known as (p, q) minimal string theories. 
We then considered the boundary states, or branes, in these (p, q) minimal strings, in 
particular focusing on a conjecture of Seiberg and Shih which claimed that all branes 
occurring in such theories were physically equivalent. This claim is challenged by noting 
that the simple form of this conjecture fails for cylinder amplitudes. We reviewed the 
suggestion that since the terms spoiling the conjecture have particular properties, they 
are non-universal and therefore these examples should not be seen as counter examples. 
However, for more complicated cylinder amplitudes the spoiling terms fail to have the 
properties necessary to be classed as non-universal. We introduced a conjecture, referred 
to as P3, which gives a natural way, based on the duality of Liouville theory, in which 
these new spoiling terms may also be considered non- universal. 

In the final Chapter |6]we considered the conjecture P3 by checking if it held for the case 
of disc with handles in the (3, 4) model. This required a matrix model representation of all 
the conformal boundary states in this model. To accomplish this we introduced a matrix 
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model describing all the Ising model on a random lattice in the presence of a boundary 
magnetic field. We demonstrated this matrix model could be solved for any potential, 
not just the cubic one necessary for the Ising model. Using this solution we computed the 
disc-with-handle amplitude and argued that the terms spoiling the Seiberg-Shih conjecture 
did not satisfy the extended conjecture. This strongly suggests that the naively different 
branes in this theory are indeed different. There exist a variety of extensions of this work 
as discussed in Section 16.31 

Finally we considered an alternative approach to testing the Seiberg-Shih conjecture 
by writing the boundary states as a sum of local operators. This is a particular elegant 
way of approaching the problem as it negates the need to consider individual amplitudes 
separately. We found that the local operator representation of the boundary states lead 
to non-trivial terms which spoil the Seiberg-Shih relation, however, these terms appear 
to be of exactly the form necessary in order for conjecture P3 to hold. We investigate 
this apparent contradiction and argue that the expansion in terms of local operators only 
holds for the identity brane. In doing so we verify certain recent results of Belavin et al. 
[26]. 

We hope that this thesis has convinced the reader of the usefulness of random graphs 
in the study of low dimensional theories of quantum gravity and furthermore that there 
is a great deal left to explore. 
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APPENDIX A 



Mellin Transforms and Asymptotics 



Standard generating functions 

We record here a number of standard results for generating functions for random walks 
on combs; the details of their calculation are given in [3]. 
On the empty comb C = oo we have 

P^x) = l-v£, (A.l) 
G^(x;n) = {l-x) n ' 2 - 7 =^ } =-. (A.3) 

°° V ' ^ V ' (1 + y/x) n - (1 - y/x) n V ' 

Note that we can promote n to being a continuous positive semi-definite real variable 
in these expressions; G^(x;n) is then a strictly decreasing function of n and P< n (x) a 
strictly increasing function of n. The finite line segment of length i has 

which it is sometimes convenient to write as 

Pe(x) = v^rtanh (moo(x)£) (A. 5) 

where 

. 1 1 + sjx . . 

^(x) = - log - - (A.6) 
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Again, £ can be promoted to a continuous positive semi-definite real variable of which 
Pt(x) is a strictly increasing function. The first return probability generating function for 
the comb with teeth of length £ + 1 equally spaced at intervals of n is given by 

Pt*n(x) = 3 ~ f (x) - \ [(3 - P t {x) - 2P< n {x)f - 4G<g>{x; nf] k . (A.7) 

and P* n (x) is obtained by setting £ = oo in this formula. Pe t * n (x) is a strictly decreasing 
function of £ and increasing function of n, viewed as continuous positive semi-definite real 
variables. 

We also need the scaling limits of some of these quantities. They are 

lim a~^G^(a^; a~^v) = ^cosech(zz^) (A. 8) 

and 



lima 2 (l — p i x (af)] = — -f 2 tanhfpf 2) 

4 + 4tanhp^ co thz/^ + tanh 2 p^] * ■ (A. 9) 



2^ 



Bump functions 



A function tp : R — )■ R is a bump function if ip is smooth and has compact support. We 
will now prove some properties concerning the Mellin transformation of a bump function 
ip, A^[V ; ](s), which has support on [a, b] where b > a > 0. 

Lemma 9 The critical strip of the Mellin transform of the nth derivative of ip, ip( n \ is 
C for all n. 

Proof Recall that the Mellin transform is defined by, M.[tp](s) = J °° tp(x)x s ~ 1 dx. Since 
ip has compact support we have, 

f b 

M{4> (n) ](s) = / ^ n \x)x s ~ 1 dx (A.10) 

J a 

and since ip is smooth, \ip^\ is bounded on [a,b] by some constant K, so, 

r b 

\M[i) {n) ]{s)\ < K x s - l dx (A.ll) 

J a 

and the RHS is finite for all s e C since b > a > 0. This also shows that M[ip^](s) is 
holomorphic for all s. 



A. 3 Asymptotic Series and Dirichlet Series 
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Lemma 10 Given n e Z+, \M[ip](a + iD)\ < -^M[\ip {n) \]{a + n) for alls e C. 

Proof Recall from the previous lemma that the critical strip of the Mellin transform of 
if) and its derivatives coincides with C. We can therefore use the integral representation 
of the Mellin transform to prove statements valid for all sgC. By integration by parts, 

M[ip (n) ){s) = - /"V n+1) (x)— dx = --M[^ (n+1) ](s + 1) (A.12) 

J a S S 

and therefore, 

M[4>){8) = / ^ n \x)x s+n dx = w _! M[^](s + n). (A.13) 
Uk=o( s + k )^ Uk=o( s + k ) 

Hence, 

\M[1>](a + iT)\ <^j'\^ n \x)\x^ n dx = ^MU^\](a + n). (A.14) 

Given a bump function \l/± e which is always positive, has support [1, 1 ± e] and is scaled 
such that its integral is one, we define the cut-off function to be, 

T)±(x) = l- / ^± e (x)dx. (A. 15) 



Lemma 11 The critical strip of A4[r)±](s) is given be Re[s] > 0. The analytic continua- 
tion of JA[r]±\(s) to all s is given by 

M[ri ± ](s) = ±M[* ±e ](s + l). (A.16) 

Proof The analytic continuation of A^[r/±](s) may be obtained by applying integration 
by parts to the Mellin transform of r)± and recalling by lemma fl9]) that .M [^i^Ks) is 
holomorphic everywhere. 

A. 3 | Asymptotic Series and Dirichlet Series 

Starting with 

oo oo 

s±(y) = j>(*M4/) = mo) + X>(*fot('v) = tio) + s±(y) (A.i7) 

1=0 1=1 
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where r]± is the smooth cut-off function introduced in IA.2l and y controls where the cut-off 
occurs we take the Mellin transform of S±\y) with respect to y, 



S^^y^dy 
= V,(s)M[ V± ](s) 



V^s)M[^ ±e ){s + l)/s (A.18) 



where is the Dirichlet series associated to the measure, 

w=Ef' (A - 19) 

It is easy to see that the fundamental strip of S±\y) is Re[s] > 0, due to the compact 
support of 7]± on the positive real axis and so the Mellin transform does indeed exist. We 
may now invert the Mellin transform to obtain, 

S±(y) = MO) + ^ r +l0 °y- s V,(s) M[ ^ ]{s + 1) ds. (A.20) 

This may be computed by rewriting the above as, 

S±(y) = A*(0) + ^-.<f y- s V,(s) M[ ^ ±e]is + 1) ds - R(y) (A.21) 

where the contour C is the rectangle composed of the points {c — ioo, c + zoo, — N + 
ioo, —N — ioo} with N > 0, c such the contour is the right of all poles and 

m=r""»-VM M ™' + 1) *. (A.22) 

J—N—ioo s 

If has slow growth in the strip — N < Re[z] < then due to lemma fflOl) . which shows 
that A^[^± e ] decays faster than any polynomial as t goes to infinity, the contributions 
from integrating along the contours c + ioo to — N + ioo and from c — ioo to — N — ioo 
are zero. Furthermore the remainder term R(y) satisfies, 

N r-N+ioo 

\R(y)\ < V / \V,(s)\\M[^ ±t ](s + l)\ds (A.23) 

and so will only contribute terms of order y to S±(y). We therefore have, 

s±{ V ) = m + ^£y' sv ^ M[xSl± f s + 1) ds ( A - 24 ) 

R(y) (A.25) 



1 + res 



A. 3 Asymptotic Series and Dirichlet Series 
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where S is the set of positions of the poles of X> M and we have used the fact that T>^(0) = 
1 — /x(0). Finally, define x±( u ) = Y^alo f Jj (^)v±(^ / u ) ■ By relating this function to S±(y) 
we may write, 



res 



P|t(a) M^( a + l) u .. a = a< 



(A.26) 
(A.27) 
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APPENDIX B 



Liouville Cylinder Amplitudes 



The full cylinder amplitudes were given in a usable form recently by [36] . We have obtained 
similar results independently and reproduce them here for the purpose of completeness. 
The cylinder amplitude between a (k, I) and (r, s) brane, Z(k, I; cri\r, s; 02), is computed 
by the integral 

POO 

Z(k,l;a 1 \r,s;a 2 ) = dTZ ghost (T)Zn ouv m e (T)Z mattei (r), (B.l) 
Jo 

where r is the modular parameter of the cylinder in the closed string channel. The 
partition function for the ghosts is well known 

^ghost(r) = r/(q) 2 , (B.2) 

where q = exp(— 2ttt) and rj(q) denotes the Dedekind //-function. The Liouville and 
matter contributions are given by 

7 ST^ S(k,l),(a,b)S( r ,s),(a,b) , v / R „\ 

^matter = > ~ X(o,6)W) \P-&) 

poo 

liouville = / dP^ ai {P)^ a2 (P) XP (cO, (B.4) 

Jo 

where E QtP is the Kac table for the minimal model, S is the relevant modular S-matrix, 
\I/(cr) is defined in (15.441) and X(o,6)(q) an d Xp(q) are the characters for the minimal 
model primary field with Kac index (a, b) and the non-degenerate Liouville primary field 
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respectively. The characters are given by, 

1 oo 

T^ o6 )(q) = — — ^ ^2pqn+qa-pb) 2 / ipq _ ^(2pqn+qa+pb) 2 /4pq^ 

xp(q) 



,P 2 



r)(q)' 



(B.5) 
(B.6) 



We first concentrate on the ^matter- Recalling that the modular S-matrix for the (p,q) 
minimal model is 



S(k,l),(r,s) = 2W (-1) 

V pq 

the explicit expression for Z mautter is, 

■ P-1 9-1 



1+lr+ks 



qkr . pis, 



sin 7T ) sm(7r- 

p 



q 



(B.7) 



Z 



matter 



~ 2 y zZSZ |(-i) a(z+r+1)+b(fc+s+1) 



a=l b=l 



sin(Tr^) sin(Tr^) sin(Tr^) sin(TT^) 
sin^Wyr^) 



^(q) , 



t (2pqn+qa—pb) 2 / ipq q(2pqn+ija+pt>) 2 /4p<j 



(B.8) 



where we have used the symmetry of the Kac table to rewrite the limits on the a and b 
summation. Also, note that the quantity in the square brackets is symmetric in a and b 
and zero if a = or b = 0. This allows us to rewrite it as, 



•Tnatter 



p-1 q-1 ^ 

^ srn ( — )sin( — )x 

<2pq L p a 

a=l L i 



b=-(q-l) 



[B.9) 



p p o o niq) ' 

v ' n=— oo 

where t = ga + pfe. Substituting these expressions into flB.2j) gives, 
1 



o 



°° dP 47 ^ cos ( 27rcTlP ) cos(27tct 2 P) 
sinh(2vr6P) sinh(^) 



/2pq 

[ s i n ( — ) sin( — )£4_i(cos — )C/ r _x(cos — )C//_i(cos — )C/ s _i(cos — ' 
^— ' L v a v v a a 



= 1 b=-{q-l) 



dr q 



(2pgn+t) 2 /4pg+P 2 



(B.10) 
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Focusing our attention on the integral over the moduli, we compute, 



dr ^ q 



(2pqn+t) 2 /Apq+P 2 



1 oo 

— y 

2tt ^ 



(2pqn + t) 2 /Apq + P 2 



sinh( 



2ttP_\ 

■pl J 



2JpqP cos 2_ cosh 



which gives, 



Z 



2tt 2 f°° dP cos(2tt(TiP) cos(27ror 2 P) sinh(fg, 



V%pq Jo P 
where F h ^ s {z) is given by, 



sinh(27r6P) sinh(^; 



#k,l,rA , — ), 



X p (k,r) X q (l,s) p—1 q—1 

Fk4,rM= E E E E 

J7=|fc-r|+l,2p=|j-s|+l,2 a=l b=-(q-l) 



. 7Ct . 7Tt 

sin — J sin — J x 
V Q 

^_i(cos^)f/ p _i(cos— ) _, 



(B.11) 



(B.12) 



q cos — — cos z 

pq 



(B.13) 



where \ p (k, r) = mm(k + r — 1, 2p — 1 — k — r). To compute (IB. 121) we first note that 
the integrand is symmetric in P and so we can extend the integral to the whole real line. 
Without loss of generality we assume o\ > a 2 , and break up the factor of cos(27rcr 1 P) in 
(IB.12p into its exponentials. For the integral containing, exp(27rza 1 P) we can close the 
integral along the real line by a semi-circle to produce the contour C + in the upper-half 
plane, for the other term we use a contour, C_ in the lower-half plane. However, the 
integral using the contour in the lower-half plane, by a change of variables from P to — P 
is equal to the integral using the contour in the upper half-plane, we therefore get, 



Z(k,l;a 1 \r, s;cr 2 ) = A c <h dP 

Jc+ 



G(P,ai,a 2 )F k>l>r , 



,2niP^ 



pq 



(B.14) 



where Ac is a numerical constant and 

i e 2 7 ri<T 1 p cos ( 27r(T2 p) sinh(g / 



G{P,a u a 2 ) 



P sinh(27r6P) sinh(^) 



(B.15) 
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The poles of F) 



k,L: 



'2vriP> 



/PQ 

The residue of F k ^ s {^) at P e S F is, 



occur at P G Sp = {?J^ '■ n E Z,n mod p, n 7^ mod q}. 



2mP in 
res {F k ^ s {—— ; P = = 



„VM . r' Kn \ ■ / 7rn \TT I nn \TT I 

-2- — sin — J sin — Jt4_i(cos — )t/ r _i(cos — J x 
2m p q V V 



TT t nn \TT t nU \ 1 

C/j_i(cos — )C/ s _i(cos — )- 



P9 



(B.16) 



The residues of the G(P) are, 

27rires G{P);P = ie 
2irires \g(P);P 
27rires \g(P);P 
2 tt /res G(P);P 



1 x 1 

— 7T<7i 



'2e 



(B.17) 



2^/pg 
inq 

2y/pql n 
inpq 1 2 



2 hp^i sin(7rn/g) 

= -(-l) n e Vp? cos h(n7r6(T2)^V — tt, 
n sm(irpn/q) 



1 =-(-!)" 
J n 

1 = — (-l) n 
2^J npg 



zigzm , , ... sin(7rn/») 
e cosh(n7r<j2/o)- 



sin(7rgn/p) ' 



(p+9+i) e cosh(n7rg6o- 2 )- 



Finally, a useful identity for computing the sums found in the computation of the cylinder 
amplitudes is, 



Eh 



t £ -at 



t=l 



cosh f3t cos <j)t = (B.18) 
2a — ln(2(cosh 2 a + cosh 2 (3 — 27 cosh a cosh (3 cos — sin 2 0)). 



APPENDIX C 



Loop equations for Disc with Handle 



To complete the calculation of the disc with handle amplitude in the SX matrix model we 
need a number of additional loop equations. In order to calculate the quantities Wsxsx 
and Wsxsxsx we need, 

x ^ A - + £ (_l_ A - ( . Y , S) _2_ + _2_ /f » WS) _l_) (o 

s - s + f (^ A ' (A ' S) ^ + ^ A " ,(XS) ^) (a2) 
x + + (a3) 

Finally the following loop equations allow for quantities of the form, Wsx-,sx etc. to be 
found, 

X -+ X + e (rV(A', S)-^ + j^KHX, S)-^) TH7( S , X) (C.4) 
S ^S + t {^K(X, S)-^ + -L^KHX, S)-^) M (S. X) (C.5) 
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